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APPEKDIZ  B 

SIKORSKY  CH-S3B  HELICOPTER  XND  ROTORVXSB  DATA 


Sikorsky  CH-53E  helicopter  characteristics  used  with  the  ROTWASH 
analysis  program  are  listed  in  table  B-1  as  obtained  from 
reference  B-1.  Figure  B-1  presents  a  three-view  drawing  of  the 
helicopter.  Table  B-2  provides  a  summary  tabulation  of 
parameters  defining  the  conditions  for  all  rotorwash  data 
presented  in  this  appendix.  Flight  test  data  used  in  correlation 
with  the  ROTWASH  calculated  data  are  obtained  from  reference  2 
and  were  measured  along  the  270-degree  azimuth  (out  the  left  side 
of  the  helicopter).  Distance  from  rotor  center  (DFRC),  gross 
weight  (GW) ,  and  rotor  height  above  ground  level  (RHAGL)  are  the 
primary  independent  variables  for  these  measured  data. 

Discussion  of  results  is  presented  in  section  3.1  of  this  report. 


REFERENCES 

B-1.  Prouty,  R.W.,  Helicopter  Performance.  Stability,  and 

Concroi .  Robert  E.  Krieger  Publishing  Company,  Malabar,  FL 
32950,  1990. 

B-2.  Harris,  D.J.,  and  R.D.  Simpson,  '*CH-53E  Helicopter  Downwash 
Evaluation,"  Naval  Air  Test  Center  Technical  Report  No.  SY- 
89R-78,  August  1,  1978. 


TABLE  B-1  CH-53E  ROTWASH  INPUT  DATA 

Parameter  Value 


Rotor  radius,  feet  39.5 

Distance  between  rotor  centers,  feet  0.0 

Airframe  download,  percent  of  rotor  thrust  5.0 

Distance  from  wheels  to  rotor  plane,  feet  17.0 

Rotor  speed,  RPM  177.0 

Rotor  tip  speed,  feet/second  732.0 

Number  of  rotor  blades  per  rotor  7 

Density  ratio  1.0 
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FIGURE  B-1  THREE-VIEW  DRAWING  OF  THE  SIKORSKY  CH-53E 


TABLB  B-2  BVALDATIOM  MATRIX  TOR  CB~S3B  TLiaBT  TBST 
NATBBMATICAL  MODEL  DATA  CORRBLATIOM 


FIGURE 

NUMBER 

GROSS 

WEIGHT 

(1^ 

DISK 

LOADING 

(Ibs/ft^) 

ROTOR 
HEIGHT 
( feet) 

DISTANCE  FROM  ROTOR  CENTER 
(DFRC)  (feet) 

B-2 

70,000 

14.28 

37.0 

31.6,  39.5,  49.4,  59.3, 

69.1,  79.0,  118.5,  177.8 

B-3 

70,000 

14.28 

77.0 

31.6,  39.5,  49.4,  59.3, 

69.1,  79.0,  118.5,  177.9 

B-4 

70,000 

14.28 

117.0 

31  .,  39.5,  49.4,  59.3,  | 

59.1,  79.0,  118.5,  177.8  E 

B-5 

56,000 

11.42 

37.0 

31.6,  39.5,  49.4,  59.3,  1 

69.1,  79.0,  118.5,  177.8  I 

B-6 

56,000 

11.42 

77.0 

31.6,  39.5,  49.4,  59.3,  I 

69.1,  79.0,  118.5,  177.8 

B-7 

56,000 

11.42 

117-0 

31.6,  39.5,  49.4,  59.3,  j 

69.1,  79.0,  118.5,  177.8 

B-8 

45,000 

9.18 

37.0 

31.6,  39.5,  49.4,  59.3, 

69.1,  79.0,  118.5,  177.8 

B-9 

45,000 

9.18 

77.0 

31.6,  39.5,  49.4,  59.3, 

69.1,  79.0,  118.5,  177.8 

B-10 

1  1  r  -  r  1  1 

45,000 

9,18 

117.0 

31.6,  39.5,  49.4,  59.3, 

69.1,  79.0,  118.5,  177.8 

NOTES: 


1)  The  values  of  DFRC  are  only  applicable  along  the  270-degree 
a::irnuth. 

2)  Ambient  winds  varied  between  0  and  3.5  knots. 


3)  Atmospheric  density  ratio  was  assumed  equal  to  1.0  since 
pressure  altitude  (which  was  near  sea  level)  was  not 
documented  in  reference  B-2.  Ambient  temperature  was 
measured  from  39  to  45  degrees  Fahrenheit  during  testing. 
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RIDL  =  5.0  PCT 
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FIGURE  B-2  CH-53E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEGRBE  A2IMDTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AND  A  GROSS  HEIGHT  OF  70,000  POUNDS  (continued) 
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1  SIKORSKY  CH.53E  HELICOPTER 

GW  =  70000  LB 

DL  =  14.28  PSF 

SIGP  (o')  =  1.0 
A2IMUTH  =  270  DEG 

DFRC  s  118,5  FT 
RHAGL  =  37.0  FT 
WHAGL  =  20.0  FT 
WIND  <  3.£  KTS 

RADIUS  s  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 

MEAN  FLIGHT  TEST  DATA 
PEAK  FUGKT  TEST  DATA 
MEAN  CALCULATED  DATA 
PEAK  CALCULATED  DATA 


SIKORSKY  CH-53E  HELICOPTER 


GW  =  70000  LB  DFRC  =  177,9  FT  RADIUS  =  39.5  FT 

DL  =  14,28  PSF  RHAGL  =  37.0  FT  ROTSEP  =  0.0  FT 

SIGP(o')  =  1,0  WHAGL  =  20.0  FT  RIDL  =  5.0  PCT 

AZIMUTH  =  270  DEG  WIND  <  3.5  KTS 
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•  PEAK  FUGHT  TEST  DATA 

-  MEAN  CALCULATED  DATA 

-  PEAK  CALCULATED  DATA 
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FIGURE  B-2  CE*>53E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEGREE  AZIKUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AND  A  GROSS  WEIGHT  OF  70^000  POUNDS  (continued) 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  70000  LB 
DL  =  14.28  PSF 
SlGP(o’)  =  1.0 
AZIMUTH  =  270  DEG 


DFRC  s  49.4  FT 
RHAGL  =  77.0  FT 
WHAGL  =  60.0  FT 
WIND  <  3.5  KTS 


RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 
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PIGURB  B-3  CH-53E  MEAN  AND  PEAK  VELOCITY  PROPILE8  FOR  BIGHT 
270-DEORBE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OP  77  FEET  AND  A  GROSS  WEIGHT  OF  70,000  POUNDS  (Continued) 
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SIKORSKY  CH-53E  HELICOPTER  I 

GW  =  70000  LB 

DL  =  14,28  PSF 

SIGP  (o’)  =  1.0 
AZIMUTH  =  270  DEG 

DFRC  =  69.1  FT 
RHAGL  =  77.0  FT 
WHAGL  =  60,0  FT 
WIND  <  3.5  KTS 

RADIUS  =  39.5  FT 
ROTSEP  =  0,0  FT 
RIDL  =  5.0  PCT 
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•  PEAK  FLIGHT  TEST  DATA 

-  MEAN  CALCULATED  DATA 

- PEAK  CALCULATED  DATA 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  70000  LB 
DL  =  14  28  PSF 
SIGP  (a)  =  1.0 
AZIMUTH  =  270  DEG 


DFRC  =  79.0  FT 
RHAGL  =  77.0  FT 
WHAGL  =  60.0  FT 
WIND  <  3,5  KTS 


RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 
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FIOURE  B-a  CH-53E  MEAN  AMD  PEAK  VELOCITY  PROFILES  FOR  BIGHT 
270-DEOREE  AZIMUTH  RADIAL  STATIONS  AT  A  KOTOR  HEIGHT 
OF  77  FEET  AND  A  GROSS  HEIGHT  OF  70,000  POUNDS  (continued) 


HEIGHT  ABOVE  GROUND  LEVEL.  FEET  HEIGHT  ABOVE  GROUND  LEVEL.  FEET 


1  SIKORSKY  CH-53E  HELICOPTER 

GW  =  70000  LB 

DL  =  14,28  PSF 

SIGP  (o’)  =  1,0 
AZIMUTH  =  270  DEG 

DFRC  = 
RHAGL 
WHAGL 
WIND  < 

118.5  FT 
=  77.0  FT 
=  60.0  FT 

3.5  KTS 

RADIUS  =  39.5  FT 
ROTSEP  =  0,0  FT 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  70000  LB 
DL  =  14,28  PSF 
SIGP  (o')  =  1.0 
AZIMUTH  =  270  DEG 


DFRC  =  177.9  FT 
RHAGL  =  77,0  FT 
WHAGL  =  60.0  FT 
WIND  <  3,5  KTS 


RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  !-•  5.0  PCT 
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FI< URE  B-3  CH-53B  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  BIGHT 
270-DEaRBE  AZINDTB  RADIAL  STATIONS  XT  A  ROTOR  HBIGHT 
OF  77  FEET  AND  A  GROSS  WEIGHT  OF  70,000  FOUNDS  (continued) 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  70000  LB 
DL  =  14.28  PSF 
SIGP  (o')  =  1.0 
AZIMUTH  =  270  DEG 


DFRC  =  59.3  FT  RADIUS  =  39.5  FT 
RHAGL  =  117.0  FT  ROTSEP  =  0.0  FT 
WHAGL  =  100.0  FT  RIDL  =  5.0  PCT 
WIND  <  3.5  KTS 
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FIGURE  B-4  CH-53B  MEAN  AMD  PEAK  VELOCITx'  PROFILES  FOR  EIGHT 
270-DEGREE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  117  FEET  AMD  A  GROSS  WEIGHT  OF  70,000  POUNDS  (contl&ued) 
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1  SIKORSKY  CH-53E  HELICOPTER 

GW  =  70000  LB 

DL  =  14.28  PSF 
SIGP(o’)  =  1.0 
AZIMUTH  =  270  DEG 

DFRC  a  69.J  FT 
RHAGL  =  117.0  FT 
WHAGL  =  100.0  FT 
WIND  <3.5  KTS 

RADIUS  =  39.5  FT 
ROTSEP  =  0,0  FT 
RIDL  =  5.0  PCT 
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PEAK  FLIGHT  TEST  DATA 
MEAN  CALCULATED  DATA 
PEAK  CALCULATED  DATA 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  70000  LB  DFRC  =  79.0  FT  RADIUS  =  39.5  FT 

Dl  =  14.28  PSF  RHAGL  =  117.0  FT  ROTSEP  =  0.0  FT 

SIGP(a’)  =  1.0  WHAGL  =  100.0  FT  RIDL  =  5,0  PCT 

AZIMUTH  =  270  DEG  WIND  <  3.5  KTS 
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•  PEAK  PLIGHT  TEST  DATA 

-  MEAN  CALCULATED  DATA 

- PEAK  CALCULATED  DATA 
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FIGURE  B-4  CH-53E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270’-DEGRBB  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  117  FEET  AND  A  CROPS  WEIGHT  OF  70,000  POUNDS  (continued) 
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FIGURE  B-4  CH-53E  MEAN  AMD  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEOREE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  117  FEET  AND  A  GROSS  WEIGHT  OF  70,000  POUNDS  (continued^ 
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SIKORSKY  CH-G3E  HELICOPTER  I 

GW  =  56000  LB 

OL  =  11.42  PSF 
SIGP(o')  s  1.0 
AZIMUTH  =  270  DEG 

DFRC  ®  31.6  FT 
RHAGL  =  37.0  FT 
WHAGL  »  20.0  FT 
WIND  <  3.5  KTS 

RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 
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PROFILE  VELOCITY,  KNOTS 

FIOURB  B-5  CH-53E  HEAM  AND  PEAK  VELOCITY  PROPILES  FOR  EIGHT 
270-DEOIlEE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AND  A  GROSS  WEIGHT  OF  56/000  POUNDS 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  56000  LB 
DL  =  11.42  PSF 
SIGP  (o’)  =  1.0 


DFRC  =  59.3  FT 
RHAGL  =  37.0  FT 
WHAGL  =  20,0  FT 
WIND  <  3.5  KTS 


RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5,0  POT 
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FIGURE  B-5  CH-53E  MEAM  AND  PEAK  VELOCITY  PROFILES  FOR  BIGHT 
270-DE0REE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AND  A  GROSS  WEIGHT  OF  56,000  POUNDS  (continued) 
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SIKORSKY  CH-53e  HELICOPTER 


GW  s  56000  LB  DFRC  =  69.1  FT  RADIUS  s  39.5  FT 

DL  =  11.42  PSF  RHAGL  =  37.0  FT  ROTSEP  =  0.0  FT 

SlGP(a')  =  1.0  WHAGL  =  20.0  FT  RIDL  =  5.0  POT 

AZIMUTH  =  270  DEG  WIND  <  3,5  KTS 


MELAN  FLIGHT  TEST  DATA 
PEAK  FUGHT  TEST  DATA 
MEAN  CALCULATED  DATA 
PEAK  CALCULATED  DATA 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  56000  LB 
DL  =  11.42  PSF 
SIGP(o')  =  1.0 
AZIMUTH  =  270  DEG 


DFRC  -  79.0  FT 
RHAGL  =  37.0  FT 
WHAGL  =  20.0  FT 
WIND  <  3.5  KTS 


RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 
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-  peak  CALCULATED  DATA 


10  20  30  40  50  60  70  80 

PROFILE  VELOCITY,  KNOTS 


90  100 


FIOURZ  B-5  CH-S3E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEaREB  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AMD  A  GROSS  WEIGHT  OF  56,000  POUNDS  (continued) 
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SIKORSKY  CH-53E  HELICOPTER 


GW  =  56000  LB 
DL  =  11.42  PSF 
SIGP  (o')  =  1.0 
AZIMUTH  =  270 


DEG 


DFRC  =  177.9  FT 
RHAGL  =  37.0  FT 
WHAGL  =  20.0  FT 
WIND  <  3.5  KTS 


RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 
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FIOURB  B-5  CH-53E  MEAM  AVID  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEaREB  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AMD  A  GROSS  WEIGHT  OF  56,000  POUNDS  (continued) 
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I  SIKORSKY  CH-53E  HELICOPTER 

GW  =  56000  LB 

DL  =  11.42  PSF 

SIGP  (o’)  =  1.0 
A2IMUTH  270  DEG 

DFRC  =  31.6  FT 
RHAGL  =  77.0  FT 
WHAGL  =  60.0  FT 
WIND  <3.5  KTS 

RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 
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SIKORSKY  OH-53E  HELICOPTER 

GW  =  56000  LB 

DL  =  11.42  PSF 

SIGP  (o’)  =  1.0 
AZIMUTH  =  270  DEG 

DFRC  s 
RHAGL 
WHAGL 
WIND  < 

39J  FT 
=  77.0  FT 
=  60.0  FT 

3.5  KTS 

RADIUS  =  39.5  FT 
ROTSEP  =  0.0  FT 
RIDL  =  5.0  PCT 
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FIGDRE  B-6  CH-S3E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  BIGHT 
270-DEGREE  A2INDTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  77  FEET  AMD  A  GROSS  HEIGHT  OF  56,000  POUNDS  (COntlnUAd) 
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1  SIKORSKY  CH-53E  HELICOPTER  I 

GW  =  56000  LB 

DL  «  11.42  PSF 

SIGI  (o’)  =  1.0 
AZIMUTH  =  270  DEG 

DFRC  =  49,4  FT 
RHAGL  =  77.0  FT 
WHAGL  ==  60.0  FT 
WIND  <  3.5  KTS 

RADIUS  =  39.5  FT 
BOTSEP  =  0.0  FT 
RIDL  =  5,0  POT 
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o  MEAN  FLIGHT  TEST  DATA 
•  PEAK  FUGHT  TEST  DATA 
— —  MEAN  CALCULATED  DATA 
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FIGURE  B-6  CH>-53E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-'DEGRBE  AZIHUTH  RADIAL  STATIONS  AT  A  ROTOR  HFIGHT 
OF  77  FEET  AND  A  GROSS  WEIGHT  OF  56,000  POUNDS  (continuad) 
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FIGURE  B-7  CH-53E  NEAR  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
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riaURE  B>7  CH-53E  MEAK  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEGREE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AND  A  GROSS  WEIGHT  OF  45,000  POUNDS  (continued) 
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FIGURE  B-7  CH-53E  MEAN  AND  FEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEGREE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  37  FEET  AND  A  GROSS  WEIGHT  OF  45,000  POUNDS  (continued) 
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FIGURE  B-8  CH-53B  MBAM  AND  PEAK  VELOCITY  PROFILES  FOR  BIGHT 
270-DE0REE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  LBIOHT 
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FIGURE  B-8  CH-53E  MEAN  AND  PEAK  VELOCITY  PR0FILE8  FOR  EIGHT 
270-DEGREE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  77  FEET  AMD  A  GROSS  WEIGHT  OF  45,000  POUNDS  (continued) 
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FIGURE  B-8  CH-53B  KEAN  300)  PEAIX  VELOCITY  PROFILES  FOR  EIGHT 
270-DEGREB  azimuth  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  77  FEET  AND  A  GROSS  HEIGHT  OF  45,000  POUNDS  (continued) 
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FIGURE  B-9  CH-53E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEGRBE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  77  FEET  AND  A  GROSS  REIGHT  OF  45,000  POUNDS 
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FIGURE  B-9  CH-53E  MEAN  AND  PEAK  VELOCITY  PROFILES  FOR  EIGHT 
270-DEGREE  AZIMUTH  RADIAL  STATIONS  AT  A  ROTOR  HEIGHT 
OF  117  FEET  AND  A  GROSS  WEIGHT  OF  45,000  POUNDS  (continued) 
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APPENDIX  D 
USER'S  QUIDS 


D.l  ROTWASH  PROGRAM  USER'S  GUIDE 

The  Style  of  the  ROTWASH  user's  guide  presented  in  this  appendix 
is  primarily  narrative.  This  format  is  designed  to  guide  the 
reader  through  a  step-by-step  explanation  on  the  use  of  each 
program  software  option.  Example  output  from  each  option  is 
presented  exactly  as  it  would  be  viewed  by  the  user  on  a  video 
terminal  for  reference  purposes.  Examples  of  user  keyboard  input 
are  presented  in  <BOLD>  text  as  an  aid  to  the  reader. 

D . 1 . 1  GETTING  STARTED 

The  ROTWASH  program  is  executed  by  typing  the  program  name  at  the 
DOS  system  prompt: 


<ROTWA8H> 

To  avoid  possible  system  errors,  the  user  should  execute  the 
program  from  the  directory  containing  the  program  (or  set  the 
appropriate  DOS  system  PATH  command) .  The  user  should  also  be 
aware  that  menu  and  data  printouts  to  the  screen  will  not  work 
correctly  unless  the  dsvlce=C: \D08\AHSI.8Y8  command  is  contained 
in  the  CONFIG.SYS  file. 

The  ROTWASH  program  responds  with  the  screen  output  presented  in 
figure  D-1.  This  output  is  the  ROTWASH  program  header  page. 


ROTWASH  PROGRAM 

ROTORCRAFT  DOWNWASH  HAZARD  ANALYSIS 

EMA  /  SYSTEMS  CONTROL  TECHNOLOGY 
***  PROGRAM  VERSION  2.1,  APRIL  1993  *** 


PRESS  <RETURN> 


FIGDRB  D-1  ROTWASH  PROGRAM  HEADER  OUTPUT 


After  typing  the  carriage  return  <RBTOitll>  or  •<E1ITER>  key,  the 
user  is  asked  to  specify  the  path  for  ROTWASH  program  input  and 
output  (I/O).  The  only  option  for  input  at  the  present  time  is 
the  terminal  keyboard.  This  is  specified  by  typing  <C01I>  for 
console  or  keyboard  (lover  case  inputs  such  as  <oon>  are  also 
permitted) .  Program  output  may  be  sent  to  one  of  three  different 
locations.  These  output  locations  are: 


Screen 

Printer 

Disk  Plotting  File 


<COM> 

<PRN> 

<PLT> 


An  example  of  the  video  screen  output  for  this  program  menu  is 
presented  in  figure  D-2  where  the  <CON>  option  has  been  chosen 
for  both  input  and  output. 


I/O  CAN  BE  DIRECTED  TO  FILES  OR  DEVICES 

VALID  DEVICES  AJRE  AS  FOLLOWS: 

<CON>  ==>  CONSOLE 
<PRN>  ==>  PRINTER 
<PLT>  ==>  GRAPHICS  FILE 


ENTER  INPUT  FILE/DEV  NAME  ===>  CON 
ENTER  OUTPUT  FILE/ DEV  NAME  ==>  CON 


FIGURE  D>2  ROTWASH  PROGRAM  INPOT/OUTPUT  CONTROL  MENU 


D.1.2  INPUT  DATA  REQUIREMENTS 

Rotorcraft  characteristics  and  atmospheric  conditions  that  are 
common  to  all  program  options  are  input  to  the  program  using  the 
master  input  data  menu.  Four  basic  configurations  of  rotor  or 
propeller  driven  aircraft  can  be  represented  using  this  menu. 
These  configurations  include  single  and  tandem  rotor  helicopters, 
tiltrotors,  and  twin-propeller  tiltwings.  This  menu  is  presented 
to  the  user  as  shown  in  figure  D-3  after  the  I/O  menu  is 
completed.  The  default  values  provided  in  the  menu  define  the 
Bell  XV-15  tiltrotor.  Design  data  describing  most  other  modern 
types  of  rotorcraft  are  provided  in  appendix  A. 


ROTWASH  USER  INPUT  DATA 

MENU 

CODE 

PARAMETER 

VALUE 

UNITS 

A 

NUMBER  OF  ROTORS  (1  OR  2) 

2 

-ND- 

B 

HUB  TO  HUB  ROTOR  SEPARATION 

32.2 

FT 

C 

ROTOR  RADIUS 

12.5 

FT 

D 

GROSS  WEIGHT 

13000.0 

LB 

E 

FUSELAGE  DOWNLOAD  FACTOR 

13.0 

PCT 

r 

ROTOR  HEIGHT  ABOVE  GROUND 

37. C 

FT 

G 

SHAFT  TILT  ANGLE  ^<20  DEG) 

.0 

DEG 

H 

AIR  DENSITY  RATIO 

1.0000 

ND 

1 

AMBIENT  WIND  (-10  TO  10  KT) 

.0 

KT 

ENTER 

DATA  ENTRY  CODE  OR  <RETURN>  TO 

CONTINUE  D 

GROSS  WEIGHT  «  13000.0 

ENTER 

NEW  VALUE  OR  <RETURN>  TO  LEAVE 

AS 

IS  13500.0 

rZGDSB  D-3  NXBTER  IHPDT  DATA  MENU 

The  first  input  data  variable  listed  on  the  menu  is  the  number  of 
rotors  (or  propellers) .  This  value  will  always  be  two  except  for 
single  main  rotor  helicopters  (the  tail  rotor  is  not  considered  a 
lifting  rotor  by  the  ROTWASH  analysis) .  The  next  four  menu 
variables  define  the  design  details  of  the  rotorcraft 
configuration.  These  variables  are  the  hub-to-hub  separation 
distance  for  twin-rotoi  configurations  (in  feet) ,  rotor  radius 
(feet) ,  rotorcraft  gross  weight  (pounds) ,  and  the  rotor-on-wing 
download  factor  (percent) .  The  download  factor  represents  the 
percent  increase  in  hover  rotor  thrust  required  to  overcome  the 
rotor-induced  vertical  drag  force  on  the  airframe  (where  thrust 
is  initially  assumed  to  equal  gross  weight) .  For  most 
helicopters,  this  value  is  less  than  5  percent.  Download  on 
tiltrotors  can  generally  be  expected  to  vary  from  8  to  13  percent 
because  of  the  large  wing  area  under  the  rotor. 

The  next  two  input  data  variables  define  the  position  of  the 
rotor  with  respect  to  the  ground  plane.  The  rotor  height  (feet) 
is  defined  as  the  distance  from  the  ground  surface  to  the  plane 
of  the  rotor.  The  mast  angle  (degrees)  is  defined  as  the  tilt  of 
the  rotor  plane  with  respect  to  the  ground  plane.  The  mast  angle 
is  defined  as  0  degrees  when  the  plane  of  the  rotor  is  parallel 
to  the  ground  plane  and  a  positive  angle  is  a  forward  tilt  of  the 
mast  (which  is  presently  limited  in  the  program  to  20  degrees) . 
Since  most  rotorcraft  hover  with  the  plane  of  the  rotor  parallel 
to  the  ground,  it  is  recommended  that  caution  be  exercised  when 
non-zero  values  are  used  for  mast  angle.  Non-zero  uses  of  the 
variable  might  involve  hover  investigations  for  tiltwing  aircraft 
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with  the  wing  tilted  forward  (where  a  pitch  fan  provides  the  tri&i 
pitching  aoment) .  Alsoi  United  rotorwash  investigations  for 
takeoff  and  landing  naneuvers  night  be  attempted  if  approxinate 
rotor  thrust  levels  are  known  for  various  segments  of  the 
naneuvers . 

The  last  two  input  data  variables  on  the  menu  define  the 
atmospheric  conditions.  The  air  density  ratio  is  defined  as  the 
ratio  of  the  desired  air  density  to  the  sea  level  standard  air 
density  (0.0023769  slugs/feet^) .  Ambient  wind  speed  (knots)  is 
specified  by  the  user  up  to  a  limit  of  10  knots.  Values  greater 
than  this  limit  are  believed  to  invalidate  several  empirically 
determined  mathematical  modeling  assumptions  (limits  are 
discussed  in  sections  2  and  3  of  Volume  I  of  this  report) . 

The  mechanics  of  using  the  menu  are  quite  simple.  The  user  types 
in  the  code  value  for  the  variable  to  be  changed  and  then  types 
<RETURN>.  The  next  prompt  asks  for  the  new  parameter  value. 

After  this  value  and  another  <RETniUI>  are  typed,  the  menu  is 
rewritten  to  the  screen  with  the  new  value.  This  simple  process 
is  continued  until  the  user  specifies  each  variable  to  Its 
desired  value.  At  this  point,  the  <RBT0RN>  key  is  typed  by 
itself.  This  menu  can  also  be  reached  from  most  of  the  other 
menus  in  the  program  whenever  the  user  decides  that  the  basic 
configuration  needs  to  be  modified.  This  is  accomplished  by 
typing  <N>  for  NEW  CASE  when  the  option  is  offered. 

D.1.3  ANALYSIS  PROBLEM  DEFINITION 

After  the  master  input  data  menu  is  completed,  the  user  specifies 
the  desired  type  of  analysis  option.  Figure  D-4  presents  the 
program  logic/comment  menu  and  the  associated  list  of  default 
values.  This  menu  has  two  groupings  of  parameters  which  need  to 
be  specified. 

The  first  parameter  on  the  menu  specifies  the  choice  of  either 
the  velocity  calculation  analysis  option  or  the  hazard  analysis 
option.  The  velocity  analysis  option  is  the  default  option  on 
the  menu.  This  option  is  otherwise  chosen  by  typing  the  code 
<A>,  then  <RETURN>,  then  <V>,  and  finally  <RBTURII>.  The  same 
process  is  used  to  specify  the  hazard  analysis  option  except  that 
<H>  is  substituted  for  <v>  as  the  parameter  value.  The  second 
menu  parameter  provides  an  interactive  toggle  switch  for  the 
option  which  writes  out  data  files  to  disk  for  graphics  programs. 
(This  parameter  is  also  offered  as  an  option  on  the  initial 
ROTWASH  menu  by  typing  <PLT>) .  There  is  no  limit  to  the  number 
of  times  this  switch  can  be  toggled.  As  long  as  the  parameter 
has  a  <Y>  or  •'yes'*  value,  the  user  must  specify  output  filenames 
before  data  files  are  written  to  disk.  The  user  is  not  allowed 
to  write  over  files  previously  written  to  disk  by  specification 
of  the  same  filename  twice. 
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ROTWASH  PROGRAM 

LOGIC/ COMMENT  MENU 

CODE 

PARAMETER 

VALUE 

A 

ANALYSIS  TYPE, 

<V>  OR  <H> 

V 

B 

GRAPHICS  FILE, 

<Y>  OR  <N> 

N 

USER  INPUT  COMMENTS 

(FOR  "PRN"  AND 

••PLT''  OUTPUT) 

< - 

50  SPACES 

- > 

C 

XV-15  CHARACTERISTICS  ARE  USED  AS 

1  INPUT  DATA 

D 

GROSS  WEIGHT  MIGHT 

BE  ONE  OF  THE 

COMMENT  STRINGS 

ENTER 

CODE  FOR  DATA  INPUT  ' 

OR  <RETURN>  TO 

CONTINUE  >==> 
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The  last  two  lines  in  the  menu  are  used  to  specify  user  comments 
in  all  data  files  that  are  written  to  disk.  These  comments  are 
also  written  out  as  header  information  on  screen  output  sent 
directly  to  the  printer.  Both  of  the  comment  lines  can  be 
changed  at  any  time  during  program  execution.  The  only 
restriction  is  that  the  character  strings  on  both  lines  be  less 
than  or  equal  to  50  characters.  The  arrowhead  symbols  above  the 
comment  lines  in  the  menu  define  a  SO-space  line  width. 

If  the  velocity  analysis  option  is  specified,  the  user  is  then 
required  to  choose  one  of  the  four  analysis  options  presented  in 
figure  D-5.  Velocity  analysis  options  reached  through  this  menu 
are  the: 

1.  simple  wall  jet  (for  both  single  and  twin-rotor 
configurations) , 

2.  interaction  plane  (twin-rotor  only), 

3.  ground  vortex  (single  rotor  only),  and 

4.  disk  edge  vortex  (single  rotor  only) . 
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SELECT  TYPE  OF  FLOW  TO  BE  ESTIMATED 


WALL  JET  PROFILE, 

TYPE 

<W> 

INTERACTION  PLANE  PROFILE, 

TYPE 

<J> 

GROUND  VORTEX, 

TYPE 

<G> 

DISK  VORTEX, 

TYPE 

<D> 

TO  EXIT  PROGRAM, 

TYPE 

<X> 

ENTER  DATA  ENTRY  CODE 

II 

II 

V 

C 
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Each  of  these  four  options  is  described  with  flowfield 
characteristics  sketches  in  subsequent  sections.  For  technical 
details  and  discussion  on  practical  limitations  of  these  options, 
the  user  is  referred  to  Sections  2  and  3  of  Volume  I  of  this 
report.  The  choice  of  one  of  these  menu  options  is  made  by 
typing  the  appropriate  code  and  <RETURM>.  If  one  of  the  five 
allowable  characters  is  not  chosen,  the  menu  will  reappear  and 
the  user  will  be  forced  to  choose  an  acceptable  option. 

If  the  hazard  analysis  option  is  specified,  the  user  is  presented 
with  the  figure  D-6  menu.  This  menu  allows  the  user  to  choose 
either: 

1.  human  overturning  force/moment  analysis,  or 

2.  particulate  cloud  analysis. 

Both  of  these  analyses  can  be  applied  to  either  single  or  twin- 
rotor  configurations.  The  mechanics  of  this  menu  operate  exactly 
like  those  of  the  velocity  analysis  option  menu. 

D.1.4  THE  WALL  JET  OPTION 

Rotorwash  velocity  profiles  are  calculated  for  single  main  rotor 
helicopter  configurations  using  the  wall  jet  option.  Velocity 
profiles  along  the  0-  and  180-degree  azimuths  for  tandem 
helicopters  and  the  90-  and  270-degree  azimuths  for  twin-rotor 
side-by-side  configurations  are  also  calculated  using  this  option 
(90  degrees  is  out  the  right  wing  on  a  tiltrotor  and  0  degrees  is 
along  the  centerline  of  the  fuselage  for  tandem  rotor 
helicopters) .  Figure  D-7  provides  a  three-dimensional  view  of 
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the  rotorvash  flowflelds  associated  with  both  rotor 
configurations.  Figure  D-8  provides  a  cross-sectional  view  of 
the  nondimensionalized  ROTWASH  wall  jet  velocity  profile  model. 
The  program  menu  associated  with  the  use  of  this  option  is 
presented  in  figure  D-9. 

The  velocity  profile  status  menu  provides  the  user  the  option  to 
specify  four  parameters  before  proceeding  with  detailed 
calculations.  The  horizontal  distance  on  the  ground  from  the 
center  of  the  rotor  to  where  the  velocity  profile  should  be 
calculated  is  the  first  parameter  specified  on  the  menu.  Figure 
D-10  is  provided  to  illustrate  this  geometry  graphically  using 
the  Bell  XV-15|  which  is  the  more  complex  example  (flight  test 
data  results  associated  with  this  figure  are  documented  in 
reference  D-3) .  To  specify  the  profile  station  position  for  the 
wall  jet  with  the  XV-15  (270-degree  radial) ,  the  user  would 
measure  the  distance  from  the  aircraft  centerline  (DFAC)  and 
subtract  16.1  feet  (the  distance  from  the  centerline  to  the 
center  of  the  rotor) .  The  remaining  two  position-related 
parameters  to  be  specified  are  the  vertical  calculation  increment 
and  the  maximum  height  above  ground  level  (AGL)  to  which  the 
profile  should  be  calculated.  The  default  values  for  \'.hese  three 
parameters  are  50  feet,  1  foot,  and  10  feet,  respectively. 
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FIGURE  D- 


8  NON -DIHEU SIGNAL  NALL  JET  VERTICAL  VELOCITY  PROFILE 
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VELOCITY  PROFILE  STATUS  MENU 


CODE 

PARAMETER 

VALUE 

UNITS 

A 

PROFILE  STATION  POSITION 

50.00 

FT 

B 

VERTICAL  INCREMENT 

1.00 

FT 

C 

MAXIMUM  PROFILE  HEIGHT 

10.00 

FT 

D 

MINIMUM  BOUNDARY  LAYER  HEIGHT 

1.50 

FT 

E 

DATA  OUTPUT  FILENAME 

DFRC . PTS 

ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  ==> 
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An  option  is  also  provided  to  allow  the  user  to  specify  a  minimum 
boundary  layer  height  based  on  flight  test  data.  The  default 
value  of  1.5  feet  for  this  value  will  insure  that  the  peak 
velocity  on  the  profile  is  always  equal  to  or  greater  than  1.5 
feet  (at  very  close  positions  to  the  rotor,  the  ROTWASH 
methodology  often  calculates  an  unrealistically  thin  boundary 
layer  if  a  minimum  limit  is  not  specified) .  Each  of  the  four 
parameters  is  input  by  typing  the  appropriate  code,  <RETDRjr>,  the 
new  input  data  value,  and  <RETnRM>  to  end  the  sequence.  The  last 
parameter  on  the  menu  is  the  filename  for  data  that  is  written  to 
disk  if  the  graphics  file  toggle  switch  is  set  to  <Y>. 

The  wall  jet  velocity  profile,  calculated  by  the  program  using 
default  inputs,  is  presented  in  figure  D-11  when  the  <RBTURM9  key 
is  typed  by  itself.  Output  from  the  analysis  describes  the  shape 
of  both  the  mean  and  peak  velocity  profiles;  an  example  is 
presented  in  figure  D-12  correlated  with  Bell  XV-15  flight  test 
data  (this  particular  example  correlates  to  a  rotor  height  of 
37.5  feet  and  a  DFAC  value  of  66.1  feet  on  figure  D-10) . 

The  output  format  provides  velocity  profile  data  in  units  of 
either  feet  per  second,  knots,  or  pounds  per  square  foot  (also 
referred  to  as  dynamic  pressure). 

If  the  specified  increment  in  vertical  height  is  small  or  the 
maximum  calculated  height  is  large,  the  quantity  of  data  to  be 
output  to  the  screen  may  exceed  the  10-line  limit  for  1  screen 
frame.  When  this  situation  occurs,  the  typing  of  <C>  at  the 
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SINGLE  ROTOR  VELOCITY  PROFILE  AT  RADIUS  -  50.0  FT 


PROFILE 

BOUNDARY 

HEIGHT  - 

11.49 

FT 

KALF-VEL. 

HEIGHT  - 

4.10 

FT 

MAX-VEL 

HEIGHT  - 

1.15 

FT 

EIGHT 

MEAN  VELOCITY 

PEAK  VELOCITY 

MEAN  Q 

PEAK  Q 

(FT) 

(FPS) 

(KN) 

(FPS) 

(KN) 

(PSF) 

(PSF) 

.00 

.000 

.000 

.000 

.000 

.000 

.000 

1.00 

42.306 

25.077 

78.945 

46.796 

2.127 

7.407 

2.00 

40.282 

23.878 

81.231 

48.151 

1.928 

7.842 

3.00 

32.043 

18.994 

74.263 

44.021 

1.220 

6.554 

4.00 

24.823 

14.714 

65. 002 

38.531 

.732 

5.021 

5.00 

18.565 

11.005 

59.944 

35.533 

.410 

4.270 

6.00 

13.240 

7.848 

51.299 

30.408 

.208 

3.128 

7.00 

8.829 

5.233 

39.909 

23.657 

.093 

1.893 

8.00 

5.319 

3.153 

27.477 

16.288 

.034 

.897 

9.00 

2.700 

1.600 

15.691 

9.301 

.009 

.293 

10.00 

.964 

.571 

6.222 

3.688 

.001 

.04C 

TYPE  <C>ONTINUE,  NEXT  <P>OINT,  <N>EW  CASE,  E<X>IT  »=•> 
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prompt  (at  the  bottom  of  the  screen)  results  in  the  next  10  lines 
of  data  being  written  to  the  screen.  If  the  code  value  <P>  is 
typed,  the  program  returns  to  the  wall  jet  analysis  menu.  The 
typing  of  code  <W>  results  in  the  program  returning  to  the  master 
input  data  menu.  If  the  code  <X>  is  typed,  the  program  returns 
to  the  DOS  system  prompt. 

Flight  test  data,  correlated  with  output  from  the  wall  jet 
option,  are  presented  in  section  3  of  Volume  I  of  this  report  for 
the  Bell  XV-15,  Bell-Boeing  MV-22,  Sikorsky  CH-53E,  Sikorsky  SH- 
60B,  and  Canadair  CL-84.  Based  on  the  correlation  conducted  for 
references  D-1  and  D-2,  it  is  generally  recommended  that  the  wall 
jet  option  be  used  for  calculation  of  velocity  profiles  at 
distances  greater  than  1.5  times  the  rotor  radius  from  the  center 
of  the  rotor.  At  distances  less  than  this  value,  the 
mathematical  model  is  not  detailed  enough  to  predict  rotorwash 
flowfield  characteristics  accurately.  This  limitation  is  not 
serious  because  distances  closer  to  the  rotor  tip  than  1.5  times 
the  rotor  radius  have  little  practical  reason  for  being  analyzed 
for  rotorwash  effects  on  the  environment.  Collision  avoidance 
with  respect  to  objects  in  close  proximity  to  the  rotorcraft  is 
the  critical  issue  at  this  close  a  distance. 
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DL  =  12.71  PSF 
SIGP(o*)  s  1.01 
AZIMUTH  =  270  DEG 
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RIDL  =  13.0  PCT 
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PEAK  FUQKT  TEST  DATA 
MEAN  CALCUIATED  DATA 
PEAK  CALCULATED  DATA 
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FIGURE  D-12  XV-15  VELOCITY  PROFILE  CORRELATION 


D.1.5  THE  TWIN-ROTOR  INTERACTION  PLANE  OPTION 

The  twin-rotor  interaction  plene  option  calculates  the  velocity 
profile  contained  in  the  plane  which  is  oriented  perpendicular  to 
the  ground  and  to  the  line  segment  which  connects  the  center  of 
both  rotor  hubs  (refer  to  figure  D-1) .  After  choosing  the 
interaction  plane  option  on  the  velocity  analysis  option  menu, 
the  user  must  specify  the  same  parameters  on  the  velocity  profile 
status  menu  that  are  required  with  the  wall  jet  option.  The  only 
difference  between  the  parameters  is  the  reference  position  for 
specification  of  the  horizontal  location  of  the  velocity  profile 
with  respect  to  the  rotor  (the  first  menu  option) .  For  the 
interaction  plane  analysis,  this  distance  is  referenced  as  0  at 
the  intersection  of  the  interaction  plane  and  the  line  connecting 
the  rotors  and  not  directly  to  the  center  of  one  of  the  two 
rotors  (i.e.,  the  input  value  for  distance  along  the  interaction 
plane  (DAIP)  for  a  tiltrotor  is  along  a  line  that  is  an  extension 
of  the  fuselage  centerline  as  seen  in  figure  D-10) .  No  velocity 
profile  differences  are  assumed  to  exist  by  the  mathematical 
model  for  points  equidistant  along  the  interaction  plane  but  on 
opposite  sides  of  the  line  connecting  the  rotors.  For  a 
tiltrotor,  this  means  that  the  calculated  velocity  profiles  both 
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directly  in  front  of  and  directly  aft  of  the  aircraft  are  the 
sane  when  equidistant  along  the  interaction  plane. 


Example  output  from  the  interaction  plane  option  is  presented  in 
figure  D-13.  This  data  format  closely  resembles  that  of  the  vail 
iet  option  except  that  both  horizontal  and  vertical  velocity 
profile  components  arc  calculated.  The  horizontal  velocity 
component  is  identified  by  the  "H**  in  the  column  following  the 
height  column  and  the  vertical  component  by  the  **7" .  The 
mechanics  for  viewing  data  un  the  screen  (if  more  data  exist  than 
will  fit  on  one  screen  frame)  and  for  transferring  to  other  menus 
are  exactly  as  are  described  for  tlie  wall  jet  option. 


TWIN  ROTOR 

INTERACTION 

PLANE  VELOCITY  PROFl 

LE  AT 

STATION  = 

50.0  FT 

HEIGHT 

MEAN  VELOCITY 

PEAK  VELOCITY 

MEAN  Q 

PEAK  Q 

(FT) 

(FPS) 

(KN) 

(FPS) 

(KN) 

(PSF) 

(PSF) 

.  00 

H 

.  000 

.000 

.000 

.000 

.000 

.000 

V 

.  000 

.000 

.000 

.000 

.000 

.000 

1 . 00 

H 

63.626 

37.715 

103.992 

61.643 

4.811 

12.852 

V 

21.760 

12.8^9 

35.555 

21.082 

.  563 

1.503 

'  0 

H 

66.665 

39.517 

108.960 

64.588 

5.282 

14.110 

V 

24.133 

14  .  305 

39.443 

23.381 

.f92 

1.849 

3 . 00 

H 

65.800 

39.004 

107.545 

63.749 

5.146 

13.746 

V 

25.  13t 

14  .  900 

41.082 

24 .352 

.751 

?.006 

O 

o 

H 

64 . 912 

38.478 

106.093 

62 . 889 

5.008 

13.377 

V 

26 . 094 

15.468 

42.649 

25.281 

.809 

2.162 

TYPE 

<C>0NTINb’E, 

NEXT  <P>OINT,  <N>EW 

CASE, 

E<X>IT  -•»> 
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Output  from  the  interaction  plane  option,  correlated  with  XV-15 
flight  trst  data,  is  presented  in  figure  D-14  as  an  example. 

These  data  are  excerpted  from  reference  D-2.  Other  original 
flight  test  data  for  both  the  Bell  XV-15,  Bell-Boeing  MV-22,  and 
Canadair  CL-84  (tiltwing)  are  dccui4onted  in  section  3  of  Volume  I 
of  this  report. 

Several  practical  recommendations  need  to  be  noted  for  users  of 
the  interaction  plane  option.  When  analyzing  tiltrotor  and 
tilewing  configurations,  the  user  must  be  careful  to  avoid 
choosing  analysis  locations  that  are  coincident  with  components 
of  the  nose  or  tail  structure  on  the  aircraft.  Also,  locations 
in  close  proximity  to  the  nose  or  tail  of  a  real  aircraft  should 
not  be  «.xpected  to  have  rotorwash  flowfield  characteristics 
identical  to  calculated  velocity  profiles.  .\t  these  locations 
unmodeled  airframe  aerodynamic  interferences  significantly 
influence  the  rotorwash  flowfield  structure.  Flight  tost  data 
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obtdined  from  the  XV-15,  MV-22,  and  CL-84  indicate  that  measured 
mean  and  peak  velocity  profiles  at  points  equidistant  along  the 
interaction  plane  (both  in  front  of  and  aft  of  an  actual 
aircraft)  often  do  not  yield  identical  results  as  might  be 
expected.  TherL',ore,  output  from  this  analysis  option  should  be 
calibrated  with  flight  test  data  whenever  possible  in  an  attempt 
to  determine  whether  positions  forward  or  aft  of  the  aircraft  may 
be  move  critical  for  analysis. 

D.l.e  GROUND  VORTEX  ANALYSIS.  OPTION 

A  ground  vortex  structure  is  formed  when  ambient  wind  and/or 
rotorcraft  translational  velocity  overcome  the  rotor-induced  wall 
jet  flowfieid.  A  diagram  of  the  ground  vortex  is  presented  in 
figure  D-15.  Due  to  the  elementary  nature  of  the  mathematical 
model  formulation  used  in  the  ROTWASH  program,  the  ground  vortex 
option  should  be  used  with  caution.  As  discussed  in  section  2  of 
Volume  I  of  this  report,  almost  no  test  data  exist  to  validate 
the  mathematical  model.  Also,  the  single  main  rotor  helicopter 
is  the  only  configuration  that  can  be  analyzed  with  the  model  as 
presently  formulated.  Since  all  examples  presented  up  to  this 
point  in  the  user's  guide  have  been  for  the  XV-15,  it  is 
necessary  to  d.2f ine  input  data  for  a  single  main  rotor  helicopter 
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GROUND  VORTEX  STRUCTURAL  CHARACTERISTICS 


ROTWASH  USER  INPUT  DATA  MENU 

CODE 

PARAMETER 

VALUE 

UNITS 

A 

NUMBER  OF  ROTORS  (1  OR  2) 

1 

-ND- 

B 

HUB  TO  HUB  ROTOR  SEPARATION 

.0 

FT 

C 

ROTOR  RADIUS 

39.5 

FT 

D 

GROSS  WEIGHT 

56000.0 

LB 

E 

FUSELAGE  DOWNLOAD  FACTOR 

5.0 

PCT 

F 

ROTOR  HEIGHT  ABOVE  GROUND 

30.0 

FT 

G 

SHAFT  TILT  ANGLE  (<20  DEG) 

.0 

DEG 

H 

AIR  DENSITY  RATIO 

1.0000 

ND 

I 

AMBIENT  WIND  (-10  TO  10  KT) 

.0 

KT 

ENTER 

DATA  ENTRY  CODE  OR  <RETURN>  TO 

CONTINUE  =^> 

riQURB  D-l«  CH-S3E  INPUT  DATA  FOR  THE  GROUND  VORTEX  EXitMPLE 


before  the  ground  vortex  option  is  explained  in  detail.  The 
Sikorsky  CH-53E  configuration  has  been  chosen  for  this  task.  The 
main  input  data  menu  for  this  configuration,  as  typed  into  the 
program,  is  presented  in  figure  D-16. 

The  ground  vortex  analysis  option  is  specified  by  selection  of 
the  character  -<0>  on  the  velocity  analysis  option  menu  as  shown 
in  figure  D-5.  The  user  is  then  required  to  complete  the 
ground/disk  vortex  input  data  menu  which  is  presented  in  figure 
D-17.  Two  of  the  parameters  specified  on  this  menu  are 
rotorcraft  configuration  parameters.  These  parameters  are  the 
rotor  tip  speed  ( feet/second)  and  the  number  of  rotor  blades. 

The  next  parameter  to  be  specified  is  the  rotorcraft 
translational  velocity  with  respect  to  the  surrounding  air  mass 
(i.e.,  an  input  of  15  knots  can  be  either  15  knots  ground  speed 
on  a  no-wind  day  or  0  knots  ground  speed  on  a  day  with  a  15-knot 
headwind) .  Each  of  these  values  is  input  with  the  keyboard  using 
the  same  techniques  that  have  been  previously  discussed. 

The  next  four  menu  parameters  define  the  position  in  three- 
dimensional  space  (feet)  where  the  velocity  profile  will  be 
calculated  (see  figure  D-18) .  The  positive  directions  for  the 
coordinate  system  are  aft  and  right  from  the  center  of  the  rotor. 
Therefore,  in  order  to  calculate  a  slice  of  the  ground  vortex 
directly  in  front  of  the  rotor,  a  negative  X-value  (longitudinal 
position)  is  input  along  with  a  aero  Y-value  (lateral  position) . 
The  Z-axis  increment  (feet)  and  maximum  calculation  height  (feet) 
parameters  define  the  number  of  points  that  are  calculated 
between  ground  level  and  the  maximum  height  of  interest. 
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GROUND/DISK  VORTEX  INPUT 

DATA  MENU 

(FOR  SINGLE  MAIN  ROTOR  HELICOPTERS  ONLY) 

CODE 

PARAMETER 

VALUE 

UNITS 

A 

ROTOR  TIP  SPEED 

733.00 

FPS 

B 

NUMBER  OF  ROTOR  BLADES 

7.00 

-ND- 

C 

TRANSLATIONAL  SPEED 

16.00 

KTS 

D 

XT  POSITION 

-60.00 

FT 

E 

YT  POSITION 

.00 

FT 

F 

ZT  CALCULATION  INCREMENT 

2 . 00 

FT 

G 

MAXIMUM  CALCULATION  HEIGHT 

20.00 

FT 

ENTER 

DATA  ENTRY  CODE  OR  <RETURN> 

TO  CONTINUE  = 

S> 

rianRB  d-i7  onoairo/DxsK  vortbz  input  data  mbmd 


After  each  input  data  parameter  has  been  specified,  the  user 
types  the  <RETURM>  key  to  initiate  the  analysis.  The  response  of 
the  program  is  to  list  three  calculated  parameters  which  are 
followed  by  a  prompt.  An  example  of  this  response  is  provided  in 
figure  D-19.  The  calculated  values  are  the  nondimensionalized 
rotor  height  above  ground  and  two  advance  ratio  parameters.  The 
program  prompt  following  the  screen  output  requires  the  user  to 
input  the  ground  vortex  strength  ratio  which  is  obtained  from  the 
graph  in  figure  D'‘20  using  the  three  calculated  parameters. 
Background  on  the  use  of  this  figure  is  discussed  in  section  2  of 
Volume  I  of  this  report.  The  limits  presented  on  the  graph  in 
figure  D-20  define  the  advance  ratio  range  wherein  the  ground 
vortex  would  be  expected  to  occur.  At  advance  ratios  much  less 
than  0.035,  the  ground  vortex  does  not  have  favorable  conditions 
for  formation.  At  advance  ratios  slightly  greater  than  0.055, 
the  ground  vortex  is  dispersed  by  the  rotor  because  the 
translational  velocity  relative  to  the  air  mass  is  too  high  for 
the  vortex  to  maintain  position.  After  the  user  has  entered  the 
ground  vortex  strength  ratio,  the  program  calculates  the  ground 
vortex  circulation  and  the  position  of  the  ground  vortex  core 
with  respect  to  the  axis  system  presented  in  figure  D-i8.  This 
output  is  also  presented  in  figure  D-19  below  the  prompt  for  the 
ground  vortex  strength  ratio. 
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ROTOR  HEIGHT  ABOVE  GROUND  H/D  .3797 
ADVANCE  RATIO  MU-STAR  .5373 
ADVANCE  RATIO  MU  .0368 


ENTER  GROUND  VORTEX  STRENGTH  RATIO 

(SKE  FIGURE  D-20)  3. 


GROUND  VORTEX  CORE  POSITION 

X-LOCATION  (XXGV)  =  -59.14  FT 

Y-LOCATION  (ZZGV)  =  9.15  FT 

GROUND  VORTEX  CIRCULATION  =  732.35  FT**2/SEC 


PRESS  <RETURN>  TO  CONTINUE 


FIGURE  D-19  GROUND  VORTEX  ANALYSIS  INTERMEDIATE  OUTPUT 


Engineering  data  from  the  ground  vortex  analysis  option  is 
obtained  by  typing  the  <RETURN>  key  after  the  vortex  position 
parameters  are  displayed.  An  example  of  the  output  format  is 
presented  in  figure  0-21.  Calculated  field  velocities  at  the 
various  points  along  the  profile  Z-axis  are  presented  In  both  a 
vectorial  XYZ  component  format  and  as  a  total  resolved  magnitude 
in  units  of  feet  per  second  and  knots.  These  same  data  are  also 
provided  to  the  user  as  dynamic  pressures  in  units  of  pounds  per 
square  feet. 

D.1.7  DISK  EDGE  VORTEX  ANALYSIS  OPTION 

The  disk  edge  vortex  analysis  option  was  developed  to  provide  a 
capability  to  estimate  the  strength  of  trailing  vortices  behind 
helicopters  in  forward  flight  as  described  by  figure  D-22.  Like 
the  ground  vortex  option,  this  option  is  limited  to  use  with  the 
single  main  rotor  helicopter  configuration.  Approximations  of 
vortex  core  size  are  not  calculated  by  the  mathematical  model  and 
must  be  estimated  using  flight  test  data.  Available  flight  test 
data  are  presented  in  section  3.6  of  Volume  I  of  this  report. 

All  examples  presented  in  the  user's  guide  for  this  analysis 
option  utilize  the  same  CH-53E  input  data  array  that  was 
described  previously. 
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Traatiiion  from  recirculation 


-20  CALCUIJtXED  GROUND  VORTEX  CIRCULATION  STRENGTH 


HEIGHT 

MEAN 

VELOCITY 

MEA’,  Q 

(FT) 

(FPS) 

(KN) 

(PSF) 

o 

o 

X 

-24.582 

-14.572 

.718 

Y 

.  000 

.000 

.000 

Z 

.000 

.000 

.000 

T 

24.582 

14.572 

.718 

2.00 

X 

-25.790 

-15.288 

.790 

Y 

.  000 

.  000 

.000 

Z 

1.013 

.601 

.001 

T 

25.810 

15.299 

.792 

4.00 

X 

-30.173 

-17.886 

1.082 

Y 

.000 

.000 

.000 

Z 

2.668 

1.700 

.010 

T 

30. 30S 

17.966 

1.092 

TYPE  <C>ONTINU 

E,  NEXT 

<P>OINT,  <N>EW  CASE,  E<X>1T  ==> 

FiaURE  D-21  GROnSTD  VORTEX  VELOCITY  FIELD  OUTPUT  DATA 


Exerucion  of  the  uisk  edge  vortex  option  is  initiated  by  choosing 
<D>  the  velocity  analysis  option  menu  shown  in  figure  D-5. 

The  g.‘:‘ond/disK  vortex  input  data  menu,  shown  in  figure  0-17, 
must  bo  completea  as  described  in  the  previous  section.  All 

sign  ccr-ventions  used  in  the  specification  of  locations  for  the 
calculation  of  field  velocities  are  the  same  as  for  the  ground 
voitex  option.  After  the  input  data  menu  is  completed,  the 
analysis  option  is  executed  by  typing  the  <RETURN>  key. 

The  initial  program  response,  figure  D-23,  is  to  write  to  the 
screen  the  calculated  values  for  the  vortex  circulation  and  the 
settling  angle  of  the  trailing  vortex  components  as  defined  in 
figure  D-22  for  a  forward-flight  velocity  of  50  knots  and  a  rotor 
height  above  ground  of  100  feet.  When  executing  this  option,  it 
is  important  that  the  user  confirm  that  the  calculated  settling 
angle  is  less  than  approximately  20  degrees.  At  settling  angles 
larger  than  this  value,  the  airspeed  of  the  helicopter  is 
probably  too  slow  to  sustain  the  formation  of  the  trailing  edge 
vortex  system  which  is  predicted  using  this  mathematical  model, 
and  results  should  be  considered  suspect.  User  specified 
velocity  field  calculations  are  presented  after  the  •<RETURM>  key 
is  typed  in  the  same  format  as  was  discussed  with  the  ground 
vortex  option.  An  example  output  for  the  CH-53E  is  presented  in 
figure  D-24. 
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Bound  Vortex 


FIGURE  D-22  HORSESHOE  VORTEX  GEOMETRY  FOR  CALCULATION  OF  FORWARD  FLIGHT 
hazard  potential 


other  features  of  the  nathematical  model  which  should  be  noted 
are  described  below. 

1.  The  calculated  trailing  vortex  strength  does  not  decay  as 
a  function  of  increasing  distance  from  the  helicopter  or 
with  increasing  time.  This  weakness  in  the  model  will 
affect  correlation  with  flight  test  data  and  could  result 
in  predicted  values  that  are  greater  than  measured  values 
for  field  velocities. 

2.  The  decay  of  the  trailing  vortex  structure  after 
impingement  with  the  ground  is  not  modeled.  Any 
prediction  of  field  velocities  behind  the  initial 
impingement  point  should  not  be  considered  valid.  The 
location  of  the  impingement  point  with  respect  to  the 
ground  and  the  location  of  points  specified  for  velocity 
field  calculations  must  be  checked  by  hand  calculation  to 
ensure  that  geometry  constraints  are  not  violated.  This 
is  accomplished  by  using  the  rotor  height  and  the  settling 
angle  to  calculate  the  horizontal  distance  behind  the 
helicopter  where  the  impingement  occurs. 

D.1.8  PERSONNEL  OVERTURNING  FORCE  AND  MOMENT  ANALYSIS 

The  personnel  overturning  force  and  moment  analysis  model  is 
formulated  for  use  with  both  the  single  main  rotor  and  twin-rotor 
configurations.  The  initial  task  of  the  model  is  to  calculate 
the  velocity  profile  for  a  specified  location.  The  calculated 
velocity  profile  is  then  integrated  over  the  projected  area  of  a 
human  body  to  obtain  estimates  of  the  applied  aerodynamic  force 
and  moment.  This  analysis  technique  is  summarized  in 
figure  D-25. 

Use  of  the  overturning  force  and  moment  option  is  initiated  by 
choosing  <H>  (for  hazard)  on  the  program  logic/comment  menu. 

This  is  followed  by  the  choice  of  <M>  (for  overturning 
force/raoment)  on  the  hazard  analysis  option  menu  as  shown  in 
figure  D-6.  The  user  then  specifies  the  parameters  listed  on  the 
overturning  force  and  moment  data  menu  presented  in  figure  D-26. 
The  first  parameter  on  this  menu  specifies  the  use  of  either  the 
wall  jet  or  the  interaction  plane  analysis  for  creation  of 
velocity  profile  data.  The  second  option  specifies  the  use  of 
either  the  "large”  (6  feet  in  height)  or  "small"  (4  feet  in 
height)  human  body  mathematical  model.  The  third  parameter 
provides  the  user  the  capability  to  specify  a  graphics  output 
filename  (assuming  this  option  has  been  toggled  ON  using  the 
program  logic/comment  menu) .  If  the  user  executes  the  option 
without  changing  the  filename  and  a  file  already  exists  with  the 
same  filename,  the  user  is  notified  and  required  to  change  the 
filename.  Three  of  the  last  four  menu  tables  define  the 
locations  that  are  to  be  analyzed  using  Lite  option.  These 
variables,  all  in  units  of  feet,  are  the  initial  station  position 
for  analysis,  the  increment  in  station  position,  and  the  final 
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DISK  VORTEX  VELOCITY  PROFILE  DATA 


X-LOCATION  (XT) 
Y -LOCATION  (YT) 


150.00  FT 
.00  FT 


VORTEX  CIRCULATION 
VORTEX  CIRCULATION 
5-M  INITIAL  CIRCULATION  = 
SETTLING  ANGLE 


5441.24  FT**2/SEC 
505.51  M**2/SEC 
232.37  M**2/SEC 
9.29  DEG 


PRESS  <RETURN>  TO  CONTINUE 


FIOURB  D-23  DISK  EDGE  VORTEX  OPTION  ZNTERNBDIATB  OUTPUT 


HEIGHT 

MEAN 

VELOCITY 

MEAN  Q 

(FT) 

(FPS) 

(KN) 

(PSF) 

.  00 

X 

-4.052 

-2.402 

.020 

Y 

.  000 

.000 

.000 

Z 

.  000 

,000 

.  000 

T 

4 . 052 

2.402 

.020 

50.00 

X 

-6,444 

-3.820 

.049 

Y 

.  000 

.000 

.000 

Z 

-27 . 295 

-16.180 

.885 

T 

28.046 

16.625 

.935 

00.00 

X 

-5.735 

-3.399 

.039 

Y 

.  000 

.  000 

.000 

Z 

-30,049 

-17.812 

1.073 

T 

30.592 

18.134 

1.112 

TYPE  <C>ONTINUE,  NEXT  <P>OINT,  <N>EW  CASE,  E<X>IT  ==> 


PIOURE  D-24  VELOCITY  FIELD  OUTPUT  DATA  FROM  THE  DISK  EDGE 

VORTEX  OPTION 
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FIOURB  D>2S  OVERTURNINQ  FORCE  AND  MOMENT  CALCDLATZON  FROCBDDRB8 


station  position,  respectively.  The  fourth  parameter  is  the  user 
option  to  specify  a  minimum  boundary  layer  height  (discussed  in 
section  D.1.4  of  the  user's  guide).  The  mechanics  for  input  of 
the  desired  values  using  this  menu  are  as  described  for  previous 
menus. 

After  the  analysis  is  executed  by  typing  the  <RBTURN>  key,  the 
calculated  results  are  written  out  in  the  format  presented  in 
figure  D-27.  Three  columns  of  data  are  written  using  this 
format.  The  first  column  identifies  either  the  distance  from 
rotor  center  (DFRC).  for  the  wall  jet  option  or  the  distance  along 
the  interaction  plane  (DAIP)  for  the  interaction  plane  option. 

The  second  and  third  columns  are  the  associated  total  force  and 
total  moment  values  calculated  for  the  projected  area  of  a  human 
body.  Example  data  for  the  XV-15  using  this  option  are  presented 
in  figure  D-28  for  reference.  At  the  bottom  of  the  screen,  the 
user  is  required  to  return  to  the  previous  menu  by  typing  <P> , 
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OVERTURNING  FORCE/MOMENT 

DATA 

MENU 

CODE 

PARAMETER 

VALUE 

UNITS 

A 

<W>ALL  JET  CR  <I>NTERACTION  PLANE 

W 

B 

<L>ARGE  OR  <S>MALL  PERSON 

L 

C 

DATA  OUTPUT  FILENAME 

OTDFRC.PTS 

D 

INITIAL  STATION  POSITION 

50 .  DO 

FT 

E 

HORIZONTAL  INCREMENT 

10.00 

FT 

F 

MAXIMUM  STATION  POSITION 

100.00 

FT 

G 

MINIMXm  BOUNDARY  LAYER  HEIGHT 

1.50 

FT 

ENTER 

DATA  ENTRY  CODE  OR  <RETURN>  TO 

CONTINUE  ==> 

riaURB  D-26  CVERTORNINQ  FORCF  lOfD  MOMEKS  DATA  MENU 


SUMMARY  OF  OVERTURNING  FORCES  AND  MOMENTS 


RADIUS 

TOTF 

TOTM 

(FT) 

(LB) 

(FT-LB) 

50.00 

42.721 

114.596 

60.00 

33 . 708 

91.718 

70.00 

25.046 

69.346 

80.00 

17.520 

49.892 

90.00 

11.511 

33.382 

100.00 

6.926 

20.269 

TYPE  <C>ONTINUE,  NEXT  <P>OINT,  <N>EW  CASE,  E<X>IT  == 


FIGURE  D-27  SIMPLIFIED  OUTPUT  FORMAT  FOR  THE  OVERTURMIMQ 

FORCE  AND  MOMENT  ANALYSIS 
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the  master  input  data  menu  by  typing  <H>,  or  to  exit  the  program 
by  typing  <X>.  By  typing  <C>,  the  user  can  view  the  next  video 
screen  of  data  if  more  than  one  screen  of  data  is  generated. 
Otherwise,  the  <C>  option  works  identically  to  the  <P>  option. 

The  user  can  examine  the  detailed  calculations  used  to  create  the 
summary  output  by  specifying  the  desired  analysis  location  as  the 
initial  station  position.  This  input  must  then  be  followed  by 
specification  of  the  increment  value  as  0.0  or  only  the  summary 
output  will  appear  (as  shown  in  figure  D-29).  The  resulting 
"large"  person  output  using  this  option  is  presented  on  two 
screen  frames  as  shown  in  figure  D~30.  The  first  output  frame 
presents  a  summary  of  the  velocity  profile  calculations  as  a 
function  of  height  above  ground  at  the  specified  station 
position.  The  second  output  frame  presents  the  associated 
calculations  for  dynamic  pressure,  overturning  force,  and 
overturning  moment.  The  last  two  columns  in  the  second  table  are 
values  of  total  force  and  moment  summed  for  the  incremental 
increase  in  height. 

In  this  example,  the  force  and  moment  values  of  42.7  pounds  and 
114.6  foot-pounds,  respectively,  at  5.75  feet  are  the  total  force 
and  moment  values  that  would  normally  be  printed  out  in  the 
summary  output  (figure  D-27).  These  values  are  checked  by 
totaling  the  individual  height-related  values  in  the  overturning 
force  and  moment  columns  (second  and  third  columns) .  An  example 
output  for  the  second  screen  of  the  "small"  person  option  is 
presented  in  figure  D-31  for  reference. 

Both  qualitative  and  quantitative  overturning  force  and  moment 
data  are  presented  in  section  5  of  volume  I  of  this  report 
correlated  with  ROTWASH  program  output  for  the  Bell  XV-15, 
Sikorsky  CH-53E,  and  the  Sikorsky  S-61.  These  calculated  data 
all  assume  a  coefficient  of  drag  for  a  human  body  of  1.1  (which 
according  to  Hoerner,  reference  D-6,  can  vary  from  1.0  to 
approximately  1.3). 

D.1.9  PARTICULATE  CLOUD  ANALYSIS  OPTION 

The  methodology  used  in  the  calculation  of  particulate  cloud  size 
is  presented  in  section  5.8  of  volume  I  of  this  report  along  with 
a  very  limited  amount  of  flight  test  data.  This  option  is 
applicable  to  both  single  main  rotor  and  twin-rotor 
configurations.  The  particulate  cloud  geometry  utilized  in  the 
analysis  option  is  presented  in  figure  D-32. 

The  particulate  clcud  option  is  initiated  with  the  typing  of  <C> 
on  the  hazard  analysis  option  menu  as  shown  in  figure  D-6.  The 
screen  that  i  .  •  itten  subsequently  presents  the  user  with  a 

prompt  for  inp  of  the  terrain  erosion  factor.  The  value  for 
this  factor  is  chosen  from  the  graph  in  figure  D-33. 
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OVERTURNING  FORCE,  POUNDS 


ROTOR  HEIGHT,  XV-15  XV-15 

FEET  FUGHT  TEST  DATA  CALCUU^TED  DATA 


0  20  40  60  80  100 


IPLANE  DISTANCE  (DAIP),  FEET 

FiaURB  D-28  BELL  X7-15  OVERTTOMINO  FORCE  AS  A  FUNCTION  OF 
DISTANCE  ALONG  THE  INTERACTION  FLANS  (0  DEGREE  AZIMUTH) 


OVERTURNING  FORCE/MOMENT 

DATA 

MENU 

CODE 

PARAMETER 

VALUE 

UNITS 

A 

<W>ALL  JET  OR  <I>NTERACTION  PLANE 

W 

B 

<L>ARGE  OR  <S>MALL  PERSON 

L 

C 

DATA  OUTPUT  FILENAME 

OTDFRC. PTS 

D 

INITIAL  STATION  POSITION 

50 . 00 

FT 

E 

HORIZONTAL  INCREMENT 

.  00 

FT 

F 

MAXIMUM  STATION  POSITION 

100 . 00 

FT 

G 

MINIMUM  BOUNDARY  LAYER  HEIGHT 

1 . 50 

FT 

ENTER 

DATA  ENTRY  CODE  OR  <RETURN>  TO 

CONTINUE  ==> 

FIGURE  D-29  MENU  SPECIFICATION  OF  DETAILED  FORCE /MOMENT  OUTPUT 
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SINGLE  ROTOR  VELOCITY 

PROFILE 

AT  RADIUS 

-  50.0 

FT 

:iGHT 

MEAN  VELOCITY 

PEAK  VELOCITY 

MEAN  Q 

PEAK  Q 

FT) 

(FPS) 

(KN) 

(FPS) 

(KN) 

(PSF) 

(PSF) 

.25 

34.705 

20.572 

64.762 

38.389 

1.431 

4.984 

.75 

40.602 

24.068 

75.767 

44.912 

1.959 

6.822 

1.25 

43.676 

25.890 

81.502 

48.312 

2.267 

7.894 

1.75 

42.515 

25.202 

82.535 

48.924 

2.148 

8.096 

2.25 

38.121 

22.597 

79.742 

47.269 

1.727 

7.557 

2.75 

34.004 

20.157 

7S.248 

45.197 

1.374 

6.909 

3.25 

30.146 

17.870 

72.134 

42.759 

1.080 

6.184 

3.75 

26.537 

15.730 

67.491 

40.007 

.837 

5.413 

4.25 

23.170 

13.734 

63.591 

37.695 

.638 

4.806 

4.75 

20.041 

11.880 

61.475 

36.441 

.Ml 

4.491 

5.25 

17.147 

10.164 

58.134 

34.460 

.349 

4.016 

5.75 

14.485 

8.586 

53.785 

31.882 

.249 

3.438 

TYPE  <RETURN>  TO  CONTINUE 


SINGLE  ROTOR 

FORCE  PROFILE  AT  RADIUS  =« 

50.0 

FT 

HEIGHT 

PEAK  Q 

FOVER 

OVERM 

TOT  F 

TOT  M 

(FT) 

(PSF) 

(LB) 

(FT-LB) 

(LB) 

(FT-LB) 

.25 

4.984 

3.016 

.754 

3.016 

.754 

.75 

6.822 

4.128 

3.096 

7 . 143 

3.850 

1.25 

7.894 

4.776 

5.970 

11.919 

9.820 

1.75 

8.096 

4.898 

8.571 

16.817 

18.391 

2.25 

7.557 

4.572 

10.287 

21.389 

28.678 

2.75 

6.909 

4 . 180 

11.495 

25.569 

40.174 

3.25 

6.184 

3.741 

12.159 

29.311 

52.333 

3.75 

5.413 

3.275 

12.282 

32.586 

64 . 615 

4.25 

4.806 

2.908 

12.357 

35.493 

76.972 

4.75 

4 . 491 

2.717 

12.907 

38.211 

89.879 

5.25 

4.016 

2.430 

12.757 

40.641 

102 . 636 

5.75 

3.438 

2.080 

11.960 

42.721 

114 . 596 

TYPE  <C>ONTINUE, 

NEXT  <P>OINT 

,  <N>EW  CASE, 

E<X>IT 

==> 

riOURE  D-30  DETAILED  OUTPUT  FORMAT  FOR  THE  "LAROE"  PERSON 

ANALYSIS  OPTION 


D-30 


SINGLE  ROTOR  FORCE  PROFILE  AT  RADIUS  =  50.0  FT 


HEIGHT 

PEAK  Q 

FOVER 

OVERM 

TOT  F 

TOT  M 

(FT) 

(PSF) 

(LB) 

(FT- LB) 

(LB) 

(FT-LB) 

.25 

4.984 

1.919 

.480 

1.919 

.480 

.75 

6.822 

2.627 

1.970 

4.546 

2.450 

1.25 

7.894 

3.039 

3.799 

7.585 

6.249 

1.75 

8.096 

3.117 

5.455 

10.702 

11.703 

2.25 

7.557 

2.909 

6.546 

13. 611 

18.250 

2.75 

6.909 

2.660 

7.315 

16.271 

25.565 

3.25 

6.184 

2.381 

7.738 

18.652 

33 .303 

3.75 

5.413 

2.084 

7.816 

20.736 

41.118 

TYPE 

<C>ONTINUE, 

NEXT  <P>0INT 

,  <N>EW 

CASE,  E<X>IT 

==> 

FIGURE  D-31  DETAILED  OUTPUT  FORMAT  FOR  THE  SECOND  SCREEN  FRAME 
OF  THE  "SMALL"  PERSON  ANALYSIS  OPTION 


D-31 


(^>l)  jopej  uo(SOJ3  uiejjai 


General  Terrain  Type 

FIGURE  D-33  APPROXIMATE  VALUES  FOR  THE  TERRAIN  EROSION 
IDENTIFIED  IN  TBE  LITERATURE 


The  output  format  for  the  particulate  cloud  option  is  designed 
for  both  single  main  rotor  and  twin-rotor  aircraft  as  shown  in 
figure  D-34.  In  this  example,  using  the  XV-15,  the  particulate 
cloud  boundaries  located  at  the  90-  and  270-degree  assimuths  (out 
the  span  of  the  wing)  are  specified  by  the  single  rotor  or  "SR" 
row  of  output  values.  The  interaction  plane  or  "IP"  row  defines 
the  cloud  boundaries  which  exist  straight  out  in  front  of  and 
directly  aft  of  the  aircraft  along  the  aircraft  centerline.  If  a 
single  main  rotor  configuration  is  being  evaluated,  the  "IP"  row 
in  the  printout  will  contain  all  zeros. 


ENTER  TERRAIN  EROSION  FACTOR  (-ND-)  ==>  0.4 


SUMMARY  OF  CLOUD  BOUNDARIES 


RC  AND 

RV  ARE  FROM 

ROTOR  CENTER 

(FT) 

RC 

RV 

ZV 

HC 

SR 

80.1 

62 . 8 

26.3 

37 .6 

IP 

113.0 

88.7 

37 . 2 

53 . 0 

QSMAX  =  13.5  PSF 


FIGURE  0-34  PARTICULATE  CLOUD  ANALYSIS  OPTION  OUTPUT 


D.2  ROTWASH  PROGRAM  DATA  OUTPUT  FILE  FORMATS 

Four  data  output  file  formats  can  be  specified  from  the  ROTWASH 
program  for  use  with  computer  graphics  programs.  Two  of  these 
output  file  formats  are  generated  by  the  wall  jet  and  interaction 
plane  velocity  profile  analysis  options.  The  other  two  formats 
are  generated  by  the  personnel  overturning  force  and  moment 
option.  These  file  formats  save  the  summary  force  and  moment 
data  for  both  the  wall  jet  and  interaction  plane  cases. 

The  first  two  lines  in  each  of  the  four  file  formats  are  user- 
specified  comments.  These  two  comment  lines  are  typed  in  through 
use  of  the  program  logic/comment  menu.  The  rest  of  the  data  in 
each  of  the  file  formats  is  either  header  information  or 
engineering  data.  The  example  files  presented  in  this  section 
are  written  for  direct  input  to  the  TECPLOT  Graphics  Program 
which  is  written  by  AMTEc  Engineering.  This  graphics  program  is 
one  of  several  IBM  PC/PC-compatible  engineering  graphics  programs 


presently  on  the  market.  Uaers  that  do  not  have  access  to  this 
program  can  easily  modify  the  ROTWASH  FORTRAN  code  for  other 
types  of  graphics  programs  by  modifying  the  appropriate  write 
statements.  Table  1  provides  the  user  with  a  cross  reference  of 
the  figure  number  for  each  of  the  four  output  types  and  the 
source  location  for  the  associated  FORTRAN  code  that  can  be 
modified  (see  program  listings  in  appendix  E) . 


TABLE  D-1  GRAPHICS  PXLB/SOtTRCB  CODE  REFERENCE  MATRIX 


FILE  OUTPUT  TYPE 

FIGURE 

SOURCE  CODE 
LOCATION 

Wall  Jet  Velocity  Profile  Output 

D-35 

Subroutine  WJVEL 

Interaction  Plan  Velocity  Profile 
Output 

D-36 

Subroutine  IPVEL 

Overturning  Force/Moment  Summary 
(Wall  Jet) 

D-37 

Subroutine 

HWJVEL 

Overturning  Force/Moment  Summary 
(Interaction  Plan) 

D-38 

Subroutine 

HIPVEL 

D-34 


XV- 15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA 
GROSS  weight  might  BE  ONE  OP  THE  COMMDIT  STRINGS 


TITLE-" VELOCITY  PROFILE, 
VARIABLES  -  X,HT 
ZONE  T  -  "MEAN  PROFILE, 
.0  .00 

25.1  1.00 

23.9  2.00 

19.0  3.00 

14.7  4.00 

11.0  5.00 

7.8  6.00 

5.2  7.00 

3.2  8.00 

1.6  9.00 

.6  10.00 

ZONE  T  -  "PEAK  PROFILE, 

.0  .00 

46.8  1.00 

48.2  2.00 

44.0  3.00 

38.5  4.00 

35.5  5.00 

30.4  6.00 

23.7  7.00 

16.3  8.00 

9.3  9.00 

3.7  10.00 

ZONE  T  -  “PEAK  Q,  PSF" , 

.0  .00 

7.4  1.00 

7.8  2.00 

6.6  3.00 

5.0  4.00 

4.3  5.00 

3.1  6.00 

1.9  7.00 

.9  8.00 

.3  9.00 

.0  10.00 


DFRC  -  SO.O  PT,  GW  -  13000  LB, 
KTS",  I-ll,  F-POINT 


KTS",  I-ll,  F-POINT 


I-ll,  F-POlNT 


WAGL  -  25.0  FT" 


riQURB  D-35  EXAMPLE  WALL  JET  OPTIOM  0RAPHIC8  FILE  FORMAT 


D-35 


)CV-15  CHAHACTERISTICS  ARE  USED  AS  IMDin« 

CROSS  HEIGHT  HIGHT  BE  OKE  OF  THE  COkSSt  ^INCS 


ZONE  T 

-  "MEAN 

.0 

.00 

37.7 

1.00 

39.5 

2 . 00 

39.0 

3 . 00 

38.5 

4.00 

37.9 

5.00 

37. < 

6.00 

36.8 

7.00 

36.3 

8.00 

35.7 

9.00 

35.1 

10.00 

ZONE  T 

-  "PEAK 

.0 

.  00 

61.6 

1. 00 

64.6 

2.00 

63.7 

3 . 00 

62.9 

4. 00 

62.0 

5.00 

61.1 

6.00 

60.2 

7.00 

59.3 

8.00 

58.3 

9.00 

57.4 

10.00 

ZONE  T  -  "PEAK  ; 

.  0 

.00 

12 . 9 

1.00 

14.1 

2.00 

13  .  7 

3.00 

13.4 

4 . 00 

13 . 0 

5.00 

12 . 6 

6.00 

12 . 3 

7.00 

11.9 

8.00 

11.5 

9.00 

11.1  10,00 


D-3.  IKTBR*OTIO«  OPTIOH  OPXPHICS  PI^p  POK^X 


C^oll  Sfx^nrcSx'’^- ..3 

=  -•«  -•  — T 

91.72 
69.35 
49 . 89 
33.38 
20.27 


60.00 
70.00 
80.00 
90.00 
100. 00 


33 . 71 
25.05 
17.52 
11.51 
6.93 


D-36 


XV- 15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA 
GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  COMMENT  STRINGS 


TITLE="TWIN  ROTOR  DAIP  DATA** 


VARIABLES 

=  DAIP, 

TOTF,TOTM 

ZONE  T  = 

"GW  =  13  000  LB,  WAGL  -  25.0  FT",  1^=6,  F^POINT 

50.00 

94.05 

288.66 

60.00 

83.94 

259.37 

70.00 

73.00 

227.47 

80.00 

62.43 

196.32 

90.00 

52 .82 

167.46 

100.00 

4  4.26 

141.31 

FIGURE  D-3a  EXAMPLE  PERBOIIMSL  OVERTDRMIMO  FORCE 
AMD  MOMENT  GRAPHICS  FILE  FORMAT  CREATED 
WITH  THE  INTERACTION  PLANE  ANALYSIS  OPTION 
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APPENDIX  E 

ROTWASH  PROGRAM  FORTRAN  77  LISTINGS 


ROTWaSH  program  listings  are  presented  in  this  appendix  for  the  ROTWASH 
main  program  and  its  24  subroutines.  The  listings  are  for  a  version  'jf 
the  program  that  is  run  on  IBM  PC/PC-compatible  computers  using  MICROSOFT 
FORTRAN  77,  Version  5.0.  The  tabular  listing  below  indexes  subroutine 
names  and  briefly  describes  functionality  for  user  reference. 


SUBROUTINE 


FUNCTION 


ROTWASH 

CLOUD 

pE-vTX 

FREAD 

GDVTX 


Main  Program  Driver  and  Initialization 
Calculates  Particle  Cloud  Boundaries 

Locates  Disk  Edge  Vortex  System  and  Calculates  Induced 
Velocity  Field 

Prompt/Validate  for  Floating  Point  Input  Data 

Locates  Ground  Vortex  System  and  Calculates  Induced  Velocit 

Field 


HAZARD 

HIPVEL 

HOMCLS 

HSVTX 


HWJVEL 


Driver  Subroutine  for  Hazard  Analysis 
Twin  Rotor  Overturning  Forces  and  Moments 
Home  Cursor  and  Clear  Screen 

Calculates  Induced  Velocity  Field  of  a  Horseshoe  Vortex 
System 

Single  Rotor  Overturning  Forces  and  Moments 


INKEY 

INPUT 

INPUTV 

INPUTX 

ICFNSH 


Menu  Input  Data  Control 

Rotorcraft  Characteristics  Input  Data  Menu 
Velocity  Profile  Status  Menu 
Ground/Disk  vortex  Input  Data  Menu 
Close  Disk  I/O  Files 


lOINIT 

IPVEL 

IREAD 

LEGAL 

LOCATE 


File  I/O  Management  Menu 

Calculates  Interaction  Plane  Velocity  Profile 
Ptompt/Val...date  for  Integer  Input  Data 
Check  Validity  of  Input  Data  Selection  Codes 
Locate  Cursor  Position 


MC.MENT 

PROPRM 

VLINE 

WALJET 

WJVEL 


Calculates  Personnel  Overturning  Forces  and  Moments 
Calculates  Radial  Wall  Jet  Velocity  Profile 
Calculates  Induced  Velocity  from  a  Line  Vortex  Field 
Defines  Initial  Wall  Jet  Position  and  Growth  Parameters 
Calculates  Single  Rotor  Velocity  Profile 


PROGRAM  ROTWASH 


1  C 

2  C 

3  PROGRAM  ROTWASH 

4  C 

5  C  ★*★★★*★★★★★*★**★★★*★*★★★★★♦*★♦★★♦★★**★★*★★★★★*★*****★★**★** 

6  C  ROTORCRAFT  DOWNWASH  HAZARD  ANALYSIS  PROGRAM 

7  C 

8  C  EMA 

9  C  SAMUEL  W.  FERGUSON 

10  C 

11  C  1  APRIL  1993 

12  C 

13  C  PROGRAM  VERSION  2.1 

14  C 

15  C  THIS  PROGRAM  WAS  DEVELOPED  USING  MICROSOFT  FORTRAN  77  FOR 

16  C  THE  DOS  OPERATING  SYSTEM.  CONSOLE  DISPLAY  CONTROL  IS 

17  C  PROVIDED  WITH  THIS  PROGRAM  IN  SEVERAL  SPECIAL  SUBROUTINES. 

18  C  *****:*f*****************’***^lE**************iiC***********ilf*^l:*** 

19  C 

20  PARAMETER (NUM  =  10) 

21  C 

22  CHARACTER*!  OKLIST (NUM) 

23  CHARACTER*!  KEY, KKEY, FLOW, VELHA2,HAZTYP 

24  CHARACTER*!  ICONT(5) 

25  CHARACTER*!  TEMCHAR 

26  CHARACTER* 12  PTSFIL(4) 

27  CHARACTER* 50  COMM (2) 

28  REAL* 4  KE 

29  C 

30  DIMENSION  CONT (9) ,CONTV(7) ,CONTX(8) 

31  C 

32  COMMON  /  CKEY/  KEY,KKEY 

33  COMMON  /CONSTS/  PI, RHO,FPSPKN,RHOD2,DRC 

34  COMMON  /HELGEO/  H, DL, YSEP,WSPD, RADIUS, SHFTAN, DXO 

35  COMMON  /INPUTC/  ICONT, COMM, PTSFIL 

36  COMMON  /INPUTD/  CONT,CONTV,CONTX, YBDLAY 

37  COMMON  /PROFIL/  RJ, ZBJ, 2HJ, ZMJ,UMJ, ZB, ZH, ZM,UM, CU, CY 

38  COMMON  /  UNITS/  lOUl, IOU4, IOU5, IOU6, IOU7, IOU8, IGRAPH 

39  C 

40  DATA  CONT/2. 0,32. 2, 12. 5, 12475. 0,13. 0,37. 5, 0.0, 1.0, 0.0/ 

41  CH53E  DATA  CONT/1 . 0, 0 . 0 , 39 . 5, 56000 . 0, 5 . 0 , 600 . 0 , 0 . 0 , 0 . 975, 0 . 0/ 

42  DATA  YBDLAY/ 1.5/ 

43  DATA  ICONT/ 'V','N','N','W','L'/ 

44  DATA  CONTV/59. 3, 1.0, 10. 0,50. 0,10. 0,100. 0,0.0/ 

45  DATA  CONTX/733. 0,7. 0,50. 0,200. 0,0. 0,20. 0,100. 0,0.0/ 

46  DATA  PTSFIL/'DFRC»Eiai*ipAIP.PTS' , 'OTDFRC.PTS' , 'OTDAIP.PTS' / 

47  DATA  COMM/ 

48  1  'XV-15  CHARACTERISTICS  ARE  USED  AS  INPUT  DATA  • , 

49  2  'GROSS  WEIGHT  MIGHT  BE  ONE  OF  THE  COMMENT  STRINGS  '/ 

50  C 

51  DATA  OKLIST  / ' W' , ' w' , ' I' , ' i' , ' G' , ' g' , ' D' , ' d' , ' X' , ' x' / 

52  C 

53  C  ************************************************************** 

54  C 

55  C  - 

56  C  INITIALIZE  I/O  SYSTEM 

57  C  - 

58  C 

59  CALL  lOINIT 

60  C 

61  C  - 

62  C  INITIALIZE  MISCELLANEOUS  CONSTANTS 

63  C  - 


E-2 


PROGRAM  ROTWASH 


64  C 

65 

66 

67 

68 

69  C 

70 

71  C 

72 

73 

74  C 

75  C 

76  C 

77  C 

78  C 

79 

80  C 

81  C 

82  C 

83  C 

84  C 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94  C 

95 

96  C 

97  C 

98  C 

99  C 

100  C 

101 
102 

103 

104  C 

105  C 

106  C 

107  C 

108  C 

109 

110 
111 
112 

113 

114 

115  C 

116  C 

117  C 

118  C 

119  C 

120  C 

121 

122  C 

123 
12', 

125  C 

126 


PI  -  ACOS(-l.O) 

RHOSL  -  0.0023769 
FPSPKN  -  1.687 
DRC  -  0.01745329252 

50  CONTINUE 

KEY  =  '  ' 

KKEY=  '  ' 


OBTAIN  INPUT  DATA  PARAMETERS  FROM  STATUS  SCREEN 


CALL  INPUT 


DEFINE  PARAMETERS  OBTAINED  FROM  SUBROUTINE  INPUT 


ROTORS  =  CONT(l) 
YYSEP  =  CONT{2) 
RADIUS  “  CONT(3) 
HELGW  -  CONT(4) 
DWNLD  -  CONT(5) 
HAGL  =  CONT(6) 
SHFTAN  =  CONT(7) 
SIGPR  =  CONT(8) 
WSPD  =•  CONT(9) 


VELHAZ  =  ICONT(l) 


ADJUST  GEOMETRY  IF  SHAFT  ANGLE  >0.0  DEGREES 


RSHFT  =  SHFTAN*DRC 
CSHFTA  =  COS (RSHFT) 

DXO  =  HAGL*TAN (RSHFT) 


NON-DIMENS lONALIZE  SOME  OF  THE  INPUT  PARAMETERS 


H  =  HAGL/RADIUa/eSlIWA 

YSEP  =  YYSEP/2 . 0/RADIUS 

EFFGW  =  HELGW* (1.0  +  (DWNLD/100 . 0) ) 

DL  =  EFFGW/ROTORS/PI/RADIUS**2 

RHO  =  SIGPR*RHOSL 

RHOD2  =  0.5*RHO 


SWITCHING  CALLS  HAZARD  PROGRAM  AND  ALLOWS  RETURN  (IF  DESIRED) 
TO  THE  MAINLINE  ROUTINE  TO  CHANGE  ROTORCRAFT  INPUT  PARAMETERS 


IF (VELHAZ .EQ. ' H' ) THEN 

KKEY  =  'H' 

CALL  HAZARD (HAZTYP) 

IF(KEY.EQ.'X' )GOTO  999 
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PROGRAM  ROT WASH 


127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 
181 
182 

183 

184 

185 

186 

187 

188 
189 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 


IF (FLOW. EQ.'X') GOTO  999 
IF(HAZTYP.EQ.'X')G0T0  999 
GOTO  50 

ELSE 

KKEY  -  'V' 

END  IF 


SELECT  FLOWFIELD  OPTION 


10  CONTINUE 


HOME  CURSOR  AND  CLEAR  SCREEN 


ICD  -  0 

CALL  HOMCLS(ICD) 

CALL  LOCATE (3,1) 

WRITE (lOUl, 11) 

11  FORMAT  (  20X, '  SELECT  TYPE  OF  FLOW  TO  BE  ESTIMATED',///, 

1  20X, 'WALL  JET  PROFILE,  TYPE  <W>',/, 

2  20X, ' INTERACTION  PLANE  PROFILE,  TYPE  <I>',/, 

2  20X, 'GROUND  VORTEX,  TYPE  <G>',/, 

2  20X, 'DISK  VORTEX,  TYPE  <D>' , /, 

3  2 OX, 'TO  EXIT  PROGRAM,  TYPE  <X>',//) 


PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT 


40  CONTINUE 

WRITE (lOUl, ' (23X, A, $) ' ) 

1  '  ENTER  DATA  ENTRY  CODE  ==>  ' 

READ (lOUl, ' (Al) ' )  FLOW 

IF  (LEGAL  (FLOW,  IOUl,OKLIST,NUM)  .EQ.DGOTO  40 


MAKE  LEGAL  LOWERCASE  INPUTS  UPPER  CASE  BEFORE  BRANCHING 


IF(FLOW.EQ.'w' ) 
IF (FLOW.EQ. '  i'  ) 
IF (FLOW.EQ. 'g'  ) 
IF (FLOW.EQ. 'd'  ) 
IF  (FLOW.EQ. 'x'  ) 


FLOW  =  'W' 
FLOW  =  'I' 
FLOW  =  'G' 
FLOW  =  'D' 
FLOW  =  'X' 


HOME  CURSOR  AND  CLEAR  SCREEN 


ICD  =  0 

CALL  HOMCLS(ICD) 
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190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 


C  BRANCH  BASED  ON  CHOSEN  OPTION,  ALSO 

C  CHECK  NUMBER  OF  ROTORS  TO  LIMIT  SOME  OPTIONS 

C  - 

C 

IF(FL0W.EQ.'X')G0T0  999 
IF (FLOW. EQ.'W' ) GOTO  12 
C 

IF (FLOW.EQ. ' I' ) THEN 

IF (ROTORS. GT. 1.0) GOTO  12 
WRITEdOUl,'  (////, 14X,A,$)  '  ) 

1  '  REQUIRES  TWO  ROTORS,  TYPE  <RETURN>  TO  CONTINUE  ' 

READ  (lOUl, '  (Al) ' )  TEMCHAR 
GOTO  10 
END  IF 
C 

IF  (FLOW.EQ. ' G' ) THEN 

IF (ROTORS.lt. 2.0) GOTO  1000 
WRITEdOUl,'  (////,  14X, A, $)' ) 

1  '  REQUIRES  ONE  ROTOR,  TYPE  <RETURN>  TO  CONTINUE  ' 

READ  dOUl,  '  (Al)  '  )  TEMCHAR 
GOTO  10 
ENDIF 
C 

IF (FLOW.EQ. ' D' ) THEN 

IF (ROTORS.lt. 2.0) GOTO  1000 
WRITEdOUl,'  (////, 14X,  A, $)') 

1  '  REQUIRES  ONE  ROTOR,  TYPE  <RETURN>  TO  CONTINUE  ' 

READdOUl,  '  (Al)  '  )  TEMCHAR 
GOTO  10 
ENDIF 
C 

GOTO  10 
C 

12  CONTINUE 
C 

C  RADIAL  WALL  JET  FLOWS 

^  ■<r********<ir#f**A***A**# 

c 

c  - 

C  ACCELERATED  SLIPSTREAM  MEAN  VELOCITY 

C  - 

C 

UN  “  SQRT(2 .0*DL/RHO) 

c  - 

C  GROUND  EFFECT  CORRECTION 

C  - - - 

C 

AKG  “  1.0  -  0 . 9*EXP (-2 . 0*H) 

C 

C  - 

C  MEAN  VELOCITY  AT  ROTOR  DISK  (RATIOED  TO  UN) 

C  - 

C 

UB  -  AKG/2 .0 
C 

C  - 

C  FIND  INITIAL  RADIUS  OF  WALL  JET 

C 
C 

C 


CALL  WALJET(H,UB,UN,UMB) 
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253 

254 

255 

256 

257 

258 

259 

260 
261 
262 

263 

264 

265 

266 

267 

268 

269 

270 

271 

272 

273 

274 

275 

276 

277 

278 

279 

280 
281 
282 

283 

284 

285 

286 

287 

288 

289 

290 

291 

292 

293 

294 

295 
298 

297 

298 

299 
3CC 

301 

302 

303 

304 

305 

306 
3C7 
306 
309 
31C 

311 

312 

313 

314 

315 


500  CONTINUE 

iriKEY.EQ.'X' )GOTO  999 
IF(KEir.EQ.'N')GOTO  50 
IF (FLOW.EQ. ' 1' )GOTO  700 

600  CONTINUE 


C 

C 

C 

C 

c 

c 

c 


c 

c 

c 

c 

c 


7CC 


c 

c 

c 

r- 


c 

c 

c 

c 


c 

c 


WALL  JET  REGION 

OBTAIN  INPUT  DATA  FOR  THE  WALL  JET  OPTION 


CALL  INPUTV(FLOW) 

RVZ  -  (CONTVd)  -  DXO) /RADIUS 
DELZ  =  CONTV (2) /RADIUS 
ZMAX  -  CONTV (3) /RADIUS 
BDLAYM  =  YBDLAY/RADIUS 


GENERATE  VELOCITY  PROFILE  AT  RVZ  IN  WALL  JET  REGION 

CALL  WJVEL  (H,  UN,  UMB,  RVZ,  RADIUS,  WSPD,  DELZ,  ZMAX,  DXO,  BDLAYM) 

GOTO  500 

CONTINUE 

INTERACTION  PLANE  UPWASH  DEFLECTION  ZONE 
OBTAIN  INPUT  DATA  FOR  THE  IPLANE  OPTION 

CALL  INPUTV(FLOW) 

XIP  -  (CONTV<l)  -  DXO) /RADIUS 
DELZ  -  CONTV (2) /RADIUS 
ZMAX  -  CONTV (3) /RADIUS 
BDLAYM  =  YBDLAY/RADIUS 

GENERATE  VELOCITY  PROFILE  AT  XIP  IN  INTERACTION  PLANE 


CALL  IPVEL (H, UN, RADIUS, UMB, XIP, YSEP, WSPD, DELZ, ZMAX, DXO, BDLAYM) 

GOTO  500 

HOC  CONTINUE 

IF  (KEY.EQ. 'X'  )GOTO  999 
IFCKEY.EQ. 'N' )GOTO  50 

ICCO  CONTINUE 

HORSESHOE  VORTEX  FLOWS 

************5*  *♦***★★** 
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316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 

339 

340 

341 

342 

343 

344 

345 
348 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 


C 

C 

C 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


c 


OBTAIN  INPUT  DATA  rOR  VORTEX  OPTIONS 


CALL  INPUTX 

OMEGAR  -  CONTXd) 

B  -  CONTX(2) 

VF  -  CONTX(3) 

XT  -  CONTX (4) /RADIUS 

YT  -  CONTX (5) /RADIUS 

DELZ  -  CONTX(6) /RADIUS 

2MAX  -  CONTX (7) /RADIUS 

VF  -  VF*FPSPKN 
AMU  -  VF/ OMEGAR 
CT  -  DL/RH0/0MEGAR**2 


ITERATE  TO  GET  INFLOW  RATIO 

ALOLD  «  SORT {CT/ 2.0) 

DO  1300  ITER-1,100 

ALNEW  «  CT/2 . 0/SQRT(ALOLD**2  +  AMU»*2) 
IF(ABS(ALNEW  -  ALOLD) . LE . 1 . OE-0 5 ) GOTO  1301 
ALOLD  “  ALNEW 
300  CONTINUE 

HOME  CURSOR  AND  CLEAR  SCREEN 


ICD  -  0 

CALL  HOMCLS (ICD) 

CALL  LOCATE (5,1) 

WRITE (lOUl, 20) 

20  rORMAT(  '*****•****•***********•*************',/ 

1  ,' ITERATIONS  EXCEEDED  FOR  INFLOW  RATIO',/ 

2  ^, *»*******•**************************') 


STOP  '  ' 

C 

1301  CONTINUE 
C 

ALAMDA  =  ALNEW 

AMUS  =  AMU/SQRT (CT/2 . 0) 

GAMT  -  OMEGAR*RADIUS*2 .0*PI*CT/B 

CHI  -  AT AN (ALAMDA/ AMU) /2 .0 

C 

C  - 

C  GAMWP  IS  FOR  UNIFORM  LOADING  IN  FORWARD  FLIGHT. 

C  0.625  FACTOR  COMES  FROM  UNPUBLISHED  FAA  FLIGHT  TEST  DATA. 

C  IF  SETTLING  ANGLE  CHI  <-8.0  DEGREES,  THEN  USE  GAMWP  AS  IS. 

C  IF  >  8.0  DEGREES,  THEN  REDUCE  GAMWP  BV  THE  LINEAR 

C  RATE  OF  6.5%  PER  DEGREE  OF  SETTLING  ANGLE.  THE  REDUCTION 

C  IS  A  SIMPLE  APPROXIMATION  FOR  THE  NEAR  TERM  THAT  IS  BASED 
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379  C  ON  AN  ANALYSIS  OF  THE  UNPUBLISHED  FAA  FLIGHT  TEST  DATA. 

380  C  MODS  MADE  FEBRUARY  1993. 

381  C  - 

382  C 

393  KE  -  0.625 

384  GAMWP  -  PI*RADIUS*OMEGAR**2*CT/VT/2.0/KE 

385  C 

386  IF(CHI.LE.0.139616)THEN 

387  GAMW  -  GAMWP 

388  ELSE 

389  GAMW  -  GAMWP*(.’..0  -  (CHI  -  0 . 139616)  *0 . 065*57 . 3) 

390  IF (GAMW.LT.O .0)  GAMW  -  0 . 0 

391  ENDIF 

392  C 

393  HOD  -  H/2.0 

394  C 

395  C  - 

396  C  HOME  CURSOR  AND  CLEAR  SCREEN 

397  C  - 

398  C 

399  ICD  =  0 

400  CALL  HOMCLS(ICD) 

401  CALL  LOCATE (3,  1) 

402  C 

403  IF(FLOW.EQ. 'D' )GOTO  1200 

404  C 

405  C  - 

406  C  GROUND  VORTEX 

407  C  - 

408  C 

409  WRITE (lOUl, 1001)  HOD,AMUS,AMU 

-'lO  ICOl  FORMAT  (  18X,  'ROTOR  HEIGHT  ABOVE  GROUND  H/D  ',2X,F8.4,/ 

411  1  , 18X, 'ADVANCE  RATIO  MU-STAR  ',2X,F8.4,/ 

412  2  , leX, 'ADVANCE  RATIO  MU  ' , 2X, F8 . 4, //) 

413  C 

414  C  - 

415  C  THE  VALUE  INPUT  HERE  REQUIRES  USE  OF  THE  CHART  IN  FIGURE  18 

416  C  OF  THE  ACCOMPANYING  DOCUMENTATION  FOB  THE  GROUND  VORTEX 

417  c  - - - - 

418  C 

419  30  CONTINUE 

420  C 

421  WRITE (lOUl, 31) 

422  31  FORMAT  (  ISX, 'ENTER  GROUND  VORTEX  STRENGTH  RATIO',/, 

423  1  25X, '(SEE  FIGURE  18)  =->  'S) 

424  C 

425  READ (lOUl, *, ERR»30)  GAMG 

426  C 

427  IF(GAMG.LT.O.C)  GAMG  -  0.0 

428  C 

429  C  - . . . 

430  C  CONTINUE  ANALYSIS 

431  C  - 

432  C 

433  GAMG  =  GAMG *G AMT 

434  C 

435  CALL  GDVTX (H, RADIUS, AMU, CT, GAMG, XT, YT,DEL2,ZMAX) 

436  C 

437  GOTO  1100 

438  C 

439  1200  CONTINUE 

440  C 

441  C  - 
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442 

443 

444 

445 

446 

447 

448 

449 

450 

451 

452 

453 

454 

455 

456 

457 


C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 


disk  vortex 

CALL  DEVTX (H, RADIUS, GAMW, CHI , XT, XT, DELZ, ZMAX) 
GOTO  1100 

NORMAL  PROGRAM  EXIT 

999  CONTINUE 
END 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

C 

(2 

C 

C 

c 

c 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


subroutine  cloud (UN, UKB) 

SUBROUTINE  CLOUD 

THIS  SUBROUTINE  MAKES  THE  CALCULATIONS  REQUIRED  IN  ESTIMATING 
THE  PARTICLE  CLOUD  BOUNDARIES  (  NO  DENSITIES  )  FOR  SINGLE  AND 
TWIN  ROTOR  CONFIGURATIONS 

*iir***ik'i)r**i(^ilr**i(^***il^A**«*#4:*«lri^*ilA*4r*******i|*ilr*i^A^^**Aik**ilr****i^i^** 

CHARACTER*!  KEY, KKET 
CHARACTER*!  ICONT(5) 

CHARACTER*12  PTSFIL(4) 

CHARACTER* 50  COMM (2) 

COMMON  /  CKEY/  KEY, KKEY 
COMMON  /CLOUDK/  QSMAX 

COMMON  /CONSTS/  P I , RHO, FPSPKN,  RHOD2,  DRC 

COMMON  /HELGEO/  H , DL, YSEP , WSPD,  RADIUS , SHFTAN, DXO 

COMMON  /INPUTCy  ICONT, COMM, PTSFIL 

COMMON  /PROFIL/  RJ, ZBJ, ZHJ, ZMJ, UMJ, ZB, ZH, ZM, UM, CU, CY 
COMMON  /  UNITS/  lOUl, I0U4, IOU5, IOU6, IOU7, IOU8, IGRAPH 

**ir*tti*tii****(Hi*****»**tiHH,t**«w»*iiHc*it**iHiHlik*k***»»***it******* 


CLEAR  SCREEN  AND  HOME  CURSOR 


ICD  -  0 

CALL  HCMCLS(ICD) 
CALL  LOCATE (3,1) 


READ  IN  THE  TERRAIN  FACTOR  CONSTANT 
(SEE  FIGURE  30  OF  USER'S  GUIDE) 


10  CONTINUE 

WRITE (lOUl, 20) 

20  FORMAT(  15X, 'ENTER  TERRAIN  EROSION  FACTOR  (-ND-) 


■>  ',$) 


READ(IOU1,*,ERR»10)  XKT 


VALIDATE  REAL  INPUT  VALUE 


IF ( XKT , LE. 0. 0. OR. XKT. GT. 500.0) GOTO  1 0 


DEFINE  CLOUD  BOUNDARY  CONSTANTS 


XKT  -  SQRT(XKT) 

QSMX  -  RHOD2* ( (SORT (QSMAX) *UN) **2) 

ERC  -  -0.437 

XUM  -  (UMB*UN)**2 

XCU  =>  CU*CU 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 
89 

O'. 

91 

32 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 
103 
10  4 
105 


Cl  -  1.0 
C2  -2.2 
C 

C  - 

C  SINGLE  ROTOR  CLOUD  BOUNDARY  CALCULATIONS 

C  - 

C 

RCR  »  RADIUS* ( (XKT/ (Cl»RHOD2*XUM*XCU) ) *»ERC) 

RVR  -  0.785*RCR 

ZVR  -  0.329*RCR 

RCVR  -  RCR  -  RVR 

AR  -  (2.0/PI) *ALOC(ZVR/RCVR) 

PHIR  -  (PI/2 .0) *ALOG (RCVR) /ALOG(ZVR/RCVR) 

AXLV  -  AR* ( (-PI/2 .0)  +  PHIR) 

XLV  -  EXP (AXLV) 

HCR  -  XLV  +  ZVR 
C 

c  - 

C  INITIALIZE  INTERACTION  PLANE  BOUNDARIES 

C  - 

c 

RCI  -  0.0 
RVI  -  0.0 
ZVI  0.0 
HCI  -  0.0 

IF(YSEP.LE.0.1)GOTO  30 
C 

C  - - - 

C  INTERACTION  PLANE  CLOUD  BOUNDARY  CALCULATIONS 


RCI  “  RADIUS* ( (XKT/ (C2»RHOD2*XUM*XCU) ) **ERC) 

RVI  -  0.785*RCI 

ZVI  -  0.329*RCI 

RCVR  -  RCI  -  RVI 

AR  -  (2 . 0/PI) *ALOG (ZVI/RCVR) 

PHIR  -  (PI/2 .0) *ALOG (RCVR) /ALOG (ZVI/RCVR) 

AXLV  =  AR* ( (-PI/2 . 0)  +  PHIR) 

XLV  =  EXP (AXLV) 

HCI  =  XLV  +  ZVI 

30  CONTINUE 


10  6  C 

107  C  - 

108  C  PRINTOUT  Of  BOUNDARY  LIMITS 

109  C  - 

110  C 

111  IF(IOU6 .NE. lOUl)  WRITE (IOU6, '('' 1"  )'  ) 

112  C 

113  C  - 

114  C  WRITE  "40  FORMAT"  IF  OUTPUT  TO  GRAPHICS  FILE 

115  C  - 

116  C 

117  IF(IOU6.EQ.6)  WRITE (10U6, 40)  COMM ( 1 ) , COMM ( 2 ) 

118  40  FORMAK  10X,A50,/,  10X,A50,  //) 

119  C 

120  WRITE (IOU6, 50) 

121  50  FC.RMAT(  //, 

122  1  20X,'  SUMMARY  OF  CLOUD  BOUNDARIES',//, 

123  2  20X, '  RC  AND  RV  ARE  FROM  ROTOR  CENTER  (FT) 

124  3  20X,'  RC  ',7X, '  RV  ',7X, '  ZV  ' , 7X, '  HC  ',/) 

125  C 

12  6  WRITE (IOU6, 60)  RCR, RVR,  ZVR,  HCR 


//, 


E-11 


SUBROUTINE  CLOUD 


127 

128  C 

129 

130 

131  C 

132 

133 

134  C 

135  C 

136  C 

137  c 

138  C 

139 

140  C 

141 

142 

143  C 


60  FORMAT(  13X, 'SR',F10.1,3ril.l) 

WRITE (IOU6, 70)  RCI, RVI, ZVI-HCI 
70  FORMAT  (  13X,'IP',Flo.i;3Fn  1) 

WRITE (IOU6, 00)  QSMX 

80  FORMAT(  /,12X,'  QSMAX  -',F7.1,'  PSF',//) 


DECIDE  NEXT  OPTION  WITH  INKEY 


CALL  INKEV 

RETURN 

END 
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1  C 

2  C 

3 

4  C 

5  C 

6  C 

7  C 

8  C 

9  C 

10  C 

11  C 

12 

13 

14 

15 

16 

17  C 

18 

19 

20 
21 
22 

23 

24  C 

25  C 

26  C 

27  C 

28  C 

29  C 

30  C 

31 

32 

33 

34  C 

35 

36 

37 

3S  C 
39  C 
4C  C 

41  C 

42  C 

43  C 

44 

45 

46 

47  C 

48 

49 

A 

51  C 

52  C 

53  C 

54  C 

55  C 

56 

57 

58 

59  C 

60 
61 
62 


SUBROUTINE  DEVTX (H, RADIUS , GAMW, CHI, XT, YT, DELZ , ZMAX) 
SUBROUTINE  DEVTX 

THIS  SUBROUTINE  LOCATES  THE  DISK  EDGE  VORTEX  SYSTEM,  AND 
DIRECTS  THE  CALCULATION  OF  ITS  INDUCED  VELOCITY  FIELD 


CHARACTER*!  TEMCHAR 
CHARACTER*!  KEY,KKEY 
CHARACTER*!  ICONT(5) 
CHARACTER*12  PTSFIL(4) 
CHARACTER*50  COMM (2) 


COMMON 

COMMON 

1 

COMMON 

COMMON 

COMMON 


/  CKEY/  KEY,KKEY 

/CHSVTX/  XLl, YLl, ZL1,XL2, YL2, ZL2, XL3, YL3, ZL3, 
XRl, YR1,ZR1,XR2, YR2, ZR2,XR3,YR3, ZR3 
/CONSTS/  PI,RHO,FPSPKN,RHOD2,DRC 
/INPUTC/  1C0NT,C0MM,PTSFIL 

/  UNITS/  lOUl, IOU4, 10U5, IOU6, IOU7,  IOU8, IGRAPH 


ASSUME  HORSESHOE  SHAPE  -  ASSIGN  LEFT  AND  RIGHT  CORNERS 


XLl  >=  0.0 
YLl  =  -1.0 
ZLl  -  H 

XRl  -  0.0 
YRl  =  1.0 
ZRl  =  H 


SET  UP  DIRECTION  POINTERS  FOR  TFUVILER  ELEMENTS 
POINT  2  IS  AT  GROUND  IMPINGEMENT 


XL2  -  XLl  +  H/TAH(CHI) 
YL2  -  YLl 
ZL2  «  0.0 

XR2  =  XRl  +  H/TAN(CHI) 
YR2  =  YRl 
ZR2  -  0.0 


POINT  THREE  EXTENDS  TRAILER  PARALLEL  TO  GROUND 


XL3  =  XL2  +  1.0 
YL3  =  YL2 
ZL3  =  ZL2 

XR3  “  XR2  +  1.0 
YR3  =  YR2 
ZR3  =  2R2 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 
7  5 

76 

77 

78 

79 

80 
31 
82 

83 

84 

85 

86 
87 
38 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
ICO 
101 
1C2 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


XT  ^  XT*BADIUS 
YT  =  YT*RADIUS 


C 

c 

c 

c 


CLEAR  SCREEN  AND  HOME  CURSOR 


ICD  «  0 

CALL  HOMCLS(ICD) 
CALL  LOCATE (4,  1) 


WRITE  OUTPUT  HEADER 


C 

C 


IF  (IOU6  .NE.  lOUl)  WRITE  (IOU6,  '  ("  1"  )  '  ) 

WRITE(I0U6, 1000)  XT<YT 

1000  FORMAT (  21X,'DISK  VORTEX  VELOCITY  PROFILE  DATA', 

1  ///,14X,'  X-LOCATION  (XT)  =  ' , 2X, F8 . 2 , 2X, ' FT' , 

2  /,14X.'  Y-LOCATION  (YT)  =  ' , 2X, F8 . 2 , 2X , ' FT' , / ) 


C 

C 


GAMW  TO  METRIC  UNITS  OF  METERS  * *2/SEC 

SA:4'.-.'ME  =  GAMW*0 . 0  92903 


5-METERS  »  16.4C42  FEET 


:re  size 


RCD5M  =  0.2*RAX)IUS/16.  4042 

GAMW5M  »  GAWWME*(1.0  -  RCD5M*ATAN ( 1 . 0 /RCD5M) ) 

WRITE  (IOU6,  1001)  GAMW,  GAM'WME,  GAMW  5M 
1001  FORMAT (  15X, 'VORTEX  CIRCULATION  =  ' , FI  0 . 2 , 2X, ' FT  * • 2 / SEC ' , , 

1  15X, ' VORTEX  CIRCULATION  =  ' , FI  0 . 2 , 2X, ' M • » 2 / SEC '  , / 

2  15X,'5-M  INITIAL  CIRCULATION  =  ' , FI  0 . 2 , 2X, ' M * *2 / SEC '  ) 


C 

c 


'  ,  no  .2,  2X,  '  DEG'  ,  /  ///; 


c 

c 

c 


CHID  =  CHI*180 . C/PI 

WRITE (IOU6, 1005)  CKID 
1005  FORMAT(  1 5X,  ' SETTLING  ANGLE 

WRITE (lOUl, ' (23X, A, $) ' ) 

1  '  PRESS  <RETURN>  TO  CONTINUE  ' 

BEAD  (lOUl, '  (Al) ' )  TEMCHAR 

CALL  HOMCLS(ICD) 

IF (IOU6 .EQ. 6)  WRITE(IOU6, 93)  COMM ( 1 ) , CCMM (2 ) 

93  FORMAT(  14X, A50, /, lOX, A50, /.G 

WRITE (IOU6, 1100) 

1100  FORMAT (  12X, 'HEIGHT' , 8X, 'MEAN  VELOCITY' , 5X, ' MEAN  Q',/, 
1  ' 0' , 12X, '  (FT) ' , 8X, '  (FPS) ' , 6X, '  (KN) ' , 5X,  '  (P3F)  ' , /) 


C 

Q 

C 

c 


SET  UP  SWEEP  OF  Z  AT  SPECIFIED  X, Y 

DELZ  AND  ZMAX  COME  FROM  A  MAINLINE  STATUS  MENU 
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127  C  - - - 

128  C 

129  XT  -  XT/RADIUS 

130  VT  -  XT/RADIUS 

131  NPTS  -  IFIX(Z1AX/DELZ)  +  2 

132  LINES  “  0 

133  C 

134  DO  200  I-1,NPTS 

135  C 

136  LINES  -  LINES  +  4 

137  2T  «  (1  -  1) *DELZ 

138  C 

139  CALL  HSVTX (XT, YT, 2T,VXF,VXF,VZF, GAMW, RADIUS) 

140  C 

141  ZZ  -  ZT*RAD1US 

142  VTF  -  SQRT{VXF*»2  +  VYF»*2  +  VZF'^*2) 

143  VXK  -  VXF/FPSPKN 

144  VYK  -  VYF/FPSPKN 

145  VZK  -  VZF/FPSPKN 

146  VTK  -  VTF/FPSPKN 

147  C 

148  QX  -  BHOD2*VXF**2 

149  QY  -  RHOD2*VYF*»2 

150  QZ  -  RH0D2*VZF**2 

151  QT  »  RHOD2*VTF**2 

152  C 

153  C  - 

154  C  KEEP  OUTPUT  PAGE  LENGTH  TO  SIZE  OF  SCREEN 

155  C  - 

156  C 

157  IF  (lOUe.EQ.IOUDTHEN 

158  IF (LINES. LE. 12) GOTO  100 

159  LINES  -  4 

160  CALL  INKEY 

161  IF(KEY.NE.'C' )GOTO  999 

162  WRITE (lOUC, 1100) 

163  ENDIF 

164  C 

165  100  CONTINUE 

166  C 

167  C  - 

168  C  REPORT  X  COMPONENT  OF  VELOCITY 

169  C  - - - 

170  C 

171  WRITE (IOU6,  1101)  ZZ, VXF,VXK,QX 

172  1101  FORMAT(  9X,F8.2,2X,'X',3F10.3) 

173  C 

174  C  - 

175  C  REPORT  Y  COMPONENT  OF  VELOCITY 

176  C  - 

177  C 

178  WRITE  ('''IU6,  1102)  VYF,VYK,QY 

179  1102  FORMAT(  1 9X, ' Y' , 3F10 . 3 ) 

180  C 

181  C  - 

132  C  REPORT  Z  COMPONENT  OF  VELOCITY 

183  C  - 

lo4  C 

185  WRITE (1CU6, 1103)  VZF,VZK,QZ 

186  1103  FORMA-^  (  19X,  '  Z '  ,  3F10 . 3) 

187  C 

188  C  - 

189  C  REPORT  TOTAL  VELOCITY 
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190 

191 
1  92 

193 

194 
1  95 
1  96 
1  99 
1  99 

1  99 

2  00 
201 
2  02 
203 


C  - - 

C 

WRITE ( IOU6, 1 104 )  VTF, VTk,  QT 
11C4  FORMAT!  1 9X,  ' T' , 3F10 . 3 ) 

200  CONTINUE 
C 

CALL  INKEY 
C 

999  CONTINUE 
c 

RETURN 

END 


I  V 


SUBROUTINE  FREAD 


1  C 

2  C 

3  SUBROUTINE  FREAD ( lOUl , PROMPT, VALUE, CONST) 

4  C 

5  Q  *  A  **************  A  *******«**•*<,«  A*  ************  * 

6  C  SUBROUTINE  FREAD  PROMPTS  USER  FOR  A  FLOATING 

7  C  POINT  DATA  ENTRY  AND  CHECKS  VALIDITY  OF  ENTRY 

0  C  ********************************************** 

9  C 

10  PARAMETER (LAST-50) 

11  C 

12  CHARACTER*50  PROMPT, SHOWIT 

13  CHARACTER* 15  ENTRY, BLANK 

14  C 

15  DATA  BLANK  /'  '/ 

16  C 

17  c  ********************************************** 

18  C 

19  C  - 

20  C  PROMPT  USER  FOR  SCALED  FLOATING  POINT  ENTRY. 

21  C  FIND  POSITION  OF  LAST  NON-BLANK  CHARACTER  IN  PROMPT, 

22  C  THEN  STORE  RIGHT  JUSTIFIED  IN  SKCWIT. 

23  c  - - - 

24  C 

25  N  =  LAST  +  1 

26  C 

27  10  IF  (N.EQ.DGOTO  20 

28  C 

29  N  =  N  -  1 
3C  C 

31  IF (PROMPT (N ;N) .EQ.'  ')GOTOlO 

32  C 

33  20  JS  -  LAST  -  N 

34  C 

35  WRITE (SHOWIT, ' (50A1) ' )  ('  ' , J-1/ US) , (PROMPT ( I : I) , I-l , N) 

36  C 

37  c  - 

38  C  NOW  ASK  USER  FOR  DATA  ENTRY 

39  c  - - - 

40  C 

41  30  WRITEdOUl,  '  (/,  1X,A,G13.6) ')  SHOWIT,  VALUE ‘CONST 

42  C 

43  WRITE (lOUl, ' (/, 8X, A, $) ' ) 

44  1  '  ENTER  NEW  VALUE  OR  <RETURN>  TO  LEAVE  AS  IS  — >  ' 

45  C 

46  READ (lOUl, ' (A) ' )  ENTRY 

47  C 

48  IF  ( ENTRY. EQ. BLANK) RETURN 

49  C 

50  READ (ENTRY,  '  (BN, F15 .0)  ' ,ERR-30)  TEMP 

51  C 

52  c  - 

53  C  CONSTANT  CAN  BE  USED  TO  SCALE  OR 

54  C  OR  CONVERT  UNITS  OF  AN  INPUT  VALUE 

55  C  - 

56  C 

57  VALUE  -  TEMP /CONST 

58  C 

59  RETURN 

60  END 

61  C 
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1 

2 

3 

4 

5 

6 
1 
8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
27 
ZP: 

29 

30 

31 

32 

33 

34 

35 

36 

37 
30 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

n 


C 

r* 

c 

c 

c 


c 


c 

c 

c 

c 

c 


c 


c 


SUBROUTINE  GDVTX (H, RADIUS, AMU, CT, GAMG, XT, XT, DELZ, ZMAX) 

*»AAl^*#Aiki^A***********4ri^*A**A****A*ilr**Aflr<llir«kA**A**l^1lilrA*»A*A 

SUBROUTINE  GDVTX 

THIS  SUBROUTINE  LOCATES  THE  GROUND  VORTEX  BASED  ON  THE 
EXPERIMENTS  BY  SUN  AND  CURTIS  (PRINCETON  UNIV. ) ,  AND  THEN 
DIRECTS  THE  CALCULATION  OF  ITS  INDUCED  VELOCITY  FIELD 

THE  OUTPUT  FROM  THIS  SUBROUTINE  SHOULD  BE  USED  CAREFULLY 
FOR  GROSS  ESTIMATION  PURPOSES  ONLY 


CHARACTER*!  ICONT(5) 

CHARACTER*!  TEMCHAR 
CHARACTER*!  KEY,KKEY 
CHARACTER*12  PTSF1L(4) 

CHARACTER*50  COMM (2) 

COMMON  /  CKEY/  KEY, KKEY 

COMMON  /CHSVTX/  Xl 1, YLl , ZLl , XL2 , YL2 , ZL2 , XL3 , YL3 , ZL3 , 

1  XR1,YR1,ZR1,XR2,YR2.  ZR2,XR3, YR3, ZR3 

COMMON  /CONSTS/  PI, RHO, FPSPKN, RHOD2, DRC 
COMMON  /INPUTC/  ICONT, COMM, PTSFIL 

COMMON  /  UNITS/  lOUl , IOU4 , 10U5 , 10U6, :OU7 , IOU8 , IGPAPH 

«r***i*r**4r***ilrik#i»llKtt********A4rilr**AA**4r*AA***<lr***l»lt*W4ri<rA****A4 


HOD 

m 

H/2 .0 

Cl 

m 

1.0  +  1 .2086*HOD**O.4374 

C2 

m 

-0.2786*KOD**0.6757 

ZGV 

m 

-lO.O-AMU  +0.6 

XGV 

m 

-(Cl  +  C2* (AMU/CT) > **2 

XXGV 

« 

XGV*  RADIUS 

2  ZGV 

m 

ZGV*RADIUS 

ASSUME  HORSESHOE  SHAPE  -  ASSIGN  LEFT  AND  RIGHT  CORNERS 


XLl  =■  XGV 
YLl  -  -1.0 
ZLl  ~  ZGV 

XRl  -  XGV 
YRl  -  1.0 
ZRl  -  ZGV 


SET  UP  DIRECTION  POINTERS  FOR  TRAILER  ELEMENTS 


XL2  =  XLl  *  1.0 
YL2  -  YLl 
ZL2  =  ZLl 

XR2  =  XRl  +1.0 
YR2  =  YRl 
2R2  =  ZRl 

XL3  »  XL2  +1.0 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 
87 
80 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 
10  7 
.08 
109 
lie 

ii: 

U2 

113 

114 

115 

116 
117 
116 

119 

120 
121 
122 

123 

124 

125 

126 


YL3  -  YL2 
ZL3  ••  Zli2 
C 

XR3  -  XR2  +  1.0 
YR3  -  YR2 
ZR3  -  ZR2 
C 

c  - 

C  WRITE  OUT  GROUNTD  VORTEX  POSITION/STRENGTH  DATA 

C  - 

C 

IFdOUe.NE.IOUl)  WRITE(IOU6,  '  ("  1"  )  '  ) 

C 

WRITE (1006,1001)  XXGV,ZZGV 

1001  FORMAT (  //, 2 4X, ' GROUND  VORTEX  CORE  POSITION',//, 

1  18X, 'X- LOCATION  (XXGV)  -  ' , IX, F8 . 2 , 2X, ' FT' , / 

2  16X, ' Y-LOCATION  (ZCGV)  =  ' , IX, F8 . 2 , 2X, ' FT' ) 

r* 

WRITE (10U6, 1002)  GAMG 

1002  FORMAT (  /, 13X, ' GROUND  VORTEX  CIRCULATION  -  ' , IX, 

1  F8 .2, 2X, ' FT**2/SEC' , //) 

C 

WRITEdOUl,'  (23X,A,S)') 

1  '  PRESS  <RETURN>  TO  CONTINUE  ' 

C 

READdOUl,  '  (M)  '  )  TEMCHAR 
C 

ICD  ••  0 

CALL  HOMCLSdCD) 

C 

C  - 

C  WRITE  OUTPUT  HEADER 
C  - 

c 

IFdOU6.EQ.6)  WRITE  dOU6,  93)  COMM  ( 1 )  ,  COMM  (2  ) 

93  FORMAT!  14X, A50, /, lOX, A50, //) 

C 

WRITEdOUe,  1100) 

1100  FORMAT(  12X, ' HEIGHT' , 0X, 'MEAN  VELOCITY' , 5X, ' MEAN  Q',/, 

1  ' 0' , 12X, '  (FT) ' , BX, '  (EPS) ' , 6X, '  (KN) ' , 5X, '  (PSF) ' ,  /  ) 

C 


C  SET  UP  SWEEP  OF  Z  AT  SPECIFIED  X,Y 

C  DELZ  AND  ZMAX  COME  FROM  MAINLINE  STATUS  SCREEN 


NPTS  ■=  IFIX  (ZMAX/DELZ)  +  2 
LINES  “  0 
C 

DO  200  1=1, NPTS 
C 

LINES  -  LINES  +  4 
ZT  =  (I  -  1) »DELZ 

r* 


C 


C 


CALL  HSVTX (XT,  YT,  ZT,  VXF, VYF, VZF, GAMG, RADIUS) 
ZZ  “  ZT*RADIUS 

VTF  “  SQRT(VXF**2  +  VYp**2  +  VZF**2) 

VXK  «  VXF/FP3PKN 
VYK  =  VYF/FPSPKN 
VZK  “  VZF/FPGPKN 
VTK  »  VTF/FPSPKN 
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127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 
1  47 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 
161 
162 

163 

164 
1  65 
166 

167 

168 

169 

170 

171 

172 

173 

174 
17  5 
176 
17  7 


QX  -  RHOD2’*'VXF*'*'2 
QY  -  RHOD2*VYF**2 
QZ  -  RHOD2*VZF**2 
QT  -  RH0D2*VTF**2 
C 

IF  (IOU6 .EQ. lOUl) THEN 
IF (LINES. LE. 12) GOTO  100 
LINES  -  4 
CALL  INKBY 

IF(KEY.NE.'C')GOTO  999 
WKiTE ( lOUb, 1100) 

ENDIF 

C 

100  CONTINUE 
C 

C  - 

C  REPORT  X  COMPONENT  OF  VELOCITY 

C  - 

C 

WRlTE(IOU6, 1101)  Z2, VXF, VXK,QX 

1101  FORMAT(  9X,F8.2,2X, 'X',3F10.3) 
C 

C  - 

C  REPORT  Y  COMPONENT  OF  VELOCITY 

C  - 

WRITE ( IOU6, 1102)  VYF,VYK,QY 

1102  FORMAT(  19X, ' Y' , 3F10 . 3) 

C 


C  REPORT  Z  CO^■?ONENT  OF  VELOCITY 

C  - 

c 

WRITE (IOU6, 1103)  VZF,VZK,QZ 

1103  FORMAT(  19X, ' Z' , 3F10 . 3) 

C 

—  _ _ _ _ _ 

C  REPORT  TOTAL  VELOCITY 

C  - 

C 

WRITE (IOU6, 1104)  VTF,VTK,QT 

1104  FORMAT(  19X, ' T' , 3f10 . 3) 

200  CONTINUE 
C 

CALL  INKEY 

c 

999  CONTINUE 


RETURN 

END 


SUBROUTINE  HAZARD 


1  C 

2  C 

3  SUBROUTINE  HAZARD (HAZTYP) 

4  C 

^  Q  ****4r*AA*****AA*A***AA*'***A***A**i<k*Al^**Alkik*A*AA***AA4lr 

6  C  SUBROUTINE  HAZARD  IS  THE  MAINLINE  DRIVER  FOR  THE 
1  C  CALCULATION  OF  SPECIFIC  HAZARDS 

0  *********************  A  ****  <t  ********  4  ****************  * 

9  C 

10  PARAMETER (NUMl  -  6) 

11  PARAMETER (nUM2  -  15) 

12  PARAMETER  (NUM3  -  4) 

13  PARAMETER (NUM4  -  4) 

14  C 

15  CHARACTER*!  OKLSTl (NUMl) 

16  CHARACTER*!  CKLST2(NUM2) 

!7  CHARACTER*!  0KLST3(NUM3) 

!8  CHARACTER*!  OKLST4(NUM4) 

19  C 

20  CHARACTER*!  KEY, KKEY, HAZTYP , HUMTYP 

2!  CHARACTER*!  CHDOL.CVALUE 

22  CHARACTER*!  ICONT(5) 

23  CHARACTER*!2  PTSFIL(4) 

24  CHARACTER*!2  TMPFIL 

25  CHARACTER* 50  COMM (2) 

26  CHARACTER*50  PROMPT 

27  C 

28  DIMENSION  COl  r ( 9 ) , CONTV (7) , CONTX ( 8) 

29  C 

30  COMMON  /  CKEY/  KEY,KKEY 

31  COMMON  /CONSTS/  PI , RHO, FPSPKN, RHOD2, DRC 

32  COMMON  /HELGEO;  H,  DL,  YSEP,  WSPD,  RADIUS,  SHFTAiJ,  DXO 

33  COMMON  /INPUTC/  ICONT, COMM, PTSFIL 

34  COMMON  /INPUTD/  CONT, CONTV, CONTX, YBDLAY 

35  COMMON  /PROFi;,/  RJ,  ZBJ,  ZHJ,  ZMJ,  UMJ,  ZB,  ZH,  ZM,  UM,  CU,  CY 

36  COMMON  /  UNITS/  lOU! , IOU4 ,  IOU5, IOU6,  IOU7 , IOU8 , IGRAPH 

37  C 

30  Q  ***************************************************** 

39  C 

40  DATA  OKLSTl  /  '  C' ,  '  c '  ,  '  M' ,  '  m' ,  '  X' ,  ' -h' / 

41  DATA  OKLST2  /'  ' , ' A' , ' a ' , ' B' , ' b' , ' C' , ' c' , ' D' , ' d' , 

42  1  'E' , 'e' , 'F' ,' f' , 'G' , 'g' / 

43  DATA  OKLST3  / ' W' , ' w' , ' I ' , ' i' / 

44  DATA  OKLST4  /' L' , ' 1' , ' S' , ' s' / 

45  C 

46  0  - 

47  C  CLEAR  SCREEN  AND  HOME  CURSOR 

48  C  - - - 

49  C 

50  ICD  »  0 

51  CALL  HOMCLS(ICD) 

52  CALL  LOCATE (4, 1) 

53  0 

54  C  - 

55  C  DETERMINE  THE  TYPE  OF  HAZARD  ANALYSIS 

56  C  OPTION  THAT  WILL  BE  EXECUTED 

57  C  - 

58  C 

59  WRITEdOUl,  10) 

60  10  FORMAT  (  25X,'  SELECT  TYPE  OF  HAZARD',///, 

61  1  18X, ' OVERTURNING  FORCE/MOMENT,  TYPE  <M>',/, 

62  2  18X, ' PARTICULATE  CLOUDS,  TYPE  <C>',/, 

63  3  18X, 'TO  EXIT  PROGRAM,  TYPE  <X>',//) 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
33 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


C 

C 

C 

c 

c 

11 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

500 


INQUIRE,  OBTAIN,  AND  CHECK  FOR  VALID  MENU  CODE 


CONTINUE 

WRITEdOUl,  '  (23X,  A,  $)  ' )  '  ENTER  HAZARD  CODE  =“> 
READdOUl,  '  (Al)  '  )  HA2TYP 

IF  (LEGAL  (HAZTYP,  IOUl,OKLSTl,NUMl)  .EQ.DGOTO  11 

CORRECT  LOWER  CASE  LETTERS  TO  UPPER  CASE 
TO  USE  AS  VALID  FLAGS  IN  PARENT  SUBROUTINE 


IF(HAZTYP .EQ. ' C' )  HAZTYP-  'C' 
IF(HAZTYP.EQ. 'm' )  HAZTYP  -  'M' 
IF(HAZTYP.EQ.'x' )  HAZTYP  -  'X' 


BRANCH  IF  EXIT  OPTION  CHOSEN 


IF(HAZTYP.EQ. 'M' )GOTO  18 
IF (HAZTYP. EQ. 'C' ) GOTO  18 
IF (HAZTYP. EQ. 'X' ) GOTO  999 

GOTO  11 

CONTINUE 

****★♦★*★*★*******#*★**★********* 

RADIAL  WALL  JET  FLOW  INFORMATION 


ACCELERATED  SLIPSTREAM  MEAN  VELOCITY 

UN  =  SQRT (2 .0*DL/RHO) 

GROUND  EFFECT  CORRECTION 

AKG  =  1.0  -  0. 9*EXP (-2.0*H) 

MEAN  VELOCITY  AT  ROTOR  DISK  (RATIOED  TO  UN) 

UB  =  AKG/2.0 

FIND  INITIAL  RADIUS  OF  WALL  JET 

CALL  WALJET (H,UB,UN,UMB) 

CONTINUE 
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127  C 

128 

129 

130  C 

131  C 

132  C 

133  C 

134  C 

135 

136  C 

137  C 

138  C 

139  C 

140  C 

141  C 

142  C 

143  C 

144  C 

145  C 

146  C 

147  C 

148  C 

149  C 

150  C 

151  C 

152  C 

153 

154  C 

155 

156 

157 

158  C 

159 

160 
161 
162 

163  C 

164  C 

165  C 

166  C 

167  C 

168  C 

169 

17C  C 

171 

172 

173  C 

174 

175  C 

176 

177 

178  C 

179 

180  C 

181 
182 
1B3 

184 

185 

186 
187 
186 

189  C 


IE (KEV.EQ. ' X' )GOTO  999 
IF (KEV.EQ. 'N' )GOTO  999 


BRANCH  IF  CLOUD  OPTION  CHOSEN 


IF (HAZTYP .EQ. 'C' )GOTO  800 
OVERTURNING  FORCES /MOMENTS 


DETERMINE: 

1.  THE  AXIS  ALONG  WHICH  THE  OVERTURNING 
FORCES /MOMENTS  WILL  BE  CALCULATED 

2.  THE  SIZE  OF  THE  PERSON  AFFECTED 

3.  THE  DISTANCES  AX  WHICH  THE  OVERTURNING 
FORCES /MOMENTS  WILL  BE  CALCULATED 


2C  CONTINUE 
ICD  =  0 

CALL  HOMCLS(ICD) 

CALL  LOCATE (2,1) 

WRITE  (lOl’i,  12) 

12  FORMAT(  20X,'  OVERTURNING  FORCE/MOMENT  DATA  MENU',//, 

1  10X,'CODE  PARAMETER  VALUE', 

2  '  UNITS',/) 


PRINT  OUT  MENU  VARIABLES  AS  BASED  ON  THE  WALL  JET 
OPTION  OR  INTERACTION  PLANE  OPTION  SWITCH  SETTING 


IF(ICONT(4) .EQ.'W' ) THEN 


WRITE (lOUl,  14) 


ICONT (4) , ICONT(5) , PTSFIL (3) , 
CONTV(4) ,CONTV(5) ,CONTV(6) , VBDLAY 


ELSE 


WRITE (lOUl, 14)  ICONT(4) , ICONT(5) ,PTSFIL(4) , 

CONTV(4) ,CONTV(3) ,CONTV(6) , YBDLAY 


END  IF 


.4  FORMAT (  11X,'A 


<W>ALL  JET  OR  <I>NTERACTION  PLANE' , 5X, A2 ,/ , 


IIX, ' B 
IIX, 'C 
IIX,  ' 
IIX, ' D 
IIX, ' E 
IIX, ' F 
IIX, ' G 


<L>ABGE  OR  <S>MALL  PERSON 
DATA  OUTPUT  FILENAME 

INITIAL  STATION  POSITION 
HORIZONTAL  INCREMENT 
MAXIMUM  STATION  POSITION 


'  ,  5X,  A2,  /, 

' , 2X, Ai2, /, 

,  /, 

'  , 5X, F7 . 2, 4X,  ' 

'  , 5X, F7 .2, 4X,  ’ .  ,  ,  , 
' , 5X, F7 .2, 4X, ' FT' , / 


MINIMUM  BOUNDARY  LAYER  HEIGHT  ' , 5X , F7 . 2 , 4 X , ' FT ' , / / ) 
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190 

193 

192 

193 

194 

195 

196 

197 

198 

199 
2C0 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 
217 
213 

219 

220 

-)  0  1 
A.  C  X 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 
239 
239 
24C 

241 

242 

243 
24  4 

245 

246 

247 

248 

249 

250 

251 


C 

C 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 


16  CONTINUE 

WRITE (lOUl, ' (8X,A, S) ' ) 

1'  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE 
READdOU*,  '  (Al)  M  CHDOL 

IF(LEGAL(CHDOL,  IOUl,OKLST2,NUM2>  .EQ.DGOTO  16 
DIRECT  OPTIONS  BASED  ON  CHOICE  FOR  "CHDOL" 


* 


IF  (CHDOL. EQ, 


)GOTO  30 


CHOOSE  WALJET  OR  IPLANE  OPTION 

IF  (CHDOL. EQ. 'A'  . OR . CHDOL . EQ . ' a' ) THEN 
40  CONTINUE 


WRITE (lOUl, ' (/, 35X, A, 1X,A2/) ' )  '  ANALYSIS  TYPE  -  ',ICONT(4) 
WRITE (lOUl, ' (37X, A, S) ' )  '  ENTER  NEW  CODE  — >  ' 

BEAD (lOUl, ' (Al) ' )  CVALUE 

IF  (LEGAL  (CVALUE,  lOUl,  OKLST3,NOM3)  .EQ.  DGOTO  40 
ICONT(4)  -  CVALUE 

IF(ICONT(4) .EQ.'w' )  ICONT(4)  -  ' W' 

IF(ICONT(4) .EQ.'i')  ICONT(4)  -  'I' 

GOTO  20 

ENDIF 


CHOOSE  LARGE  OR  SMALL  PERSON 

IF (CHDOL. EQ. ' B' . OR . CHDOL . EQ .' b' ) THEN 
41  CONTINUE 

WRITEdOUl,  '  (/,  35X,A,  1X,A2/) ')  '  PERSON  TYPE  -  MCONT(5) 
WBITEdOUl,  '  (37X,A,  $) ')  '  ENTER  NEW  CODE  »«>  ' 

READ (lOUl, ' (Al) ' )  CVALUE 

IF(LEGAL(CVALUE,  I0U1,0KLST4,NUM4)  .EQ.DGOTO  41 
ICONT(5)  -  CVALUE 

IF(ICONT(5) .EQ. '1' )  ICONT(5)  -  ' L' 

IF(ICONT(5) .EQ. 's' )  ICONT(5/  =  'S' 

GOTO  20 

ENDIF 

CHOOSE  GRAPHICS  FILENAME 
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253  C 

254  C 

255 

256  C 

257 

258 

259 

260 
261 
262 

263 

264  C 

265 

266 

267  C 
26B 

269  C 

270 

271 

272 

273  C 

274 

275  C 

276  C 

277  C 

278  C 

279  C 
28C 

281  C 

282 

283 

284  C 

285 

286 

287  C 

280 

289  C 

290  C 

291  C 

292  C 

293  C 

294 

295  C 

296 

297 

298  C 

299 

300 

301  C 

302 

303  C 

304  C 

305  C 

306  C 

307  C 

308 

309  C 

310 

311 

■312  C 

313 

314 

315  C 


IF (CHDOL . EQ . ' C' . OR . CHDOL . EQ . ' c ' ) THEN 

IF(ICONT(4) .EQ.'W' ) THEN 

WRITE (lOUl,' (/, 25X, A, 1X.A12/) ' ) 

1  '  FILENAllE  -  ',PTSFIL(3) 

ELSE 

WRITEdOUl,'  (/,25X,A,  1X,A12/)' ) 

1  '  FILENAME  -  ',PTSFIL(4) 

END  IF 

WRITEdOUl,  '  (20X,  A,  $)  ' ) 

1  '  ENTER  NEW  FILENAME  (xxxxxxxx . xxx)  — >  ' 

READdOUl,  '  (A12)  '  )  TMPFIL 

IF(ICONT(4) .EQ.'W')  PTSFIL(3)  -  TMPFIL 
IF(ICONT(4) .EQ.' I')  PTSFIL(4)  -TMPFIL 
GOTO  20 

END  IF 


CHOOSE  INITIAL  STATION  POSITION 


IF  (CHDOL. EQ. 'D'  . OR . CHDOL . EQ .' d' ) THEN 

PROMPT  -  ' INITIAL  STATION  POSITION  - 
CALL  FREAD  dOUl, PROMPT, CONTV ( 4 ) ,1.0) 

IF(CONTV(4) .LT.O .0)  CONTV{4)  -0.0 
GOTO  20 

ENDIF 


CHOOSE  HORIZONTAL  INCREMENT 


IF  (CHDOL . EQ . ' E'  . OB . CHDOL . EQ . ' e' ) THEN 

PROMPT  -  'HORIZONTAL  INCREMENT  -  ' 
CALL  FREAD (lOUl, PROMPT, CONTV (5) ,1.0) 

IF (CONTV(5) .LT.O . 0)  CONTV(5)  -  0.0 
GOTO  20 

ENDIF 


CHOOSE  HORIZONTAL  INCREMENT 


IF (CHDOL . EQ . ' F' . OR . CHDOL . EQ . ' f ' ) THEN 

PROMPT  -  'MAXIMUM  STATION  POSITION  =  ' 

CALL  FREAD (lOUl, PROMPT, CONTV (6) ,1.0) 

IF (CONTV(6) .LT.CONTV(4) )  CONTV(6)  »CONTV(4) 
GOTO  20 
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316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 

334 

335 

336 

337 

338 
33  9 

340 

341 

342 

343 

344 

345 

346 

347 

348 

349 

350 

351 

352 

353 

354 

355 

356 

357 

358 

359 

360 

361 

362 

363 

364 

365 

366 

367 

368 

369 

370 

371 

372 

373 

374 

375 

376 

377 

378 


ENDIF 

C 

IF  (CHDOL . EQ . ' G'  . OR . CHDOL . EQ . ' g' >  THEN 
C 

PROMPT  -  'MINIMUM  BOUNDARY  LAYER  HEIGHT  -  ' 

CALL  FREAD (lOUl, PROMPT, YBDLAY, 1.0) 

C 

IF(YBDLAY.LT.O.O)  YBDLAY  -  0.0 
GOTO  20 
C 

ENDIF 

C 

GOTO  20 
C 

30  CONTINUE 
C 

ICD  -  0 

CALL  HOMCLS(ICD) 

C 

IF(ICONT(4) .EQ.'I')GOTO  700 
C 

600  CONTINUE 
C 

Q  - 

C  WALL  JET  REGION 

C 

C  OBTAIN  DATA  FOR  THE  HWJVEL  OPTION 

C  - 

c 

RVZ  “  (CONTV(4)  -  DXO) /RADIUS 

DELH  -  CO.NTV(5) 

HMAX  -  CONTV ( 6 ) 

HUMTYP  -  ICONT(5) 

BDLAYM  -  YBDLAY/RADIUS 

C 

C  - 

C  GENERATE  VELOCITY  PROFILE  AT  RVZ  IN  WALL  JET  REGION 

C  - 

C 

CALL 

HWJVEL (H, UN, UMB, RVZ , RADIUS , WSPD, DELH, HMAX, HUMTYP , DXO, BDLAYM) 
C 

GOTO  500 
C 

700  CONTINUE 
C 

C  - 

C  INTERACTION  PLANE  UPWASH  DEFLECTION  ZONE 

C 

C  OBTAIN  INPUT  DATA  FOR  THE  HIPVEL  OPTION 

C  - 

C 

XIP  =  {CONTV(4)  -  DXO) /RADIUS 
DELH  -  CONTV (5) 

HMAX  -  CONTV (6) 

HUMTYP  =  ICONT(5) 

BDLAYM  =  YBDLAY/RADIUS 
C 

C  - 

C  GENERATE  VELOCITY  PROFILE  AT  XIP  IN  INTERACTION  PLANE 

C  - 

C 

CALL  HIPVEL (H, UN, RADIUS , UMB, XIP, YSEP, WSPD, DELH, HMAX, 
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379 

380 

381 

382 

383 

384 

385 

386 

387 

388 

389 

390 

391 

392 

393 

394 

395 

396 

397 

398 

399 
•100 
401 


HUMTYP, DXO, BDLAYM) 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GOTO  500 
800  CONTINUE 

CAIrCUIATE  PARTICULATE  CLOUD  BOUNDARIES 

CALL  CLOUD (UN, UMB) 

GOTO  500 

NORMAL  PROGRAM  EXIT 
999  CONTINUE 

return 

END 


C 
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1 

2 

3 

A 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

4  6 

47 

48 

49 

50 

51 

52 

53 

54 

55 

5  6 

57 

58 

59 

60 
61 
62 
63 


C 

C 


c 

c 

*  Ik 

c 

c 

c 

c 

*  * 

c 


SUBROUTINE  HIPVEL (H, UN, RADIUS, UMB, XIP, YSEP, WSPD, DELH, 
►  HMAX,HUMTYP,DXO,BDLAyM) 


******tk*ik*******<k*****«********D*lt**lk***tk**1k****lkik**)k*lk*****lk* 

SUBROUTINE  HIPVEL  GENERATES  THE  VELOCITY  PROFILE  AND  THE  FORCES 
AND  0''>/ERTURNING  MOMENTS  FOR  A  HUMAN  BEING  ALONG  THE  INTERACTION 
PLANE  FOR  THE  TWIN  ROTOR  CASE 


CHARACTER*!  TEMCHAR 
CHARACTER*!  KEY, KKEY, HUMTYP 
CHARACTER*!  ICONT(5) 

CHARACTER*12  PTSFIL(4) 

CHARACTER* 50  COMM (2) 

COMMON  /  CKEY/  KEY, KKEY 
COMMON  /CONSTS/  PI, RHO, FPSPKN, RHOD2, DRC 
COMMON  /INPUTC/  ICONT, COMM, PTSFIL 
COMMON  /PERSON/  QP(12),DSET 

COMMON  /PROFIL/  RJ, ZBJ, 2HJ, ZMO, UMJ, ZB,  ZH, ZM, UM, CU, CY 
COMMON  /  UNITS/  lOUl , IOU4 , IOU5. IOU6,  10U7 , IOU8 , IGRAPH 


C 

C 

■*  * 

c 


C 

c 

C 

C 

c 


**************■>*************•****•***************1 


r********** 


rCD  »  0 

CALL  HOMCLS(ICD) 


INPUT  FOR  DELH  AND  HMAX  COMES  FROM  INPUTV  STATUS  MENU 


C 

C 

C 

C 

C 


1 


c 

C 


DSET  -  DELH 

IF(DSET.EQ.O. )DELH  -  HMAX 

DF.LH  -  DELH/RADIUS 
HMAX  =  HMAX/RADIUS 

NHPTS  =  IFIX((HMAX  -  XIP) /DELH)  +  1 
IF (DSET. EQ. 0. ) GOTO  33 


WRITE  OUTPUT  HEADER  (FOB  PLOT  FILE,  SEE  BELOW) 


IF(IOU6  .NE.  lOUl)  WRITE  (IOU6,  '  ("1"  ) '  ) 

IF (IOU6 .EQ. 6)  WRITE (IOU6, 93)  COMM ( 1) ,  COMM  2 ) 

93  FORMAT(  lOX, A50 , / , lOX, A50, / / ) 

WBITEdOUe,  1001) 

001  FORMAT(  12X,'  SUMMARY  OF  OVERTURNING  FORCES  AND  MOMENTS',//, 

1  19X, 'RADIUS' , 6X, 'TOTF'  ,  6X,  'TOTM'  ,  /  , 

2  20X, '  (FT) ' , 7X, '  (LB) ' , 5X, '  (FT-LB) ' , /) 

33  CONTINUE 


WRITE  OUT  GRAPHICS  FILES  IF  SWITCH  IS  SET  BY  USER 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 
97 
88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 
'13 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


C 

C 

C 

C 

C 

c 

c 


c 


c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


IFdGRAPH.EQ.DTHEN 


OPEN  GRAPHICS  FILE  AND  WRITE  FILE  HEADER 


OPEN (I0U8,FILE-PTSF1L (4) /STATUS-' NEW' ,ERR-2000) 

WRITE (IOU8, 83)  COMM(l) ,COMM(2) 

83  FORMAT(  10X,A50, /,  10X.A50, //) 

WPITEdOUa,  80) 

80  FORMAT (  IX, ' TITLE-"TW1N  ROTOR  DAIP  DATA"') 

WRITE(I0U8, 81) 

81  FORMAT (  IX, 'VARIABLES  -  DAIP, TOTF,  TOTM' ) 

WRITE (IOU8, 88) 

88  FORMAT (  IX, ' ZONE  T  -  "GW  -  xxxxx  LB,  WAGL  -  xx  FT",', 
*  '  I-X,  F-POINT') 


ENDIF 


BEGIN  LOOP  INCREMENTING  THE  RADIAL  POINTS  AT  WHICH 
THE  OVERTURNING  MOMENT  CALCULATIONS  ARE  MADE 


DO  565  K  -  1,NHPTS 


TF  IS  INTERACTION  PLANE  AMPLIFICATION  FACTOR 
(SEE  NOTE  IN  IPVEL.FOR  FOR  VERSION  2.1) 


TF  -  1,65  -  (0.65) *EXP(-0.5*XIP) 


GET  PARAMETERS  AT  BASE  RADIUS  FOR  'BOUNDARY  LAYER' 


RIPO  -  SQRT(XIP**2  +  YSEP**2) 


' PROPRM'  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 
OF  A  RADIAL  WALL  JET  (WITHOUT  INTERACTION  PLANE) 


CALL  PROPRM(H,UMB,RIP0) 

ZIPB  -  ZB 
ZIPM  -  ZM 
ZIPH  -  ZH 

RIPM  =■  SQRT(XIP**2  +  (YSEP  +  ZIPM)**2) 


CALL  PROPRM (H,UMB, RIPM) 


UMM  =  UM 


OUTPUT  HEADER 
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127  C  - 

126  C 

129  IF{DSET.NB.0.)GOTO  78 

130  C 

131  XXIP  -  RADIUS*XiP 

132  XIPOUT  -  XXIP  +  0^0 
333  C 

134  IFdOUe.NE.IOUl)  WRITE(I0U6/' ("1")') 

135  C 

136  IF(10U6.EQ. 6)  WRITE ( I0U6, 93)  CONM(l) , COMM (2) 

137  C 

138  WRITE (10U6, 1000)  XIPOUT 

139  1000  FOBMAT(  2X, 'TWIN  ROTOR  INTERACTION  PLANE  VELOCITY  PROFILE', 

140  1  '  AT  DISTANCE  -  ',r7.1,'  FT*,//) 

141  C 

142  WRITE  (IOU6,  1002) 

143  1002  FOFMATC  3X, 'HEIGHT' , 6X, 'MEAN  VELOCITY' , 7X, ' PEAK  VELOCITY' , 6X, 

144  1  'MEAN  Q',4X,'PEAK  Q',/, 

145  2  3X,'  <FT)',7X,' (FPS)',6X,' <KN)',5X,' (FPS)',6X,' (KN)',5X, 

14€  3  '  (PSF)' ,5X,' (PSF)',/) 

147  78  CONTINUE 

148  C 

149  C  - 

150  C  'AN'  IS  ACTUALLY  '-  1. 0/7.0' 

151  C  - 

152  C 

153  Ai;  -  0.142857142 

154  C 

155  DELZ  -  0.5/RADIUS 

156  NPTS  -  12 

157  C 

158  UO  500  I  -  1,NPTS 
J.59  C 

160  ZIP  -  DELZ*  (I  -  1)  +  (0.25/RAX)IUS) 

161  C 

162  C  - 

163  C  GET  MAX  WALL  JET  VELOCITY  AT  EFFECTIVE  RADIOS 

164  c  - 

165  C 

166  RIP  -  SQRT(XIP**2  +  (YSEP  +  ZIP)**2) 

167  C 

168  CALL  PROPRM(H,UMB,RIP) 

169  C 

170  VN  -  UN 

171  VZ  -  UM 

172  C 

173  C  - 

174  C  INTERACTION  PLANE  'BOUNDARY  LAYER' 

175  C 

176  C  CODE  MODIFIED  IN  MAY  1992  FOR  USER  SPECIFIED 

177  C  MINIMUM  BOUNDARY  LAYER  THICKNESS  (BDLAYK) 

178  C  - 

179  C 

180  ZIPl  -  BDLAYM 

181  C 

182  IF (ZIP. LT.21PM.OR.ZIP.lt. ZIPl) THEN 

183  C 

184  IF(ZIPi.LT.ZIPM)THEN 

185  C 

106  VZ  -  UMM* (ZIP/ZIPM) **AN 

187  C 

188  ELSE 

189  C 
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190  V2  -  UMM* (ZIP/ZIPl)  **AN 

191  C 

192  ENDIF 

193  C 

194  ENDIF 

195  C 

196  C  - 

197  C  DEVELOPED  INTERACTION  PLANE  JET 

198  C  - 

199  C 

200  VH  -  TF*VZ*X1P/RIP 

201  W  =  TF*V2*  (YSEP  +  ZIP) /RIP 

202  C 

203  Z2  “  ZIP*RADIUS 

204  C 

205  C  - 

206  C  MEAN  HORIZONTAL  VELOCITIES  AND  DYNAMIC  PRESSURE 

207  C  - 

208  C 

209  VHMF  -  VH*UN 

210  VHMK  =  VHMF/FPSPKN 

211  C 

2:2  C  - 

213  C  PEAK  VELOCITIES  (BOTH  FT’/SEC  AND  KNOTS) 

214  C 

215  C  EQUATION  FOR  VMFD3I  UPDATED  FROM  Ist  TO 

216  C  2nd  ORDER  POLYNOMIAL  FOR  VERSION  2.1 

217  C  - - - - 

218  C 

219  VMFD3I  =  C. 712887  +  0.304369’XIP  -  0  . 018496*XIP*X:P 

220  C 

221  IF(VMFD3I .LT. 1 .2)  VMFD3I  -  1.2 

222  C 

223  VHPF  -  VMFD3I*VHMF 

224  VHPK  »  VHPF/FPSPKN 

225  C 

226  C  . - - - - - . . 

227  0  THE  EFFECT  OF  WIND  IS  TO  ADD  (DOWNWIXD  SIDE)  OB  SUBTRACT 

228  C  (UPWIND  SIDE)  ' XWK'  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 

229  C  THE  HORIZONTAL  FROFILE  VELOCITY  (EMPIRICAL,  CH-53E  BASED) 

230  C  - - - - - - 

231  C 

232  XKW  =  (-0 . 5*H)  +  2.5 

233  C 

234  IF  (X. W , LT . 1 . 0) XKW  -  1.0 

235  C 

236  WSPD2  '  WSP'J*XKW 

237  VHMK  =  VHMK  +  W0PD2 

238  VHMF  =  VHMK*FPSIKN 

239  VHPK  =  VHPK  +  WS?d2 

240  VHPF  =  VHPK*FPSPKN 

241  C 

242  C  - 

243  C  DYN/vMIC  PRESSURE 

244  C  - 

245  C 

246  QHK  =  RH0D2*VHMF* *2 

247  QP(I)  =  RHCD2*VHPF**2 

248  C 

249  IF(DSET.NE.O. )GOTO  77 

250  C 

251  C  - - 

252  C  BEPOPT  HORIZONTAL  COMPONENTS 
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253  C  - 

254  C 

255  WRITE (IOU6, 1003)  ZZ, VHMF, VHMK, VHPF, VHPK, QHM, QP ( I) 

256  1003  FORMAT  (  F8 . 2, 2X, 6F10 . 3) 

257  77  CONTINUE 

258  C 

259  500  CONTINUE 

260  C 

261  IF (DSET.NE. 0 . )GOTO  520 

262  WRITE (lOUl, 73) 

2  63  73  FORMAT  (  ) 

264  C 

265  WRITE  aOUl,  '  (19X,A,  $)  '  ) 

266  1  '  TTPE  <RETURN>  TO  CONTINUE 

267  C 

2  68  READdOUl,  '  (Al)  '  )  TEMCHAR 

269  C 

270  ICD  -  0 

271  CAT.L  HOMCLS(ICD) 

272  IF  (IOU6.NE.IOU1)  WRITE  ( IOU6, '("  1 ")' ) 

273  C 

274  IF(IOU6.EQ.6)  WRITE ( IOU6, 93)  COMM ( 1) , COMM (2 ) 

275  C 

276  WRITE (IOU6, 1007)  XIPOUT 

277  1007  FORMAT (  12X, ' TWIN  ROTOR  FORCE  PROFILE  AT  DISTANCE  - 

278  1  F7.1,'  FT',//) 

279  WRITE (IOU6,  1008) 

280  1008  FORMAT(  2X, ' HEIGHT' , 6X, ' PEAK  Q' , 6X, ' FOVER' , 7X, ' OVERM' , 7X, 

281  1  'TOT  F',7X, 'TOT  M',/, 

282  2  3X, '  (FT) ' , 8X, '  (PSF) ' , 7X, '  (LB) ' , 6X, '  (FT-LB) ' , 7X, 

283  3  ' (LB) ', 6X, ' (FT-LB) ', /) 

284  520  CONTINUE 

285  C 

286  C  - - - 

237  C  CALL  SUBROUTINE  TO  CALCULATE  THE 

288  C  FORCES  AND  MOMENTS  ON  A  HUMAN  BEING 

289  C  - - - 

290  C 

291  CALL  MOMENT (NPTS, HUMTYP, TOTF, TOTM) 

292  C 

293  IF(USET.EQ.0.)GOTO  545 

294  C 

295  HH  -  X1P*RADIUS 

29  6  HHCUT  =  HH  +  DX'^ 

297  C 

298  WRITE (ICU6, 1014)  HHOUT, TOTF, TOTM 

299  1014  FORMAT(  1 CX, F0 . 1 , 2F1 C . 3) 

3C0  C 

301  IF (IGRAPH . EQ. 1) THEN 

302  C 

303  WRITE (IOU8, 90)  HHOUT, TOTF, TOTM 

304  90  FORMAT(  IX, F7 . 2 , IX, F7 . 2, IX, FS . 2 ) 

305  C 

506  ENDIF 

307  C 

308  545  CONTINUE 

309  C 

310  XIP  =  XIP  +  DELH 

311  C 

312  565  CONTINUE 

313  C 

314  C  - 

3:5  C  CLOSE  AN  OPEN  GRAPHICS  FILE 
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316  C  - 

317  C 

318  ir (IGRAPH.EQ.l)THEN 

319  C 

320  CLOSE (IOU8, STATUS-' KEEP') 

321  C 

322  ENDir 

323  C 

324  CALL  INKEY 

325  C 

326  GOTO  999 

327  C 

328  C  - 

329  C  THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 

330  C  OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 

331  C  - 

332  C 

333  2000  CONTINUE 

334  C 

335  CALL  HOMCLS(O) 

336  WRITB(IOU1,2001) 

337  2001  FORMAT(  ////,0X, 

338  1  '  ***  ERROR  ***  PLEASE  CHOOSE  A  NEW  OUTPUT  FILENAME', 

339  2  ///,0X, '  TYPE  <RETURN>  TO  CONTINUE  ' , $) 

340  READ (lOUl, ' (Al) ' )  TEMCHAR 

341  KEY  -  'P' 

342  C 

343  999  CONTINUE 

344  C 

345  RETURN 

346  END 

347  C 

348 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 
2.5 
26 

27 

28 
29 
3C 

31 

32 

33 

34 

35 

36 

37 

38 

3  9 

40 

41 

42 

43 

44 

45 

46 

4  7 

48 

49 

50 

51 

52 

53 

54 

55 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 


c 

c 

c 


SUBROUTINE  HOMCLS (CODE) 

SUBRCUTINK  HOMCLS 

THIS  SUBROUTINE  HOMES  THE  CURSOR  AND  CLEARS  THE  TER  IINAL 
SCREEN  (CODE-0)  OR  HOMES  THE  CURSOR  ONLY  (CODE-1) 

****!»  *r*»rt**<lr***  A  ★*>**★***  Jttlr*l^*****^*#**/'iri>**!*r*-*^**  'A**** 

.COMMON  /  UNITS/  lOUl,  IOU4,  IOU5,  IOU6,  1CU7,  1CU8 ,  IGRiVPH 

INTEGER *4  CODE 

CHARACTER*4  ED 
CHARACTER*!  EED{4) 

EQUIVALENCE  (ED, EED  (1) ) 

CHARACTER*  3  EF, 

C.HAKACTER*!  EFEO) 
fclQUIVALENCE  (EE, FEE (1) ) 

IF  (CODE  .EQ, .',)  GOTO  2  0 


HOKE  CURSOR  /iND  CLEAR 
ANSI  CONTROL  SEQUENCE; 

SCREEN 

;  ED  -  ESC(2J 

EED(l)  -  CHAR(27) 

EEC  (2)  -  CHAR(91» 

FED (3)  -  CHAR  (50) 
SED(4)  «  CHAR (74) 

WRITE (lOUl, *)  ED 

CONTINUE 

HOME  CURSOR  ONLY 

ANSI  CONTROL  SEQUENCE; 

;  EE  -  ESCi'H 

EEE(l)  -  CHAR(27) 
EEL  (2)  -  CHAR  (91) 
EEE(3)  -  CHAR (72) 

WRITE (lOUl, *)  EE 

RETURN 

END 


I 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
Ifc 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 


62 

63 


C 

C 

SUBROUTINE  HSVTX (XT,  YT,  ZT,  VX,  VY,  VZ, GAMMA, RADIUS) 
C 

c 


c 

c 

c 

c 


SUBROUTINE  HSVTX 

THIS  SUBROUTINE  DIRECTS  THE  CALCULATION  OF  THE  INDUCED  VELOCITY 
FIELD  DUE  TO  A  HORSESHOE  VORTEX  SYSTEM  OF  UNIT  STRENGTH.  POINT 


C  (LEFT  =  LI,  RIGHT  =  Rl )  DEFINE  THE  EXTENT  OF  THE  BOUND  PORTION 

OF 

C  THE  HORSESHOE.  THE  TRAILERS  START  AT  POINT  1  AND  EXTEND 

THROUGH 

C  POINT  2,  AND  THEN  ON  TO  POINT  3.  THIS  ALLOWS  TWO  ELEMENTS  FOR 

C  EACH  TRAILER  SO  THAT  IT  CAN  'BEND'  TO  ACCOUNT  FOR  GROUND 

CONTACT . 

C 

*ifc'*<rA*******^*i*t1r'***»*<^******r***********<>*i<t*****Tk****H****yk*****1»t^* 

/-* 

COMMON  /CHSVTX/  XLl ,  YLl ,  2L1 ,  y.L2 ,  YL2 ,  ZL2 ,  XL3 ,  YL3 ,  ZL3 , 

1  XRl , YRl , ZRl , XR2 , YR2 , ZR2 ,  XR3 ,  YR3 ,  ZR3 

COMMON  /CVLINE/  IF  I , XA, YA, ZA, XB, YB, ZB, XC, YC, 2C, Q1 , Q2 , Q3 


C 

c 


c 

C 


AT  SPECIFIED  (X,Y,Z)  TARGET  POINT  IN  VICINITY 
OF  ROTOR,  CALCULATE  THE  VECTOR  VELOCITY 


VX  =  0.0 
VY  =  0,0 
VZ  =  0.0 

XC  =  XT 
YC  =  YT 
ZC  =  ZT 


LEFT  TRAILER  CONTRIBUTION,  POINT  1  TO  POINT  2 


46 

IFI  ^ 

=  0 

47 

XA  - 

XLl 

48 

YA  = 

YLl 

49 

ZA  = 

ZLl 

50 

XB  = 

XL2 

51 

YB  = 

YL2 

52 

ZB  = 

ZL2 

53 

C 

54 

CALL 

VLINE 

55 

C 

56 

VX  = 

VX  -  Q1 

57 

VY  = 

VY  -  Q2 

58 

VZ  = 

VZ  -  Q3 

5  9 

Vs," 

D  0 

61 

V. 

LEFT 

TRAILER  IMAGE 

C 

C 
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64 

ZA  -  -ZA 

65 

ZB  -  -ZB 

66 

C 

67 

CALL  VLINE 

68 

C 

69 

VX  -  VX  +  Q1 

70 

VY  -  VY  +  Q2 

71 

VZ  -  VZ  +  Q3 

72 

c 

73 

c 

74 

c 

LEFT  TRAILER  CONTRIBUTION,  POINT  2  TO  POINT  3 

75 

c 

76 

c 

77 

IFI  -  1 

78 

XA  -  XL2 

79 

YA  -  YL2 

80 

ZA  -  ZL2 

81 

XB  -  XL3 

82 

YB  -  YL3 

83 

ZB  -  ZL3 

84 

c 

85 

CALL  VLINE 

86 

c 

87 

VX  -  VX  -  Q1 

88 

VY  -  VY  -  Q2 

89 

VZ  -  VZ  -  Q3 

90 

c 

91 

c 

92 

c 

LEFT  TRAILER  Il-IAGE 

93 

c 

94 

c 

9  5 

ZA  -  -ZA 

96 

ZB  -  -ZB 

97 

c 

98 

CALL  VLINE 

99 

c 

100 

VX  -  VX  +  Q1 

101 

VY  -  VY  +  Q2 

102 

VZ  -  VZ  +  Q3 

103 

c 

104 

c 

105 

c 

SPANWISE  VORTEX  CONTRIBUTION 

106 

c 

107 

c 

108 

IFI  »  0 

109 

XA  -  XLl 

no 

YA  -  YLl 

111 

ZA  -  ZLl 

112 

XB  -  XRl 

113 

YB  -  YRl 

114 

ZB  -  ZRl 

115 

c 

116 

CALL  VLINE 

117 

c 

118 

VX  =  VX  +  Q1 

119 

VY  -  VY  +  Q2 

120 

VZ  =  VZ  +  Q3 

121 

c 

122 

c 

123 

c 

SPANWISE  VORTEX  IMAGE 

124 

c 

125 

c 

126 

ZA  =  -ZA 
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12T 

128  C 

129 

130  C 

131 

132 

133 

134  C 

135  C 
13  6  C 

137  C 

138  C 

139 

140 

141 

142 

143 

144 

145 

146  C 

147 

148  C 

149 

150 

151 

152  C 

153  C 

154  C 

155  C 
15>'  C 

157 

158 

159  C 

160 

161  C 

162 

163 

164 

165  C 

166  C 

167  C 

168  C 

169  C 

170 

171 

172 

173 

174 

175 

176 

177  C 

178 

179  C 

180 
181 
182 

183  C 

184  C 

185  C 
•186  C 

187  C 

188 
189 


ZB  -  -ZB 

CALL  VLINE 

VX  -  VX  -  Q1 
VTC  -  VY  -  Q2 
VZ  -  VZ  -  Q3 


RIGHT  TRAILER  CONTRIBUTION,  POINT  1  TO  POINT  2 


IFI  -  0 
XA  -  XRl 
YA  -  YRl 
ZA  -  ZRl 
XB  -  XR2 
YB  -  YR2 
ZB  -  ZR2 

CALL  VLINE 

VX  -  VX  4-  Q1 
VY  -  VY  +  Q2 
VZ  -  VZ  Q3 


RIGHT  TRAILER  IMAGE 


ZA 

«K 

-ZA 

ZB 

■■ 

-ZB 

CALL 

VLINE 

VX 

VX  - 

Q1 

VY 

VY  - 

Q2 

VZ 

- 

VZ  - 

Q3 

RIGHT  TRAILER  CONTRIBUTION,  POINT  2  TO  POINT  3 

IFI  - 

'  1 

XA  - 

XR2 

YA  « 

YR2 

ZA  - 

ZR2 

XB  - 

XR3 

YB  - 

YR3 

ZB  - 

ZR3 

CALL 

VLINE 

VX  - 

VX  +  Q1 

VY  - 

VY  +  Q2 

VZ  “ 

VZ  +  Q3 

RIGHT 

'  TRAILER  IMAGE 

ZA  “ 

-ZA 

ZB  - 

-ZB 
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190 

191 

C 

CALL 

VLINE 

192 

193 

C 

vx  - 

VX  -  Q1 

194 

VY  - 

VY  -  Q2 

195 

VZ  m 

VZ  -  Q3 

196 

197 

c 

c 

198  C  DIMENSIONALIZE 


199  C 


200 

C 

201 

202 

203 

GDR-  GAMMA/RADIUS 

VX  -  VX*GDR 

VY  -  VY*GDR 

204 

VZ  -=  VZ*GDR 

205 

C 

206 

RETURN 

207 

END 

208 

C 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 
3  4 
35 
3  6 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


C 

C 

SUBROUTINE  HWvTVEL  (H,UN,UMB,  RVZ,  RADIUS,  WSPD.DELH,  HMAX, 

•  HUMTVP,DXO,BDLAYM) 

C 

Q  *4tilr**A*lk*****ifc*4fA**1ft***4t*******ir*'*ik**4rik4t4r4r*A4r4r4i*4nfrA:lr* 

C  SUBROUTINE  HWJVEL  GENERATES  THE  VELOCITY  PROFILE  AND 

C  THE  FORCES  AND  OVERTURNING  MOMENTS  FOR  A  HUMAN  BEING 

C  AT  A  GIVEN  RADIUS 

c 

CHARACTER*!  ICONT(5) 

CHARACTER*!  TEMCHAR 
CHARACTER*!  KEY, KKEY, HUMTYP 
CHARACTSR*12  PTSFIL(4) 

CHARACTER* 50  COMM (2) 

C 

COMMON  /  CKEY/  KEY, KKEY 
COMMON  /CONSTS/  P I , RHO, FPSPKN, RHOD2 , DRC 
COMMON  /INPUTC/  ICONT, COMM, PTSFIL 
COMMON  /PERSON/  QP(12),DSET 

COMMON  /PROFIL/  RJ, ZBJ, ZHJ, ZMJ, UMJ, ZB, ZH, ZM, UM, CU, CY 
COMMON  /  UNITS/  lOUl, IOU4, IOU5, 10U6, IOU7, 10U8, IGRAPH 
C 

Q  AA*************************************************** 

C 

ICD  -  0 

CALL  HOMCLS(tCD) 

C 

C  - 

C  INPUT  FOR  DELH  AND  HMAX  COMES  FROM  INPUTV  STATUS  MENU 

C  - 

C 

DSET  »  DELH 
C 

IF(DSET.EQ.O.)  DELH  »  HMAX 
C 

DELH  =  DELH/RIJlIUS 
HMAX  =  HMAX /RADIUS 

NHPTS  =  IFIX((HMAX  -  BVZ) /DELH)  +  1 
C 

IF (DSET .EQ. 0 .) GOTO  50 

IF  ( IOU6  .NE  .  lOU!)  WRITE  (IOU6,  '  ("!")') 

C 

IF(I0U6.EQ.6)  write (IOU6, 93)  COMM ( 1 ) , COMM (2 ) 

93  FORMAK  lOX,  A50  ,  /  ,  lOX,  A50 ,  / /) 

C 

WRITE(I0U6, 100!) 

1001  FORMAT(  12X, '  SUMMARY  OF  OVERTURNING  FORCES  AND  MOMENTS',//, 

1  19X,  ' RADIUS' , 6X, ' TOTE' , 6X, ' TOTM'  , /, 

2  2  0X, '  (FT)  ' , 7X, '  (LB) ' , 5X, '  (FT -LB) ' , /) 

50  CONTINUE 


C  WRITE  OUT  GRAPHICS  FILE  IF  SWITCH  IS  SET  BY  USER 

C  - 

/*» 

IF { IGRAPH .EQ. 1) THEN 

C 

V..  »  —  ^  —  —  —  —  —  «  —  —  •  —  —  —  —  —  — 

C  OPEN  GRAPHICS  FILE  AND  WRITE  FILE  HEADER 

C  - 

n 
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64  OPEN (IOU8,FILE-PTSFIL( 3) , STATUS-' MEW» ,ERR-20 00) 

65  C 

66  WRITE (lOUe, 83)  COMM(l) ,COMM(2) 

67  83  FORMAT(  lOX, A50, /, lOX, A50, //) 

68  C 

69  WRITE (IOU8,  801 

70  80  FORMAT (  IX, ' TITLE-"SINGLE  ROTOR  DFRC  DATA"') 

71  C 

72  WRITE (IOU8, 81) 

73  81  FORMAT(  IX, ' VARIABLES  -  DFRC, TOTF, TOTM' ) 

74  C 

75  WRITE (IOU8,  88) 

76  88  FORMAT (  IX, ' ZONE  T  -  "GW  -  xxxxx  LB,  WAGL  -  xx  FT",', 

77  *  '  I-x,  F-POINT') 

78  C 

79  ENDIF 

80  C 

81  C  - 

82  C  BEGIN  LOOP  INCREMENTING  THE  RADIAL  POINTS  AT  WHICH 

83  C  THE  OVERTURNING  MOMENT  CALCULATIONS  ARE  MADE 

84  C  - 

85  C 

86  DO  565  K  -  1,NHPTS 

87  C 

88  C  - 

89  C  'PROPRM'  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 

90  C  OF  A  RADIAL  WALL  JET  (WITHOUT  INTERACTION  PLANE) 

91  c  - 

92  C 

93  CALL  PROPRM(H,UMB,RVZ) 

94  C 

95  ZETAM  -  ZM/ZB 

96  ZETAH  -  ZH/ZB 

97  C 

98  C  - 

99  C  CALCULATION  OF  THE  NON -DIMENSIONAL I ZED  MINIMUM  ALLOWED 

100  C  BOUNDARY  LAYER  THICKNESS  SO  THAT  THE  BOUNDARY  LAYER  CAN 

101  C  BE  ADJUSTED  IF  THE  ZM  POSITION  IS  PHYSICALLY  TOO  LOW 

102  C  (BDLAYM,  IN  FEET,  COMES  FROM  A  MENU  INPUT  PARAMETER) 

103  C  - 

104  C 

105  ZETAl  -  BDLAYM/ZB 

106  C 

107  C  - 

103  C  BOUNDARY  LAYER  REGION  EXPONENT 

109  C  'AN'  IS  ACTUALLY  '-  1. 0/7.0' 

110  C  - 

111  C 

112  AN  =  0.142857142 

113  C 

114  C  - 

115  C  SHEAR  LAYER  REGION  EXPONENT,  TO  MEET  EDGE  CONDITIONS 

116  C  (FROM  FIGURE  7,  USAAVLABS  TECHNICAL  REPORT  68-52,  JULY  1968) 

117  C  - 

118  C 

119  ALPW  -  ALOGd.O  -  1 . 0/SQRT  (2 . 0)  ) /ALOG  (  (ZH  -  ZM)/(ZB  -  ZM)  ) 

120  C 

121  VN  -  UN 

122  VMN  -  UM 

123  C 

124  C  - 

125  C  PRINT  DETAILED  REPORT  IF  DSET  =0.0  INSTEAD  OF  SIMPLE  REPORT 

126  C  - 
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127  C 

128  IF(DSET.NE.O.)GOTO  78 

129  C 

120  RRVZ  -  RVZ*RADIUS 

131  RVZOUT  -  RRVZ  +  DXO 

132  C 

133  IFdOOe.NE.IOUl)  WRITE  (IOU6, '("  1"  )' ) 

134  C 

135  1F(10U6.EQ.6)  WRITE ( IOU6, 93)  COMM ( 1) , COMM (2 ) 

136  C 

137  WRITE (10U6, 1000)  RVZOUT 

138  1000  FORMAT(  lOX, 'SINGLE  ROTOR  VELOCITY  PROFILE  AT  RADIUS  - 

139  1  F7.1,'  FT',//) 

140  WRITE (IOU6,  1005) 

141  1005  FORMAT (  2X, ' HEIGHT' , 5X, 'MEAN  VELOCITY' , 7X, ' PEAK  VELOCITY' , 6X, 

142  1  'MEAN  Q',4X,' PEAK  Q',/, 

143  2  3X, '  (FT) ' , 5X, '  (FPS) ' , 6X, '  (KN) ' ,5X, '  (FPS) ' , 6X, '  (KN) ' , 5X, 

144  3  '  (PSF) ' ,5X,'  (PSF)',/) 

145  78  CONTINUE 

146  C 

147  C  - 

148  C  SET  UP  ABILITY  TO  CALCULATE  AT  0.5  FT. 

149  C  INCREMENTS  UP  THE  VELOCITY  PROFILE 

150  C  - 

151  C 

152  DELZ  -  0.5/RADIUS 

153  NPTS  -  12 

154  C 

155  DO  500  I  =  1,NPTS 

156  C 

157  Z  -  DELZ*(I  -  1)  +  (0.25/RADIUS) 

158  ZETA  -  Z/ZB 

159  C 

160  IF ( ZETA . LT . ZETAM . OR . ZETA . LT . ZETAl ) THEN 

161  C 

162  C  - 

163  C  Z  IS  WITHIN  BOUNDARY  LAYER 

164  C 

165  C  NOTE  THAT  THE  BOUNDARY  LAYER  CALCULATIONS  NOW  USE 

166  C  THE  MINIMUM  THICKNESS  PARAMETER  AND  THE  PEAK  TO 

167  C  MEAN  VELOCITY  PARAMETER  IS  THE  MAXIMUM  VELOCITY 

168  C  HEIGHT  RATIO  (AT  ZM) .  ADDED  MAY  1992  FOR  V2 . 1 . 

169  C  - 

170  C 

171  VZM  -0.0 

172  C 

173  IF (ZETAM. GT. 0.0) THEN 

174  C 

175  VZM  -  ( ZETA/ ZETAM) ** AN 

176  C 

177  IF (ZETAl .GT. ZETAM) THEN 

178  C 

179  VZMl  -  (1.0  -  ((ZETAl  -  ZETAM)/ (1.0  -  ZETAM) ) **ALPW) **2 

180  VZM  -  VZMl* (ZETA/ ZETAl ) **AN 

181  C 

182  ENDIF 

183  C 

184  VMTOPK  -  1.04653  +  0.373894*RVZ  -  0 . 0422525*RVZ*RVZ 

185  C 

186  IF (VMTOPK.lt. 1.2)  VMTOPK  -  1.2 

187  C 

188  ENDIF 

189  C 
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190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 

215 

216 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GOTO  400 

END  IF 

Z  IS  WITHIN  SHEAR  LATER 

THE  PEAK  TO  MEAN  VELOCITY  RATIO  EQUATIONS  ARE 
SUBSTANTIALLY  IMPROVED  OVER  THOSE  USED  PRIOR  TO 
KAY  1992.  EQUATIONS  ARE  NOW  USED  FOR  BOTH  THE 
MAXIMUM  VELOCITY  HEIGHT  (ZM)  AND  THE  1/2  VELOCITY 
HEIGHT  (ZH) .  VALUES  BETWEEN  ARE  INTERPOLATED  AND 
VALUES  ABOVE  ZH  USE  THE  ZH  RATIO* ( ZETA/ZETAH) . 

THESE  2nd  ORDER  EQUATION  SUBSTANTIALLY  IMPROVED 
CORRELATION  WITH  MODEL  AND  FLIGHT  TEST  DATA 
DURING  THE  MAY  1992  EFFORT  FOR  V2 . 1 . 

VZM  -  0.0 

IF (Z .LE. ZB) THEN 

VZM  -  (1.0  -  ( (ZETA  -  ZETAM)/(1.0  -  ZETAM) ) **ALPW) * *2 
IF(ZETA.GE .ZETAH) THEN 

VMTOPK  -  (1.48086  +  0.569177*RVZ  -  0 . 0692514*RVZ*RVZ) 
1  * (ZETA/ZETAH) 

IF (VMTOPK.lt. 1.2)  VMTOPK  -  1.2 

ELSE 

VMPKMX  “  1.04653  +  0.373894*RVZ  -  0 .0422525*RVZ*RVZ 
VMPK12  -  1.48086  +  0.569177*RVZ  -  0 . 0692514*RVZ*RVZ 
FRAC  “  (ZETA  -  ZETAM) / (ZETAH  -  ZETAM) 

IF ( ZETAl . GT . ZETAM) THEN 
FRAC  -  (ZETA  -  ZETAl )/ (ZETAH  -  ZETAl) 

ENDIF 

VMTOPK  =  FRAC*VMPK12  +  (1.0  -  FRAC) * VMPKMX 
IF (VMTOPK.lt. 1.2)  VMTOPK  -  1.2 
ENDIF 

ENDIF 

400  CONTINUE 

VZN  -  V2M*VMN 

DIMENSIONAL  HEIGHT 

ZZ  -  Z*RADIUS 
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253  C  - 

254  C  MEAN  VELOCITIES 

255  C  - 

256  C 

257  VMF  -  VZN*VN 

258  VMK  -  VMF/FPSPKN 

259  C 

260  C  - 

261  C  PEAK  VELOCITIES 

262  C  - 

263  C 

264  VPF  -  VMF*VMTOPK 

265  VPK  -  VPF/FPSPKN 

266  C 

267  C  - 

268  C  THE  EFFECT  OF  WIND  IS  TO  ADD  (DOWNWIND  SIDE)  OR  SUBTRACT 

269  C  (UPWIND  SIDE)  ' XWK'  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 

270  C  THE  HORIZONTAL  PROFILE  VELOCITY  (EMPIRICAL,  CH-53E  BASED) 

271  C  - 

272  C 

273  XKW  -  (-0.5*H)  +  2.5 

274  C 

275  IF (XKW.LT. 1 .0)  XKW  »  1.0 
27  6  C 

277  WSPD2  -  WSPD*XKW 

270  VMK  -  VMK  +  WSPD2 

279  VMF  =  VMK*FPSPKN 

280  VPK  -  VPK  +  WSPD2 

281  VPF  ”  VPK*FPSPKN 

282  C 

283  C  - 

284  C  DYNAMIC  PRESSURE 

285  C  - 

286  C 

287  QM  -  RHOD2’*VMF**2 

288  QP(I)  -  RH0D2*VPF**2 

289  C 

290  IF(DSET.NE.0.)GOTO  77 

291  C 

292  WRITE  (lOr:. 1002)  ZZ, VMF,VMK,VPF, VPK,QM,QP (I) 

293  1002  FORMATC  F8.2,6F10.3) 

294  77  CONTINUE 

295  C 

296  500  CONTINUE 

297  C 

298  IF(DSET.NE.0.)GOTO  520 

299  C 

300  WRITE (lOUl, 73) 

301  73  FORMAT (  ) 

302  C 

303  WRITE (lOUl, '  (15X, A, S) ' ) 

304  1  '  TYPE  vRETURN>  TO  CONTINUE  ' 

305  C 

306  READdOUl,  '  (Al)  '  )  TEMCHAR 

307  C 

308  ICD  -  0 

309  CALL  HOMCLS(ICD) 

310  IF  (IOU6.NE. lOUl)  WRITE (IOU6, '("  1' ')' ) 

311  C 

312  IF (IOU6 .EQ. 6)  WRITE ( IOU6 , 93)  COMM ( 1 ) , COMM ( 2 ) 

313  C 

314  WRITE ( IOU6, 1007)  RVZOUT 

315  1007  FORMAT(  12X, 'SINGLE  ROTOR  FORCE  PROFILE  AT  RADIUS  = 


E-43 


SUBROUTINE  HWJVEL 


316  1  F7.1,'  FT',//) 

317  C 

318  WRITE (IOU6, 1008) 

319  1008  FORMAT(  2X, ' HEIGHT' , 6X, ' PEAK  Q' , 6X, ' FOVER' , 7X, ' OVERM' , 7X, 

320  1  'TOT  F' ,7X,'TOT  M', /, 

321  2  3X, '  (rT)',8X, '(PSF)',7X,'  (LB)',6X, '  (FT-I,B)',7X, 

322  3  ' (LB) ' ,6X,' (FT-tB)',/) 

323  C 

324  520  CONTINUE 

325  C 

326  C  - 

327  C  CALL  SUBROUTINE  TO  CALCULATE  THE 

328  C  FORCES  AND  MOMENTS  ON  A  HUMAN  BEING 

329  C  - 

330  C 

331  CALL  MOMENT  (NPT.'J.HUMTYP,TOTF,TOTM) 

332  C 

333  rF(DSET.EQ.0.)GOTO  545 

334  C 

335  HH  -  RVZ*RADIUS 

33  6  HHOUT  -  HH  f-  DXO 

337  C 

338  WRITE (IOU6, 1014)  HHOUT,TOTF, TOTH 

339  1014  FORMAT(  18X,  F8 .2,  T-FIO  .  3) 

340  C 

341  IFdGRAPH.EQ.DTHEN 

342  C 

343  WRITEdOUe,  90)  HMOUT,  TOTE  ,  TOTM 

344  90  FOPMAT(  IX, F7 . 2, IX, F7 . 2, IX, F8 . 2) 

345  C 

346  ENDIF 

347  C 

348  545  CONTINUE 

349  C 

350  RVZ  -  RV2  +  DELH 

351  C 

352  565  CONTINUE 

353  C 

354  C  - - 

355  C  CLOSE  Ai'l  OPEN  GRAPHICS  FILE 

356  C  - 

357  C 

358  IFdGRAPH.EQ.DTHEN 

359  C 

CLOSE  dOU8,  STATUS-' KEEP') 

361  C 

362  EUDIF 

363  C 

364  CALL  INKEY 

365  C 

366  GOTO  999 

367  C 

368  C  - 

335  C  THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 

370  C  OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 

371  C  - 

372  C 

373  2000  CONTINUE 

374  C 

375  CALL  HOMCLS(O) 

376  WRITEdOUl,  2001) 

377  2001  FORMAT (  ////,BX, 

373  1  '  ***  ERROR  ***  PLEASE  CHOOSE  A  NEW  OUTPUT  FILENAME', 


E-44 


SUBROUTINE  HWJVEL 


■xifi  2  ///  8X  '  <RETURiJ>  TO  CONTINUE  '  ,  $) 

380  REAd'(IOU1,MA1)M  TEHCHAR 

381  ><EV  -  'P' 

382  C 

383  999  CONTINUE 

384  C 

385  return 

386  Ein> 

387  C 

388 


SUBr.OUTlNE  INKEY 


1  C 

2  c 

3  SUBROUTINE  INKEY 

4  C 

5  Q  A***************************************************** 

6  C  SUBROUTINE  INKEY 

Q  Hr  *  *  *  *  Ik  *  #  *  *  «^*  *  *  **  *  *  *  *  *  *  **  *  W  A  **  A  *  *  *  *  llr  *  *  itr  *  ik  A  *  *  *  *  *  *  *  rtr  *  ** 

0  C 

9  PARAMETER (NUM  -  8) 

10  C 

11  CHARACTER*!  KEY,KKEY 

12  CHARACTER*!  OKLIST (NUM) 

13  C 

14  COMMON  !  CKEY/  KEY, KKEY 

15  COMMON  /  UNITS/  lOUl, IOU4, 1005, 10U6, IOU7, ICU8, IGRAPH 

16  C 

27  Q  */,****  Jr  ff«***A****«.#***-'*****«*****1lr*ilr**fl****#*****ik**** 

18  C 

19  DATA  OKIilST  /' C' , '  c' ,  '  P' , 'P' /  '  N' , '  n' ,  '  X' , '  x' / 

20  C 

21  C  - 

22  C  INQUIRE,  OBTAIN,  AND  CHECK  ?OR  VALID  MENU  OPTION 

23  C  - - - 

24  C 

25  10  CONTINUE 

26  C 

27  WRITE (lOUl,  20) 

28  20  FORMAT(  //,7X, '  TYPE  <C>ONTINUE,  NEXT  <P>OINT,  <N>EW  CASE.' 

29  1  ,'  E<X>IT  '->  ',$) 

30  READdOUl,  '  (Al)  '  )  KEY 

31  C 

32  IF (LEGAL (KEY, IOUl,OKLIST, NUM) .EQ.l) GOTO  10 

33  C 

34  c  - 

35  C  CORRECT  LOWER  CASE  LETTERS  TO  UPPER  CASE 

36  C  TO  USE  AS  VALID  FLAGS  IN  PARENT  SUBROUTINE 

37  C  - 

30  C 

39  IF(KEY.EQ.'c' )  KEY  -  'C' 

40  IF(KEY.EQ. 'p' )  KEY  -  'P' 

41  IF(KEY.EQ.'n' )  KEY-  'N' 

42  IF(KEY.EQ. 'X' )  KEY-  'X' 

43  C 

44  c  - 

45  C  CLEAR  SCREEN  Ai^D  HOME  CURSOR 

46  C  - 

47  C 

40  ICD  -  0 

49  CALL  HOMCLS(ICD) 

50  C 

51  RETURN 

52  END 

53  C 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 
10 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
4fl 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


C 

C 

SUBROUTINE  INPUT 
C 

C  SUBROUTINE  INPUT 
C 

C  THIS  SUBROUTINE  PRESENTS  THE  INPUT  STATUS  MENU 

C  AND  MANIPULUTES  THE  DATA  FOR  PROGRAM  USE 

c 

PARAMETER (NUMl  -  19) 

PARAMETER (NUM2  -  9) 

PARAMF.TER(NUM3  -  4) 

PARAMETER (NUM4  -  4) 

C 

CHARACTER* 1  CHDOE 
CHARACTER*!  GENTRY 
CHARACTER‘1  OKLSTl(NUMl) 

CHARACTER*!  OKLST2(NUM2) 

CHARACTER*!  0KLST3(NUM3) 

CHARACTER*!  OKLST4(NUM4) 

CHARACTER*!  ICONT(5) 

C 

CHARACTER* 12  PTSPIL(4) 

C 

CHARACTER*50  COMM (2) , LENTRY 
CHARACTER* 50  PROMPT 
C 

DIMENS  ION  CONT ( 9 ) , CONTV ( 7 ) , CONTX { 0 ) 

C 

COMMON  /INPUTC/  ICONT,  CO:4M,  PTSFIL 
COMMON  /INPUTD/  CONT, CONTV, COKTX, YBDEAY 
COMMON  !  UNITS/  lOUl, iOU4, IOU5, 1006. IOU7 , IOU8, IGRAPH 
C 

Q  ;4rA********AA***tt*4**<t*ilir****A**i*****»A*rt**'4r#’***|lt**'<l****<lr 

c 

C  - 

r  SET  DATA  TO  CHECK  FOR  ILLEGAL  DATA  INPUT 


DATA  OKLSTl 

1 

DATA  OKLST2 
DATA  OKLST3 
DATA  OKLST4 
C 

10  CONTINUE 
C 

C  - 

C  CLEAR  SCREEN  AND  HOME  CURSOR 

C  - 

C 

CALL  HOMCLS(O) 

CALL  LOCATE (2,1) 

C 

C  - - - - - 

C  WRITE  FIRST  ENGINEERING  DATA  MENU 

C  - 

C 

IROTOR  -  IFlX(CONT(I) ) 

C 

WRITEdOUl,  '  (24X,A/)  '  )  '  ROTWA5H  USER  INPUT  DATA  MENU' 


/' V' , ' V , ' H' , ' h' / 
/' Y'  ,  'y' , 'N' , ' n' / 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 
7P 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 
:  15 
116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


C 

C 

c 

c 

c 


c 

c 

c 

c 

c 


c 


c 

c 

c 

c 


c 

c 


c 

c 


WRITEdOUl,  '  (T7,A,T25,A,T51,A,T61,A/) ')  '  CODE','  PARAMETER', 
1  '  VALUE','  UNITS' 

WRITEdOUl,  '  (8X,A,T50,  IS,8X,A)  '  > 

1  '  A  NUMBER  OF  ROTORS  (1  OR  2)  ' , IROTOR, ' -ND-' 

WRITEdOUl,  '  (8X,A,T50,F8.1,  6X,A) ' ) 

1  '  B  HUB  TO  HUB  ROTOR  SEPARATION' , CONT (2) ,' FT' 

WRITEdOUl,  '  (8X,A,T50,F8.1,  6X,A) ' ) 

1  '  C  ROTOR  RADIUS 

WRITEdOUl,  '  (8X,A,T50,F8.1,6X,A)') 

1  '  D  GROSS  WEIGHT 

WRITEdOUl,  '  (8X,A,TS0,F8.1,6X,A)') 

1  '  E  FUSELAGE  DOWNLOAD  FACTOR 

WRITEdOUl,  '  (8X,A,T50,F8.1,6X,A>') 

1  '  F  ROTOR  HEIGHT  ABOVE  GROUND 

WRITEdOUl,  '  (8X,A,T50,F8.1,6X,A)') 

1  '  G  SHAFT  TILT  ANGLE  (<20  DEG) 

WRITEdOUl,  '  (8X,A,T50,F8.4,6X,A)') 

1  '  H  AIR  DENSITY  RATIO 

WRITEdOUl,  '  (8X,A,T50,F8.1,6X,A//)  ') 

1  '  I  AMBIENT  WIND  (-10  TO  10  KT) ' , CONT ( 9) , ' KT' 


' ,CONT(3), 'FT' 

' , CONT (4) , 'LB' 

' , CONT (5), 'PCT' 
' , CONT (6) , 'FT' 

' , CONT (7) , 'DEG' 
' , CONT (8), 'ND' 


PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 


20  WRITE (lOUl,  '  (8X,A, $) ' ) 

1'  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  — >  ' 
READ(IOUl, ' (Al) ' )  CHDOL 

IF  (LEGAL  (CHDOL,  lOUl ,  OKLSTl , NUMl )  .EQ.DGOTO  2  0 

IF  (CHDOL. EQ. '  ' ) GOTO  30 

IF  (CHDOL . EQ . ' A'  . OR . CHDOL , EQ . ' a ' ) THEN 

PROMPT  -  'ROTORS  -  ' 

CALL  IREAD (lOUl, PROMPT, IROTOR) 

IF ( IROTOR.lt. 1)  IROTOR  -  1 
IF (IROTOR. GT. 2)  IROTOR  -  2 
CONT(l)  -  FLOAT (IROTOR) 

GOTO  10 

ENDIF 

IF (CHDOL . EQ . ' B' . OR . CHDOL . EQ . ' b' ) THEN 

PROMPT  -  'ROTOR  SEPARATION  -  ' 

CALL  FREAD ( lOUl , PROMPT, CONT (2 ) ,1.0) 

GOTO  10 

ENDIF 

IF(CHDOL.EQ. 'C' . OR . CHDOL . EQ .' c' ) THEN 

PROMPT  “  'ROTOR  RADIUS  -  ' 

CALL  FREAD  (lOUl,  PRO.MPT,  CONT  (3)  ,  1.0) 

GOTO  10 

ENDIF 
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127 

128  C 

129 

130 

131 

132  C 

133 

134  C 

135 

136  C 

137 

138 

139 

140  C 

141 

142  C 

143 

144  C 

145 

146 

147 

140  C 

149 

150  C 

151 

152  C 

153 

154 

155  C 

156 

157 

158 

159  C 

160 

161  C 

162 

163  C 

164 

165 

166 

167  C 

168 

169  C 

170 

171  C 

172 

173 

174  C 

175 

176 

177 

178  C 

179 

180  C 
101 

182  C 
103 

184  C 
105  C 
■  106  C 

107  C 

108  C 

109 


ir (CHDOL.BQ.'D' . OR . CHDOL . EQ . ' d' ) THEN 

PROMPT  -  'GROSS  WEIG,  ■  ■  ' 

CALL  FREADdOUl^PROM'  :  ONT<4)  ,1.0) 
GOTO  10 

END  IF 

IF (CHDOL . EQ . ' E' . OR . CHDOL . EQ . ' ) THEN 

PROMPT  -  'DOWNLOAD  FACTOR  -  ' 

CALL  FREAD (lOUl, PROMPT, CONT (5), 1.0) 
GOTO  10 

END  IF 

IF (CHDOL. EQ. 'F' . OR. CHDOL. EQ.' f' )THEN 

PROMPT  -  'HEIGHT  ABOVE  GROUND  -  ' 

CALL  FREAD  ( lOUl ,  PROMPT,  CONT  (0,1.0) 
GOTO  10 

END  IF 

IF (CHDOL. EQ. 'G' .OR .CHDOL .EQ .' g' ) THEN 

PPOMPT  “  'SHAFT  TILT  -  ' 

CALL  FREAD (lOUl, PROMPT, CONT (7),  1.0) 

IF  (CONT  (7)  .LT.  0.0)  CO:  '’(7)  -  0.0 

IF(CONT(7) .GT.20.0)  CONT(7)  -  20.0 
GOTO  10 

ENDIF 

I F ( CHDOL . EQ .' H '. OR . CHDOL . EQ .' h ' ) THEN 

PROMPT  -  'DENSITY  RATIO  -  ' 

CALL  FRE,-D  ( ICUl ,  PROMPT,  CONT  ( B )  ,  1 . 0 ) 

got;  10 

ENDIF 

IF (CHDCL.EQ.' I'  . OR. CHDOL. EQ. '  i' )TBEN 

PROMPT  -  'WIND  VELOCITY  -  ' 

CALL  FREADdOUl, PROMPT, CONT(9)  ,1.0) 

IF(C0NT(9) .LT.-IO .0)  CONT(9)  -  -10.0 
IF (CONT { 9) .GT .  10.0)  CONT(9)  -  10.0 

GOTO  10 

ENDIF 

GOTO  10 

30  CONTINUE 


CLEAR  SCREEN  AND  HOME  CURSOR 


CALL  HOMCLS(O) 
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190  CALL  LOCATE <2, 1) 

191  C 

192  C  - 

193  C  WRITE  SECOND  ENGINEERING  DATA  MENU 

194  c  - 

195  C 

196  WRITE (lOUl, 40) 

197  40  FORMAT ('  ' , 18X, '  ROTWASH  PROGRAM  LOGIC/COMMENT  MENU',//, 

198  1  10X,'CODE  PARAMETER  VALUE'/) 

199  C 

200  WRITE (lOUl, 50)  ICONT (1) , ICONT (2) , COMM ( 1) , COMM (2 ) 

201  50  FORMATC  12X, 'A  ANALVSIS  TYPE,  <V>  OR  <H>' , 5X, A4, / ,  • 

202  1  12X, 'B  GRAPHICS  FILE,  <Y>  OR  <N>' , 5X, A4, ///, 

203  2  15X, 'USER  INPUT  COMMENTS  (FOR  "PRN"  AND  "PLT"  OUTPUT)', 

204  3  //,  17X,'< -  50  SPACES  - >' , 

205  4  /,12X, 'C',4X,A50,/, 

206  5  12X, 'D' , 4X,A50,//) 

207  C 

208  C  - 

209  C  PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 

210  C  - 

211  C 

212  60  WRITE (lOUl, ' (8X, A  $) ' ) 

213  1'  ENTER  CODE  FOR  DATA  INPUT  OR  <RETURN>  TO  CONTINUE  — >  ' 

214  C 

215  READdOUl,  '  (Al)  '  )  CHDOL 

216  C 

217  IF (LEGAL (CHDOL, IOU1,OKLST2,NUM2) .EQ.l) GOTO  60 

218  C 

219  IF (CHDOL. EQ.'  ' ) GOTO  70 

220  C 

221  C  - 

222  C  CHOOSE  VELOCITY  OR  HAZARD  ANALYSIS  OPTION 

223  C  - 

224  C 

225  IF (CHDOL. EQ. 'A' .OR.CHDOL.EQ. ' a' ) THEN 

226  C 

227  80  CONTINUE 

228  C 

229  WRITE (lOUl, ' (/, 35X,A, 1X,A2/) ' )  '  ANALYSIS  TYPE-  ',ICONT(l) 

230  WRITEdOUl,'  (40X,A,$)')  '  ENTER  NEW  VALUE  — >  ' 

231  READdOUl,  '  (Al)  '  )  CENTRY 

232  C 

233  IF (LEGAL (CENTRY, 10U1,0KLST3,NUM3) .EQ.l) GOTO  80 

234  C 

235  ICONT (1)  -  CENTRY 

236  IFdCONT(l)  .EQ.'v')  ICONTd)  -'V' 

237  IFdCONTd)  .EQ.'h')  ICONTd)  -  'H' 

238  GOTO  30 

239  C 

240  ENDIF 

241  C 

242  C  - 

243  C  CHOOSE  OUTPUT  TO  A  GRAPHICS  FILE  OR  NOT 

244  C  - 

245  C 

24  6  IF (CHDOL. EQ.'B' . OR. CHDOL. EQ. ' b' )  THEN 

247  C 

248  90  CONTINUE 

249  C 

250  WRITEdOUl,  '  (/,35X,A,1X,A2/)' )  '  GRAPHICS  FLAG  -  ',  ICONT  (2) 

251  WRITEdOUl,'  (40X,  A,  $)')  '  ENTER  NEW  VALUE  ==>  ' 

252  READdOUl,'  (Al)')  CENTRY 


E-50 


SUBROUTINE  INPUT 


253  C 

234  Zr(I.EOAL(CENTRy,  XOUl,OKLST4,NUM4)  .EQ.DGOTO  90 

255  C 

256  ICOMT{2)  -  GENTRY 

257  IF(ICONT(2) .EQ.'y' )  IC0NT(2)  -  'Y' 

258  IF(ICONT(2) .EQ.'n' )  1C0NT(2)  -  'N' 

259  IF(ICONT(2) .EQ.'Y' )  IGBAPH  -  1 

260  ZF(ZCONT<2) .EQ.'N' )  lGRAPH-0 

261  GOTO  30 

262  C 

263  ENDZF 

264  C 

265  C  - 

266  C  CHOOSE  COMMENT  STRINGS  FOR  OUTPUT  DATA 

267  C  - 

268  C 

269  1F(CHD0L.EQ.'C' .OR.CHDOL.EQ. ' c' ) THEN 

270  C 

271  MRITEdOUl,  '  (/,  5X,A,  A50/) '  )  '  COMMENT  STRING  -  '  ,  COMM  ( 1 ) 

272  MRZTEdOUl,'  (5X,A,  $)')  '  ENTER  STRING  — >  ' 

273  READdOUl,  '  (A50)  '  )  LENTRY 

274  C 

275  COMM(l)  -  LENTRY 

276  GOTO  30 

277  C 

278  ENDZF 

279  C 

280  IF(CHDOL.EQ.'D' .OR .CHDOL. EQ . '  d' ) THEN 

201  C 

202  WRITEdOUl,'  (/,5X,A,A50/)' )  '  COMMENT  STRING  -  ',COMM(2) 

283  WRZTEdOUl, '  (5X,A,  $) '  )  '  ENTER  STRING  — >  ' 

284  READ<I0U1, ' (A50) ' )  LENTRY 

285  C 

286  COMM(2)  -  LENTRY 

287  GOTO  30 

288  C 

289  ENDZF 

290  C 

291  GOTO  30 

292  C 

293  70  CONTINUE 

294  C 

795  RETURN 

296  END 

297  C 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  INPUTV (FLOW) 

«****«lr4r***ilr***A4r*******Aifc***A**4t*ir***ik*ilrAr4r**A*A*A*ir************ 

SUBROUTINE  INPUTV 

THIS  SUBROUTINE  PRESENTS  THE  INPUT  STATUS  MENU  AND  MANIPULATES 
DATA  FOR  FOR  THE  WALJET  AND  IPLANE  PROGRAM  OPTIONS 

4r*ilr*4r^*#*4rAAAr*****A*iikikr**,‘ir**AA4k4r*****Ai^***i^fk*****ilr****ilr**A**l^*** 

PARAMETER (NUM  -  11) 

CHARACTER*!  OKLIST(NUM) 

CHARACTER*!  CHDOL, FLOW 
CHARACTER*!  ICONT(5) 

CHARACTER*!2  TMPFIL 
CHARACTER*!2  PTSFIL(4) 

CHARACTER*50  COMM (2) 

CHARACTER*50  PROMPT 

DIMENSION  CONT(9) ,CONTV(7) ,CONTX(8> 

COMMON  /INPUTC/  ICONT, COMM, PTSFIL 

COMMON  /INPUTD/  CONT,  CONTV,  CONTX,  VTBDLAT 

COMMON  /  UNITS/  lOUl, IOU4, IOU5, IOU6, IOU7, IOU8, IGRAPH 

A************************************************************** 

DATA  OKLIST  /'  ' , ' A' , ' a' ,  ' B' , ' b' , ' C' , ' c' , ' D' , ' d'  , ' E' , ' e' / 


CLEAR  SCREEN  AND  HOME  CURSOR 


ICD  "  0 

20  CONTINUE 

CALL  HOMCLS(ICD) 

CALL  LOCATE (1, 1) 

WRITE (lOUl, 12) 

12  FORMAT (  22X,'  VELOCITY  PROFILE  STATUS  MENU',///, 

1  lOX, 'CODE  PARAMETER  VALUE', 

2  3X, 'UNITS' , /) 


PRINT  OUT  MENU  VARIABLES  AS  BASED  ON  THE  WALL  JET 
OPTION  OR  INTERACTION  PLANE  OPTION  SWITCH  SETTING 


IF(FLOW.EQ.'W'  )THEN 

WRITE (lOUl,  14)  (CONTV(I) , I-!, 3 )  ,  YBDLAY, PTSFIL ( ! ) 

ELSE 

WRITE  (lOUl,  14)  (CON .■'■(I)  ,1  =  1,3),  YBDLAY,  PTSFIL  (2) 

END  IF 

14  FORMAT(  12X, 'A  PROFILE  STATION  POSITION  ',5X,F7.2, 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
IOC 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 
119 
12C 
121 
122 

123 

124 

125 

126 


1  4X,'FT',/, 

2  12X,'B  VERTICAL  INCREMENT  ' , 4X, F7 . 2 , 4X, 

3  12X,'C  MAXIMUM  PROFILE  HEIGHT  ' , 4X, F7 . 2 , 4X, 

4  12X, ' D  MINIMUM  BOUNDARY  LAYER 

HEIGHT' , 1X,F7 .2, 4X, 'FT' ,11, 

5  12X,'E  DATA  OUTPUT  FILENAME  ',7X,A12,//) 

C 

C  - 

C  PROMPT  FOR,  OBTAIN,  AND  CHECK  FOR  LEGAL  INPUT  DATA 

C  - - - - 

c 

10  CONTINUE 
C 

WRITEdOUl,  '  {8X,A,  $)  '  ) 

1  '  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  — > 
C 

READ (lOUl, ' (Al) ' )  CHDOL 
C 

IF  (LEGAL  (CHDOL,  10U1,0KLIST,NUM)  .EQ.DGOTO  10 
C 

IF (CHDOL. EQ. '  ' ) GOTO  30 
C 

IF (CHDOL.EQ. ' A'  . OR .CHDOL . EQ . '  a' ) THEN 
C 

PROMPT  -  'PROFILE  STATION  POSITION  -  ' 

CALL  FREADdOUl,  PROMPT.  CONTV(l)  ,1.0) 

C 

IF(CONTV(l) .LT.O .0)  CONTV(l)  -  0.0 
GOTO  20 
C 

ENDIF 

C 

IF (CHDOL . EQ . ' B' . OR . CHDOL . EQ . ' b' ) THEN 
C 

PROMPT  -  'VERTICAL  INCREMENT  -  ' 

CALL  FREAD  dOUl,  PROMPT,  CONTV  (2)  ,1.0) 

C 

IF(CONTV(2)  .LT.O  .0)  CONTV(2)  -  0.0 
GOTO  20 
C 

ENDIF 

C 

IF (CHDOL . EQ . ' C' . OP . CHDOL . EQ . ' c' ) THEN 
C 

PROMPT  -  'MAXIMUM  PROFILE  HEIGHT  -  ' 

CALL  FREAD  dOUl,  PROMPT,  CONTV  (3)  ,1.0) 

C 

IF  (CONTV  (3)  .LT.O  .0)  CONTV  (3)  »().') 

GOTO  20 
C 

ENDIF 

C 

IF (CHDOL.EQ. ' D' . OR . CHDOL . EQ .' d' ) THEN 
C 

PROMPT  -  'MINIMUM  BOUNDARY  LAYER  HEIGHT  =  ' 

CALL  FREADdOUl,  PROMPT,  YBD  LAY,  1.0) 

IF (YBDLAY.lt. 0.0)  YBDLAY  =0.0 
GOTO  20 
C 

ENDIF 

C 

C  - 


'FT',/, 
'FT' ,//, 
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127 

128 
129 

C 

CHOOSE  GRAPHICS  rZEENAME 

C 

130 

ir (CHDOL . EQ . ' E' . OR . CHDOL . EQ . ' a' ) THEN 

131 

C 

132 

IF ( FLOW. EQ.'W') then 

133 

WRITE dOUl,' (/,25X,A,1X,A12/)') 

134 

1  '  FILENAME  -  ',PTSFIL(1) 

135 

ELSE 

136 

WRITEdOUl,'  (/,25X,A,  1X,A12/)') 

137 

1  '  FILENAME  -  ',PTSFIL(2) 

138 

END  IF 

139 

C 

140 

WRITEdOUl,  '  (20X,A,  S)  ' ) 

141 

1  '  ENTER  NEW  FILENAME  (xxxxxxxx.xxx) 

142 

C 

143 

READdOUl,  '  (A12)  '  )  TM.7FIL 

144 

C 

145 

1F(FL0W.EQ.'W' )  PTSFILd)  -TMPFIL 

146 

IF(FLOW.EQ.'I')  PTSFIL(2)  -TMPFIL 

147 

GOTO  20 

148 

C 

149 

ENDIF 

150 

c 

151 

GOTO  20 

152 

c 

153 

30  CONTINUE 

154 

c 

155 

ICD  -  0 

156 

CALL  HOMCLSdCD) 

157 

c 

158 

RETURN 

159 

END 

160 

c 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


C 

c 

SUBROUTINE  INPUTX 
C 

0  ***1kA*4t****^***4rAA*******ik****4t*ik*#Atk**1l*ibikA»*A:ikA*ir4r* 

C  SUBROUTINE  INPUTX  PRESENTS  THE  INPUT  STATUS  MENU  AND 

C  MANIPULATES  DATA  FOR  PROGRAM  USE  WITH  THE  GROUND  AND 

C  DISC  VORTEX  OPTIONS 

Q  ************iii**************************************** 

C 

PARAMETER(NUM  -  15) 

C 

CHARACTER*!  OKLIST(NUM) 

CHARACTER*!  CHDOL 
CHARACTER*!  ICONT(5) 

CHARACTER*12  PTSFIL(4) 

CHARACTER*50  COMM (2) 

CHARACTEB*50  PROMPT 
C 

DIMENSION  CONT (9) , CONTV ( 7 ) , CONTX ( 8 ) 

C 

COM!40N  /INPUTC/  ICONT,  COMM,  PTSFIL 
COMMON  /INPUTD/  CONT, CONTV, CONTX, TBDLAY 
COMMON  /  UNITS/  lOUl , IOU4 , IOU5, IOU6 . IOU7 , IOU8 , IGRAPH 
C 

c 

DATA  OKLIST  /'  '  ,  '  A' ,  '  a' ,  '  B' , '  b' , '  C' ,  '  c' ,  '  D' ,  '  d' , 

1  'E','e','F','f','G','g'/ 

C 

c  - - - - - 

C  CLEAR  SCREEN  AND  HOME  CURSOR 

C  - 

C 

ICD  =  0 
C 

20  CONTINUE 


C 


c 


c 

c 

c 

c 

c 

c 


CALL  HOMCLS(ICD) 
CALL  LOCATE (1,1) 


WRITE (lOUl, 10) 

10  rORMAT(  18X,  '  GROUND/DISK 
1  15X, '  (FOR  SINGLE 

1  8X,'CODE 

2  4X,  "UNITS' , /) 


VORTEX  INPUT  DATA  MENU',/, 

MAIN  ROTOR  HELICOPTERS  ONLY)',///, 
PARAMETER  VALUE' , 


WRITE (lOUl,  11)  (CONTX  (I) ,1  =  1,7) 

11  FORMAT (  10X,'A  ROTOR  TIP  SPEED  ',5X,F8.2, 

1  4X,'FPS',/, 


2 

lOX,  'B 

NUMBER  OF  ROTOR  BLADES 

,  5X,  F6 . 2, 

4X, 

'  -ND- 

3 

lOX, ' C 

TRANSLATIONAL  SPEED 

,  5X,  F8 . 2, 

4X, ' KTS' 

4 

ICX, 'D 

XT  POSITION  ' 

,  5X,F8.2, 

4X, ' FT' , 

5 

lOX, 'E 

YT  POSITION 

,  5X, F8 . 2, 

4X, ' FT' , 

6 

ICX, ' F 

ZT  CALCULATION  INCREMENT  ' 

,  5X,F8.2, 

4X. 

'FT'  , 

7 

lOX,  'G 

MAXIMUM  CALCULATION  HEIGHT' 

, 5X, F8 .2, 

4X, 

'  FT'  , 

INQUIRE,  OBTAIN,  AND  CHECK  FOR  VALID  MENU  OPTION 


40  CONTINUE 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 
i09 

no 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


C 

C 

c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 


c 

c 

c 


c 

c 

c 


c 

c 

c 


c 


c 


WRITE (lOUl, '  (8X,A,$)  '  ) 

1  '  ENTER  DATA  ENTRY  CODE  OR  <RETURN>  TO  CONTINUE  — >  ' 
READdOUl, '  (Al)  '  )  CHDOL 

IF  {LEGAL (CHDOL,  IOUl,OKLISr,NUM)  .EQ.DGOTO  40 

IF (CHDOL. EQ.'  ' ) GOTO  30 

IF (CHDOL . EQ . ' A' . OR . CHDOL . EQ . ' a ' ) THEN 

PROMPT  -  'ROTOR  TIP  SPEED  -  ' 

CALL  FREADdOUl, PROMPT, CONTXd)  ,1.0) 

IF(CONTX(l)  .LT.0.0)  CONTXd)  -  0.0 
GOTO  20 

ENOIF 

IF (CHDOL. EQ.'B' . OR . CHDOL . EQ b' ) THEN 

PROMPT  -  'NUMBER  OF  ROTOR  BLADES  -  ' 

CALL  FREAD ( lOUl , PROMPT, CONTX (2 )  ,  1 . 0 ) 

1F(C0NTX(2) .LT.2.0)  CONTX(2)  -  2.0 
GOTO  20 

ENDIF 

IF (CHDOL. EO. 'C' . OR. CHDOL. EQ. 'c' ) THEN 

PROMPT  -  'TRANSLATIONAL  SPEED  -  ' 

CALL  FREADdOUl,  PROMPT,  CONTX  (3)  ,1.0) 

IF  (CONTX  (3)  .LT.  0.0)  CONTXd)  -  0.0 
GOTO  20 

ENDIF 

IF(CHDOL.EQ. 'D' . OR . CHDOL . EQ .' d' ) THEN 

PROMPT  -  'XT  POSITION  -  ' 

CALL  FREAD ( lOUl , PROMPT, CONTX ( 4 )  ,  1 . 0 ) 

GOTO  20 

ENDIF 

IF (CHDOL. EQ. 'E' . OR . CHDOL . EQ .' e' )THEN 

PROMPT  -  'YT  POSITION  -  ' 

CALL  FREAD (lOUl, PROMPT, CONTX (5) , 1 .0) 

GOTO  20 

ENDIF 

IF (CHDOL. EQ. ' F' . OB . CHDOL . EQ . ' f ' ) THEN 

PROMPT  -  'ZT  CALCULATION  INCREMENT  -  ' 

CALL  FREADdOUl,  PROMPT,  CONTX  (6),  1.0) 

IF(CONTX(6) .LT.0.0)  CONTX(6)  -  0.0 
GOTO  20 
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127 

ENDIF 

128 

C 

129 

IF  (CHDOIi.EQ. 'G'  ,OR .CHDOL  . EQ  . ' g'  >  THEN 

130 

C 

131 

PROMPT  -  'MAXIMUM  CALCULATION  HEIGHT 

132 

CALL  FREADdOUl,  PROMPT, CONTX (7)  ,1.0) 

133 

C 

134 

IF{CONTX(7) .LT.0.0)  CONTX(7)  -  0.0 

135 

GOTO  20 

136 

C 

137 

ENDIF 

138 

C 

139 

GOTO  20 

HO 

C 

141 

30  CONTINUE 

142 

C 

143 

ICD  -  0 

144 

CALL  HOMCLS(ICD) 

145 

C 

1<16 

RETURN 

147 

END 

148 

149 
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1  C 

2  C 

3 

4  C 
3  C 
6  C 
1  C 

8  C 

9 

10  c 

11  c 

12  C 
1.3  C 
14  r: 
li)  C 

16  C 

17 

18 

19  C 

20 
21 

22  C 


SUBROUTINE  lOFNSH 

SUBROUTINE  lOFNSH  CLOSES  FILES  OPENED  FOR  DISK  I/O 

****4rAii^****ifr*****ilrA*****A**i**llfl4r****llnik*A**flrAik4r***illri^ 

COMMON  /  UNITS/  lOUl, I0U4, IOU5, IOU6, IOU7, IOU8, IGRAPH 


KEEP  STANDARD  I/O  FILES  IF  ON  DISK,  ELSE  DELETE 


IB  (IOU5 . EQ. 8)  CLOSE (IOU5, STATUS-' KEEP' ) 
IFdOUe.EQ.  6)  CLOSE  (IOU6,  STATUS-' KEEP'  ) 

RETURN 

Er'D 
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1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
IT 
18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 
46 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
GO 
61 
62 
63 


C 

c 

SUBROUTINE  lOINIT 
C 

Q  **#*ik***W*A***<tA***'ir4r***********4rA-****iArAtkilr*’**4rllt*i*r*'lri^1^ 

C  SUBROUTINE  lOINIT  DISPLAYS  THE  OPENING  BANNER 
C  AND  OPENS  THE  FILES  FOR  DISK  I/O,  FILENAMES 

C  ARE  PROMPTED  FROM  THE  TERMINAL 

Q  1tr4rA****1kA***lk^***4t*A****4^AA*iDt1lr*ik**i^A*A-/b*4r'/^*i*A**4rilt*«t* 

c 

CHARACTER*3  IPFILE, OPFILE 
CHARACTER*!  TEMCHAR 
C 

COMMON  /  UNITS/  lOUl, IOU4, IOU5, IOU6, IOU7, IOU8 , IGRAPH 
C 

c 

C  - 

C  ASSIGN  DEFAULT  VALUES  TO  I/O  UNIT  POINTERS 

C  - 

C 

lOUl  -  0 

10U4  -  0 

10U5  -  0 

10U6  =  6 

IOU7  -  0 

IOU8  -  8 

IGRAPH  -  0 

C 

C  . - . — . . 

C  HOME  CURSOR  AND  CLEAR  THE  DISPLAY 

C  - - - - - 

C 

ICD  >«  0 

CALL  HOMCLS(ICD) 

C 

C  - 

C  DISPLAY  BANNER 

C  - 

C 

CALL  LOCATE (5,1) 

C 

C  - 

C  ORIGINAL  ROTHAZ  PROGRAM  WAS  VERSION  1.0 

C  FIRST  VAX  MENU  VERSION  OF  ROTHAZ  WAS  VERSION  1.1 

C  PROGRAM  ROTWASH  REPLACES  ROTHAZ  AT  VERSION  2.0 

C  - 

C 

WHITE (lOUl, 10) 

10  FORMAT (  27X, 'ROTWASH  PROGRAM',// 

1  /,  17X, ' ROTORCRAFT  DOWNWASK  HAZARD  ANALYSIS',/ 

2  /,18X,'EMA  /  SYSTEMS  CONTROL  TECHNOLOGY', 

3  /,15X,'***  PROGRAM  VERSION  2.1,  APRIL  1993  ***' 

4  ///////) 

C 

WRITE (lOUl, ' (A, $) ' ) 

1  '  PRESS  <RETURN> 

READdOUl,  '  (Al)  '  )  TEMCHAR 
C 

C  - 

C  HOME  CURSOR  AND  CLEAR  THE  DISPLAY 

C  - 

C 
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65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 
63 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 
109 

no 

111 

112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


C 

C 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 


c 

c 

c 


ICD  -  0 

CALL  HOMCi.S(lCD) 
CALL  LOCATE (3,1) 


PROMPT  FOR  I/O  FILES  ANP  READ  USER  RESPONSE 


WRITE (1001,12) 


12  FORMAT ( 

14X,  ' 

I/O  CAN 

BE  DIRECTED  TO  FILES  OR  DEVICES 

1 

19X,  ' 

VALID  DEVICES  ARE  AS  FOLLOWS:',//, 

2 

22X,  ' 

<CON> 

— >  CONSOLE'  ,  / , 

2 

22X,  ' 

<PRN> 

— >  PRINTER',/, 

4 

22X,  ' 

<?LT> 

— >  GRAPHICS  FILE',///) 

16  CONTINUE 

WRITE (lOUl, 13) 

13  FORMAT (  16X, '  ENTER  INPUT  FiLE/DEV  NAME  — >  ',5) 
11  FORMAT (A) 

READdOUl,  11)  IPFILE 


CHECK  FOR  DATA  ENTRY  ERROR 


IFdPFILE.NE. 'CON'  .  AND  .  IPFILE .  NE .' COn'  )  THEN 
WRITEdOUl,  15) 

15  FORMAT (/, 16X, '  **  INPUT  ERROR,  PLEASE  REENTER  **',//) 

GOTO  16 
END  IF 

17  CONTINUE 

WRITEdOUl,  14) 

14  FORMAT (  16X,'  ENTER  OUTPUT  FILE/DEV  NAME  — >  ' , S) 
READdOUl,  11)  OPFILE 


CHECK  FOR  DATA  ENTRY  ERROR 


IF (OPFILE. NE.' CON' .AND .OPFILE . NE .' con' .AND. 

1  OPFILE.NE.'PRN' . AND . OPFILE . NE .' prn' .AND. 

2  OPFILE. NE. 'PLT' . AND . OPFILE . NE .' pit ') THEN 

WRITE  dOUl,  15) 

GOTO  17 
END  IF 


REDIRECT  INPUT/OUTPUT  FILES  IF  REQUESTED 


IFdPFILE. EQ. 'CON'  .OR.  IPFILE. EQ. 'con' )  IOU5  -  lOUl 
IF(OPFILE. EQ. 'CON' .OR. OPFILE. EQ. 'con' ) IOU6  -  lOUl 
IF(OPFILE. EQ. 'PRN' .OR. OPFILE. EQ. 'prn' ) IOU6  -  6 

IF(OPFILE. EQ. 'PLT' .OR. OPFILE. EQ. 'pit') IGRAPH  -  1 
IF(OPFILE. EQ. 'PLT' .OR. OPFILE. EQ. 'pit') IOU6  -  lOUl 
IF{OPFILE. EQ. 'PLT'  .OR. OPFILE. EQ. 'pit') IOU8  -  8 


OPEN  STANDARD  I/C  FILES 


E-60 


12*7  C 
12ij  C 

129  C 

130  C 

131  C 

132 

133 

134  C 

135  C 

136  C 

137  C 

138  C 

139 

140 

1^1  C 

142 

143 

144  C 


IOU5  -  STO.  INPUT 
IOU6  -  STD.  OUTPUT 


[F (IOU5 .BQ . 5)  OPEN (10U5, F IBB- IPFItE, 


STATUS- 


HOME  CURSOR  MTO  CLEAR  THE  OISPLAY 


ICD  -  0 

CALL  HOMCLSdCD) 

RETURN 

END 


SUBROUTINE  lOT.NIT 


■  'OLD' ) 
^'NBW' ) 
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1  C 

2  C 

3  SUBROUTINE  IPVEL (H,  UN, RADIUS, UMB, XIP, VSEP, WSPD, DELZ, 

4  1  ZMAX,DXO,BDLAYM) 

5  C 

g  ^  l^**lkA**«lrib*Ar)k***iAt***4rip*******ilr**^**#*******«llk******4r*A***** 

1  C  SUBROUTINE  IPVEL  GENERATES  THE  VELOCITY  PROFILE  V(X, Z)  AT 

8  C  XVZ  ALONG  THE  INTERACTION  PLANE  FOR  THE  TWO  ROTOR  CASE 

^  Q  *  *  *  Ik  ilr  *  A  *  *  llr  *  1/^  Af  ifr  *  *  *  A  *  *  *  Ar  *  *  Hr  *  4r  *  A  *  *  *  *  *  #  Ik  llr  A  4r  *  *  *  *  *  A  *  1^  *  Itr  itr  *  Hr  *  * 

10  c 

11  CHARACTER* 1  TEMCHAR 

12  CHARACTER*!  rCONT(5) 

13  CHARACTER*!  KEY, KKEY 

14  CHARACTER*12  PTSFIL(4) 

15  CHARACTER* 50  COMM (2) 

16  C 

17  DIMENSION  ZZ ( 60 ) , VHMF ( 60 ) , VHMK ( 60) , VHPF ( 60) , VHPK ( 60 ) , 

18  1  QHM(60)  ,QHP(60)  ,VVMF(60)  ,VVMK(60)  ,  VVTPF(60)  , 

19  2  VVPK(60) ,QVM(60),QVP(60) 

20  C 

21  COMMON  /  CKEY/  KEY,KKEY 

22  COMMON  /CONSTS/  PI, RHO, FPSPKN, RHOD2 , URC 

23  COMMON  /INPUTC/  ICONT, COMM, PTSFIL 

24  COMMON  /PROFIL/  RJ, ZBJ, ZHJ, ZMJ, UMJ, ZB, ZH, ZM, UM, CU, CY 

25  COMMON  /  UNITS/  lOUl, ICU4, IOU5, IOU6, IOU7, IOU8, IGRAPH 

26  C 

27  c  ********************************************************** 

23  C 

29  c  - 

30  C  TF  IS  THE  INTERACTION  PLANE  AMPLIFICATION  FACTOR 

31  C  (ORIGINALLY  DEVELOPED  BY  M.  GEORGE  IN  USAAVLABS  TR  68-52) 

32  C 

33  C  ORIGINAL  EQUATION  FOR  TF  FACTOR  WAS: 

34  C 

35  C  TF  -  1.55  -  0.55*EXP(-1.35*XIP) 

36  C 

37  C  REPLACED  WITH  MODIFIED  EXPRESSION  (SEE  BELOW)  DURING 

38  C  CORRELATION  EFFORT  CF  MAY  1992  FOR  VERSION  2.1 

39  c  - 

40  C 

41  TF  "  1.65  -  0.65*EXP (-0.5*XIP) 

42  C 

43  C  - 

44  C  OBTAIN  PARAMETERS  AT  BASE  RADIUS  FOR  THE  'BOUNDARY  LAYER' 

45  C  - 

46  C 

47  RIPO  -  SQRT(XIP**2  +  YSEP**2) 

4  8  C 

49  c  - 

50  C  'PROPRW'  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 

51  C  OF  A  RADIAL  WALL  JET  (WITHOUT  INTERACTION  PLANE) 

32  c  - 

55  C 

54  CALL  PROPRM(H,UMB,RIP0) 

55  C 

56  ZIPS  -  ZB 

57  ZIPK  -  2M 

58  ZIPH  =  ZH 

59  C 

60  RIPM  -  SQR"'(XIP**2  +  (YSEP  +  ZIPM)**2) 

61  C 

62  CALL  PROFRM(H,UMB, RIPM) 

63  C 
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64  UMM  -  UM 

65  C 

66  C  - 

67  C  INCREMENT  AND  MAXIMUM  HEIGHT  ARE  DELZ  ANO  ZMAX 

68  C  - - - 

69  C 

VO  NETS  -  IFIX(ZMAX/DELZ)  +  2 

71  C 

72  IF (NETS .GT .60)  NPTS  -  60 

73  C 

74  C  - 

75  C  DIMENSIONALIZE  VELOCITY  PROFILE  PARAMETERS 

76  C  - 

77  C 

78  XXIt-  »  RADIUS*XIP 

79  ZZB  -  ZIPB*RADIUS 

80  ZZH  “  ZIPH*RADIUS 

81  ZZM  -  ZIPM*RADIUS 

82  C 

83  c  - 

84  C  OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILE  HEADER 

85  c  - 

86  C 

87  ICD  -  0 

3  9  CALL  HOMCT  SdCD) 

89  C 

90  IF  ilOU6  .NE.  lOUl)  WRITE  (10U6, '("  1"  )' ) 

91  C 

92  IF (IOU6 .EQ. 6)  WRITE ( 10U6, 93)  COMM ( 1) , COMM (2) 

33  93  FORMAK  lOX,  A50, /,  10X,A50, //) 

94  C 

95  XIPOUT  »  XXIF  +  DXQ 

96  C 

97  WRITE (IOU6, 1000)  XIPOUT 

98  1000  FORhiATI  2X, 'TW7N  ROTOR  INTERACTION  PLANE  VELOCITY  PROFILE', 

99  1  'AT  STATION  -  ',F7.1,'  FT',//) 

100  c 

101  WRITE  (ICiU6,  1002) 

102  1002  FORMAT (  2X, ' HEIGHT' , 8X, ' MEAN  VELOCITY' , 7X, ' PEAK  VELOCITY' , 6X, 

103  1  ' MEAN  Q' , 4X, ' PEAK  Q' , /, 

1C4  2  3X, '  (FT) ' , 8X, '  (FPS) ' , 6X, '  (KN) ' , 5X, '  (FPS) ' , 6X, '  (KN) ' , 5X, 

105  3  '  (PSF) ' , 5X, '  (PSF) ' , /) 

106  C 

107  LINES  =  0 

108  C 

109  C  - 

110  C  'AN'  IS  ACTUALLY  ' 1. 0/7.0' 

11]  C  - 

112  C 

113  AN  =  0.142857142 

114  C 

115  C  - 

116  C  CALCULATE  THE  VELOCITY  PROFILE  POINTS  FOR  OUTPUT 

117  C  'NPTS'  IS  THE  NUMBER  OF  VERTICAL  STATION  POINTS 

118  C  - 

119  C 

120  DO  500  1  “  1,NPTS 

121  C 

122  LINES  =  LINES  +  2 

123  ZIP  =  DELZ*FLOAT(I  -  1) 

124  C 

125  C  - 

126  C  GET  MAX  WALL  JET  VELOCITY  AT  EFFECTIVE  RADIUS 
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127  C  - 

128  C 

129  RIP  -  SQRT<XIP**2  +  (YSEP  +  ZIP)**2) 

130  C 

131  CALL  PR0PRM(H,UMB,R1P) 

:32  C 

133  VN  -  UN 

134  VZ  -  UM 

135  C 

136  C  - 

137  C  INTERACTION  PLANE  'BOUNDARY  LAYER' 

138  C 

139  C  CODE  MODIFIED  IN  MAY  1992  FOR  USER  SPECIFIED 

140  C  MINIMUM  BOUNDARY  LAYER  THICKNESS  (BDLAYM) 

141  c  - 

142  C 

143  ZIPl  =  BDLAYM 

144  C 

145  IF (ZIP .LT. ZIPM.OR.ZIP .LT. ZIPl) THEN 

146  C 

147  IF<ZIP1.LT.ZIPM)THEN 

148  C 

149  VZ  -  UMM* (ZIP/ZIPM)  **AN 

150  C 

151  ELSE 

152  C 

153  VZ  -  UMM* (ZIP/ZIPI)  **AN 

154  C 

155  ENDIF 

156  C 

157  ENDIF 

158  C 

159  c  - 

160  C  DEVELOPED  INTERACTION  PLANE  JET 

161  C  CONTAINS  BOTH  HORIZONTAL  AND  VERTICAL  VELOCITY  COMPONENTS 

162  C  - 

163  C 

164  VH  -  TF*VZ*XIP/BIP 

165  VV  =  TF*VZ*(YSEP  +  ZIP) /RIP 

166  C 

167  ZZ(I)  =  ZIP*RADIUS 

168  C 

169  C  - 

170  C  MEAN  VELOCITIES  (BOTH  FT/SEC  AND  KNOTS) 

171  c  - 

172  C 

173  VHMF(I)  =  VH*UN 

174  WMFd)  •  VV*UN 

175  VHMK(I)  -  VHMF (I) /FPSPKN 

176  WMKd)  WMFd) /FPSPKN 

177  C 

178  C  - 

179  C  PEAK  VELOCITIES  (BOTH  FT/SEC  AND  KNOTS) 

180  C 

181  C  EQUATION  FOR  VMFD3I  UPDATED  FROM  Ist  TO 

182  C  2nd  ORDER  POLYNOMIAL  FOR  VERSION  2.1 

183  C  - 

18  4  C 

185  VMFD3I  =  0.712887  +  0.304369*XIP  -  0 . 018496*XIP*XIP 

106  C 

187  IF (VMFD3I.lt. 1.2)  VMFD3I  =  1.2 

188  C 

189  VHPF(I)  =  VMFD3I*VHMF (I) 
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190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 
20'' 
208 

209 

210 
211 
212 

213 

214 
2\5 
2  i  b 

217 

218 

219 

220 
221 
222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 

240 

241 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 


C 

C 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 


c 


wpr(i)  -  VMrD3i*wMr (I) 
VHPK(I)  -  VHPP(I) /FPSPKN 
WPK(l)  -  WPFd) /FPSPKN 

IF(VHMF(1)  .EQ.OJGOTO  55 


THE  EFFECT  OF  WIND  IS  TO  ADD  {DOWNWIND  SIDE)  OR  SUBTRACT 
(UPWIND  SIDE)  'XWK'  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 
THE  HORIZONTAL  PROFILE  VELOCITY  (EMPIRICAL,  CH-53E  BASED) 


XKW  -  (-0.5*H)  +  2.5 

IFOCKW.LT.l.O)  XKW  -  1.0 

WSPD2  -  WSPD*XKW 

VHMK(l)  -  VHMK(I)  +  WSPD2 

VHMF(I)  -  VHMK (I ) ‘FPSPKN 

VHPK(I)  -  VHPK(I)  +  MSPD2 

VHPF(I)  -  VHPK(I) ‘FPSPKN 

CONTI  N-JE 


DYNAMIC  PRESSURES 


QHM(I)  -  RHOD2*VHMF(I)“2 
QVM(I)  -  RH0D2*WMF(1)“2 
QHP(I)  -  RH0D2*VHPF(I)  “'2 
QVP(I>  -  RH0D2‘WPF(I)“2 

IF ( IOU6 . EQ . lOUl) THEN 

IF  (LINES.lt. 12) GOTO  450 
LINES  -  2 
CALL  INKEY 

IF (KEY.NE.'C' )GOTO  999 
MRITEdOUe,  1002) 

ENDIF 


C 


450 

C 

C 

C 

c 

c 


1003 

C 

C 

C 

c 

c 


1004 

C 


CONTINUE 


REPORT  HORIZONTAL  COMPONENTS  OF  VELOCITY  PROFILE 


WRITE  (10U6,  1003)  2Z  (I)  ,  VHMF(I)  ,VHMK(I)  ,  VHPFd)  , 
VHPK(I)  ,QHM(1)  ,QHP  (I) 

FORMAT  (  F8.2,2X,'H' ,6F10.3) 


REPORT  VERTICAL  COMPONENTS  OF  VELOCITY  PROFILE 


WRITE (IOU6, 1004) 
FORMAT  (  10X,'V' 


WMF  (I)  ,  WMKd)  ,  WPFd)  ,  WPK  (1)  , 
QVMd)  ,QVP(I) 

,6F10.3) 


500  CONTINUE 


C 


E-65 


SUBROUTINE  IPVEL 


253  C  - - - 

254  C  WRITE  OUT  GRAPHICS  FILES  IF  SWITCH  IS  SET  BY  USER 

255  C  - 

256  C 

257  IF(IGRAPH.EQ.1)THEN 

258  C 

259  C  - 

260  C  OPEN  GRAPHICS  FILE 

261  C  . . . . - 

262  C 

2  63  OPEN (IOU8, FILE=PTSF1L(2)  ,  STATUS-' NEW'  ,  ERR-2000) 

264  C 

265  WRITE (IOU8, 89)  COMM (1) , COMM (2) 

266  89  FORMAT(  lOX, A50, /, lOX,  A50, //) 

267  C 

268  WRITE (IOU8. 001  XIPOUT 

269  8C  FORMAT (  IX, ' TITLE-"VELOCITY  PROFILE,  DAIP  =',F5.1,'  FT,' 

270  *  'GW  -  xxxxx  LB,  WAGE  -  xx . x  FT"') 

271  C 

272  C  - - - 

273  C  PRINT  OUT  MEAN  VELOCITY,  PEAK  VELOCITY,  AND  PEAK 

274  C  DYNAMIC  PRESSURE  PROFILES  VERSUS  PROFILE  HEIGHT  (AGL) 

275  C  - 

276  C 

277  WRITE (IOU8, 38) 

278  88  format:  IX,  '  VARIABLES  ■=>  X,HT') 

279  C 

280 

281  81 

282 

283 

284  83 

285  82 


286 

287 

C 

WRITE :iOU8, 84) 

283 

289 

290 

84 

format:  IX, 'ZONE  T  -  "PEAK  PROFILE,  KTS",  1-xx,  F-POINT' 
DO  85  I  =  l.NPTS 

WRITE  :iOU8, 83)  VHPK :i)  ,  ZZ (I) 

291 

292 

C 

85 

CONTINUE 

2  93 

294  66 

295 

296 

297  87 
2  98  C 

299  C 

300  C  CLOSE  GRAPHICS  FILE 

301  C  - 

302  C 

303  CLOSE (IOU8, STATUS-'KEEP' ) 

304  C 

305  END IF 

306  C 

307  CALL  INKEY 

308  C 

309  GOTO  999 

31C  C 

311  C  - 

312  C  THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 

313  C  OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 

314  C  - 

315  C 


WRITE {IOU8, 86) 

format:  lx, 'zone  T  =  "PEAK  Q,  PSF”,  I=XX,  F='POINT' ) 

DO  87  I  =  1,NPTS 

WRITE  :iOU8, 83)  qhp:i),zz:i) 

CONTINUE 


WRITE  :iOU8, 81) 

format:  ix, 'ZONE  T  =  "MEAN  PROFILE,  KTS",  I=xx,  F=POINT' ) 

DO  82  I  =  1,NPTS 

WRITE ::ou8, 83)  vhmk :r) , zz :i) 

format:  ix, F6 . 1, IX, F6 .2) 

CONTINUE 
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316 

317 

318 

319 

320 

321 

322 

323 

324 

325 

326 

327 

328 

329 

330 

331 

332 

333 


2000  CONTINUE 

CALL  HOMCLS(O) 


C 

2001 


WRITE (lOUl, 2001) 
format (  /////8X, 

^  r  ***  error 

2  //,8X,' 


***  PLEASE  CHOOSE  A  NEW  OUTPUT  FILENAME' 
TYPE  <RETURN>  TO  CONTINUE  ' / $) 


READ  (lOUl,  '  (Al)  ■' )  TEMCHAR 


C 


KEY  -  'P' 


C 

999  CONTINUE 
C 

RETURN 

END 


C 
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1  C 

2  C 

3  SUBROUTINE  IREAD ( lOUl , PROMPT, IVALUE) 

4  C 

6  C  SUBROUTINE  IREAD  PROMPTS  USER  FOR  AN  INTEGER 

7  C  DATA  ENTRY  AND  CHECKS  VALIDITY  OF  ENTRY 

0  C  ********************************************** 

9  C 

10  PARAMETER (LAST-50) 

11  C 

12  CHARACTER*50  PROMPT, SHOWIT 

13  CHARACTER*15  ENTRY, BLANK 

14  C 

15  DATA  BLANK  /'  '/ 

16  C 

17  Q  ********************************************** 

18  C 

19  c  - 

20  C  PROMPT  USER  FOR  INTEGER  ENTRY.  FIND  POSITION 

21  C  OF  LAST  NON-BLANK  CHARACTER  IN  PROMPT, 

22  C  THEN  STORE  RIGHT  JUSTIFIED  IN  SHOWIT 

23  c  - 

24  C 

25  N  =  LAST  +  1 

26  C 

27  10  IF  (N.EQ.  DGOTO  2C 

28  C 

29  N  =  N  -  1 

30  C 

31  :F(PR0MPT(N:N) .EG.'  ')GOTO10 

32  C 

33  20  JS  =  LAST  -  N 

34  C 

35  WRITE (SHOWIT, ' (SOAl) ' )  ('  ' , J-1 , JS) , (PROMPT ( I : I ) , 1=1 , N) 

.36  C 

37  C  - 

38  C  NOW  ASK  USER  FOR  DATA  ENTRY 

39  C  - - 

40  C 

41  30  WRITE ( lOUl, ' (/, IX, A, 13) ' )  SHOWIT, IVALUE 

42  C 

43  WRITE (ICUl, '  (/, 8X, A, $)  '  ) 

44  1  '  ENTER  NEW  VALUE  OR  <RETURN>  TO  LEAVE  AS  IS  ==>  ' 

45  C 

46  READdOUl,  '  (A)  '  )  ENTRY 

47  C 

48  IF  (ENTRY. EG. BLANK) RETURN 

49  C 

50  READ (ENTRY,  '  (BN,  17)  ERB=3C)  ITEMP 

51  C 

52  IVALUE  =  ITEMP 

53  C 

.5  4  RETURN 

55  END 

56  C 
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1  C 

2  C 

3  FUNCTION  LEGAL (CHDOL, IOUl,OKLIST,NUM) 

4  C 

g  Q  *4rA4r*A*4rA*ft*ilr*4r*4f****AA'Ar»*i^********lbA*%A 

6  C  FUNCTION  LEGAL  DETERMINES  IF  THE  VALUE 

7  C  FOR  CHDOL  IS  A  VALID  INPUT.  THIS  VALUE 

8  C  IS  CHECKED  AGAINST  THE  LIST  OF  LEGAL 

9  C  VALUE  IN  ARRAY  OKLIST (NUM) 

Q  ****  h  **********  ±  *********  *************  * 

11  C 

12  CHARACTER*!  CHDOL, OKLIST (NUM) 

13  C 

^  Q  **************************************** 

15  C 

16  LEGAL  -  0 

17  C 

18  DO  10  1=1, NUM 

19  IF(CHDOL.EQ.OKLIST(I))  RETURN 

20  10  CONTINUE 

21  C 

22  LEGAL  ■=  1 

23  C 

24  WRITEdOUl,  '  (/,T9,A,A1,A/) ')  '  **•  CHDOL, 

25  1  '  IS  NOT  A  VALID  INPUT  ***' 

26  C 

27  RETURN 

2  8  END 

29  C 
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1 

2 

3 

A 

5 

6 
7 
ft 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 
61 
62 
63 


C 

C 

c 

c 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

r* 


c 


c 


c 

c 

c 


c 

c 

c 


c 


SUBROUTINE  LOCATE (IROW, ICOL) 

**ik«t*ARr4t*ikikA*lb*A*4r****A«******^*i^4rA**th’*#r*%*i^lkik**ft**#>r 

SUBROUTINE  LOCATE  LOCATES  THE  CURSOR  POSITION 

A**************************************************** 

COMMON  /  UNITS/  lOUl, IOU4, IOU5, IOU6, lOUV, IOU8, IGRAPH 

CHARACTER’S  CUP 
CHARACTER*!  ECUP(8) 

EQUIVALENCE  (CUP, ECUP (1) ) 

CHARACTER* 10  FMT 
CHARACTER*!  EFMT(IO) 

EQUIVALENCE  (FMT, EFMT (1) ) 

DATA  FMT  /  '  ("  ",  A?  \) '  / 

*^l^iiltl^l***********^t**^l»»»*****»**1l******^i^l^i^l^i1^t^l^i*^l^l** 


ANSI  CONTROL  SEQUENCE:  CUP  -  ESC [’ ROW' COLUMN' H 


IBl  ==  IROW/10 

IR2  =■  IROW  -  IRl’lO 

ICl  »  ICOL/10 

IC2  “  ICOL  -  IC1*10 

ECUP (1)  -  CHAR(27) 
ECUP (2)  -  CHAR(9l) 
IPOS  »  3 


IF(IR1.GT.0)THEN 

ECUP  (IPOS)  »  CHARdRl  +  48) 

IPOS  =  IPOS  +  1 

END  IF 

ECUP  (IPOS)  =■  CHAR(IR2  +  48) 

IPOS  “  IPOS  +  1 


ECUP (IPOS)  =  CHAR(59) 
IPOS  =  IPOS  +  1 


IF (ICl .GT. 0) THEN 

ECUP (IPOS)  =  CHAR (ICl  +  48) 

IPOS  =  IPOS  +  1 

END  IF 

ECUP(IPOS)  -  CHAB(IC2  +  48) 

IPOS  =  IPOS  +  1 


EC'JP(IPOS)  =  CHAR(72) 


EFMT (7)  =  CHAR (IPOS  +  48) 


WRITE  (lOUl,  *)  CUP 


RETURN 

END 
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1  C 

2  C 

3  SUBROUTINE  MOMENT (NPTS, HUMTYP, TOTF, TOTM) 

4  C 

5  Ik  **********  ir  ************************************  **4t  *  * 

6  C  SUBROUTINE  MOMENT  CALCULATES  THE  TOTAL  OVERTURNING 

1  C  FORCE  AND  MOMENT  ON  A  MAN  OR  YOUNG  PERSON  AND  PRINTS 

8  C  OUT  THE  RESULTS. 

9  ***************************************************** 

10  C 

11  CPAX  DIMENSION  PAXMAN(12) 

12  CHARACTER*!  HUMTYP 

13  C 

14  COMMON  /PERSON/  QP(12),DSET 

15  COMMON  /  UNITS/  lOUl , IOU4 , IOU5 , 10U6, IOU7 , IOU8 , IGRAPH 

16  C 

]_7  Q  ***************************************************** 

18  C 

19  C  - 

20  C  IN  THE  SUBROUTINE  TO  CALCULATE  FORCES  AND  MOMENTS: 

21  C 

22  C  "CDt"  IS  THE  COEFFICIENT  OF  DRAG  OF  THE  PERSON: 

23  C  CDP  -  1.1  FOR  HAZARD  ANALYSES 

24  C  CDP  -  1.0  FOR  CORRELATION  WITH  NATC  FLIGHT  TEST  DATA 

25  C  "WIDTHP"  IS  THE  WIDTH  OF  THE  PERSON  WHERE: 

26  C  WIDTHP  IS  'L'  TYPE  IF  1.1  FT 

27  C  WIDTHP  IS  'S'  TYPE  IF  0.7  FT 

28  C  - 

29  C 

30  WIDTH?  -  1.1 

31  CDP  =1.1 

32  C 

33  (3  - - - 

34  C  INITIALIZE  AREAS  FOR  NATC  MAN  FROM  ROTORWASH  FLIGHT  TES'^S 

35  C  (UNITS  ARE  FEET2) 

36  c  - 

37  C 

38  CPAX  CDP  -  1.1 

39  CPAX  PAXMAN(l)  -  0.41 

4  0  CPAX  PAX.MAN(2)  =0.37 

41  CPAX  PAXMANO)  -  0.365 

42  CPAX  FAXMAN(4)  -  0.42 

43  CPAX  PAXMAN(5)  -  0.5425 

44  CPAX  PAXMAN(6)  -  0.685 

45  CPAX  PAXMAN(7)  -  0.8 

46  CPAX  PAXMAN(8)  -  0.845 

47  CPAX  PAXMANO)  -  0.7875 

48  CPAX  PAXMAN(IO)  -  0.625 

49  CPAX  PAXMAN(ll)  -  0.3 

50  CPAX  PAXMAN(12)  »  0.00625 

51  C 

52  Q  - 

53  C  INITIALIZE  INTEGRATION  STEP  SIZE  AND  ZERO  SUMMATION  VARIABLES 

54  C  - 

55  C 

56  DELZ2  =0.5 

57  TOTM  =0.0 

58  TOTF  =0.0 

59  C 

60  C  - 

61  C  CHOOSE  HU14AN  SIZE  TYPE 

62  C  - 

63  C 
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5.4  IFCHUMTVP.EQ. 'S' >THEN 

tZ  WIDTH?  -  0.7 

66  NPTS  -  8 

67  END  IF 
69  C 

£9  c  - 

7  0  C  INTEGRATE  DWAMIC  PRESSURE  PROFILE  OVER 

71  C  THE  HEIGHT  OF  THE  PERSON  CHOSEN 

72  C  - 

73  C 

74  DO  10  J  -  1,NPTS 

75  C 

76  FOVER  -  QP ( J) *DELZ2*WIDTHP*CDP 

"7  CPAX  FOVER  -  QP ( J) *PAXMAN < J) *CDP 

78  ZZ  -  O.S-<J  -  1)  +  0.25 

79  OVERM  -  FOVER* ZZ 

80  TOTF  -  TOTF  +  FOVER 

81  TOTM  ■>  TOTM  +  OVERM 

82  C 

83  C  - 

84  C  PRINT  OUT  RESULTS 

65  C  - 

8  6  C 

87  i:'(DSET.NE.0.)GOTO  10 

80  C 

89  WRITE  (10U6, 20)  ZZ,QP  (J)  ,F07F.R,  OVERM,  TOTF,  TOTM 

90  20  FOi'MAT(  F3 . 2,  5  (2X,  FIO  .  3)  ) 

91  C 

92  1C  .  ■  fIHUE 

93  C 

94  RETUflN 

95  END 

96  C 

97 
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SUBROUTINE  PROPRM 


1  C 

2  C 

3  SUBROUTINE  PROPRM(H,UMB, RVZ) 

4  C 
b  C 

7  C  SUBROUTINE  PROPRM 

8  C 

9  C  THIS  SUBROUTINE  CALCULATES  THE  VELOCITY  PROFILE  V(R,Z) 

10  PARAMETERS 

11  C  OF  THE  RADIAL  WALL  JET  FOR  THE  NON-INTERACTING  ROTOR  CASE 

12  C 

14  C 

15  CO^LMON  /PROFIL/  RJ,  ZB J,  7-H J,  ZMJ,  UMJ,  ZB,  ZH ,  2M,  UM,  CU,  CY 

16  C 

17  C 

^0  ****************************************************************** 

19  C 

20  IF (RVZ .GE.RJ)GOTO  600 

21  C 

22  C  - 

23  C  RVZ  .LT.  RJ  — >  TRANSITION  REGION 

24  C 

25  C  TRANSITION  REGION  EQUATIONS,  EMPIRICALLY  APPLIED  BY 

26  C  JDK  PRIOR  TO  VERSION  2.1,  WERE  SIGNIFICANTLY  IMPROVED 

27  C  BY  THE  FOLLOWING  CHANGES  FOR  V2 . 1  IN  MAY  1992. 

28  C 

2S  C  OLD  OR  REPLACED  EQUATIONS: 

30  C 

31  C  UM  =  UMJ*RVZ 

32  C  IF (RVZ .GT. 1 . 0) 

33  C  ZH  =  2HJ*RVZ 

34  C  IF(RVZ.GT.l.O) 

35  C  ZM  =  ZMJ*RVZ 

36  C  IF (RVZ .GT . 1 . 0) 


.57 
38 
.3  9 
40 

C 

C 

c 

UM  ■•=  UMJ*  (RVZ/RJ)  **0 .5 

w 

c 

BOUNDARY  GROWTH  IN  TRANSITION 

REGION 

43 

c 

44 

c 

SEE  NOTE  ABOVE 

45 

46 

c 

OLD  OR  REPL.ACED  EQUATIONS; 

4  7 

c 

48 

c 

ZBO  -  1.5 

49 

c 

IF (H . LT , 1 . 5)  ZBO  -  H 

5  0 

ZH  =  (ZHO  -  ZHJ) /RJ**2* (RJ 

-  RVZ ) *  *  2  + 

ZHJ 

c  2 

c 

53 

ZBO  -  H**0.5 

54 

z;io  =■•  ZBO/2  ,  5 

^  *■. 

^  .) 

c 

5b 

ZH  =  (ZHO  -  ZHJ) /RJ* *2* (RJ  - 

RVZ) **1.5  + 

ZHJ 

57 

ZU  -  2,5*ZH 

b  ? 

ZM  ■•=  0.33‘ZH 

<-9 

{  ; 

.-,0 

GOTO  700 

0  i 

62 

c 

.  -  ~  w 

('3  C  RVZ  .GE.  RJ  '>'■=>  DEVELOPED  WALL  JET  REGION 


UM  -  UMJ 
ZH  -  ZHJ 
ZM  ■“  2MJ 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 


C 

C  SEVERAL  rOEFFICIENTS  IN  THE  GROWTH  EQUATIONS 
C  WERE  MOD  .TIED  IN  MAY  1992  FOLLOWING  THE 
C  CORRELAT.'.ON  EFFORT  (SEE  NOTES  IN  SUBROUTINE 

C  WALJET  FCr.  DETAILS)  . 

C 

C  OLD  EQUATIONS: 

C 

C  UM  -  CU*RVZ** (-1 . 143) *UMB 

C  ZH  -  CY^RVZ** (1.028) 

C  ZM  -  0.1944*ZH 

C  - 

c 

60 0  CONTINUE 
C 

UM  -  CU*RV2** (-1.0 *UMB 
ZH  -  CY*RVZ** (1 .0) 

ZB  -  2.8*ZH 
ZM  -  0.28*ZH 

C 

700  CONTINUE 
C 

RETURN 

END 

C 
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SUBROUTINE  VLINE 


1  C 

2  C 

3  SUBROUTINE  VLINE 

4  C 

5  C 

^  «*4r*A4kAtkAr**AAA***4r***^***4rikA*4rA*****'*‘Arik'**-***4r4rilr*4t*A4r**4r4tA**!kAAf;*** 

7  C  SUBROUTINE  VLINE 

8  C 

9  C  THIS  SUBROUTINE  APPLIES  THE  BIOT-SAVORT  LAW  TO 

10  C  CALCULATE  THE  VELOCITY  INDUCED  BY  A  LINE  VORTEX 

11  C 

12  C  XA,YA, ZA  -  STARTING  POINT  OF  VORTEX 

13  C  XB,YB,  ZB  -  ENT-ING  POINT,  OR  DIRECTION  POINTER 

14  0  XC,YC,ZC  -  TARGET  POINT  WHERE  VELOCITY  IS  INDUCED 

15  C 

16  C  IFI  ■  0  VORTEX  IS  FINITE,  FROM  POINT  A  TO  POINT  B 

17  C  IFI  -  1  VORTEX  IS  SEMI-INFINITE  FROM  POINT  A  THROUGH  B 

18  0 

19  C 

20  ********************  V,  ************* ******************************** 

21  0 

22  COMMON  /CONSTS/  PI,RHO,FPSPKN,RHOD2,DRC 

23  COMMON  /CVLINE/  IFI , XA, YA,  ZA, XB, YB,  ZB, XC, YC, ZC, Ql, Q2 , Q3 

24  C 

25  C 

25  *«**-*ir**2(*4r***4r*4r'*r*^7^*^(k-^A*Ar*i#A***A*****4r******«1k********ir'****ilr*** 

27  C 

28  A  -  (XA-XC)**2  +  (YA-YC)**2  +  (ZA-ZC) **2 

29  B  -  2.0*(  (XA-XB) * (XC-XA)  +  (YA-YB) » (YC-YA>  +  ( ZA-ZB) * ( ZC-ZA) 

30  ) 

31  C  -  (XA-XB) **2  t  (YA-YB) **2  +  (ZA-ZB) **2 

32  C 

33  Cl  -  (YC-YB)*ZA  +  (YA-YC) *ZB  +  (YB-YA) *ZC 

34  C2  -  (ZC-ZB)*XA  +  (2A-ZC)*XB  +  (ZB-2A) *XC 

35  C3  -  (XC-XB)*YA  4  (XA-XC) *YB  +  (XB-XA) *YC 

36  C 

37  Q  -  4.0*A*C  -  B**2 

38  C 

39  c  - 

40  C  CHECK  FOB  COLINEAR  TARGET  POINT 

41  C  - 

42  C 

43  QB  -  0.0 

44  IF(ABS(Q) .LT.1.0E-06)GOTO  100 

45  C 

46  C  - 

47  C  FINITE  LENGTH  VORTEX 

48  C  - 

49  C 

50  IF(IFI.EQ.0)THEN 

51  QB  -  1.0/Q* ( (2.0*0  4  B)/SQRT(A  4  B  4  C)  -  B/SQRT (A) ) /2 . 0/PI 

52  ENDIF 
=^3  C 

54  C  - 

55  C  SEMI-INFINITE  VORTEX 

56  C  - 

57  C 

58  IFdFI.EQ.DTHEN 

59  QB  “  1 .0/Q* (2 .0*SQRT (C)  -  B/SQRT (A) ) /2 . 0/PI 

60  ENDIF 

61  C 

62  100  CONTINUE 

63  C 
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SUBROUTINE  VLINE 


S4  C  - 

65  C  VELOCITY  COMPONENTS 

6  6  C  - - 

67  C 

68  Q1  -  C1*QB 

69  Q2  -  02*08 

70  03  -  C3*0B 

71  C 

72  RETURN 

7  3  END 

74  C 
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SUBROUTINE  WALJET 


1  C 

2  C 

3  SUBROUTINE  WAWET  (H,UB,UN,UMB) 

4  C 

5  C 

g  ir*4r4r***4tAi**4r*t^*A«*ik4rA**i<k*AA***i^****Atfi#****##A*****A*A#tiiir*A**i(^*'ArA#A 

7  C  SUBROUTINE  WALJET 

8  C 

9  C  THIS  SUBROUTINE  CALCULATES  THE  STARTING  POSITION  OF  THE 

10  C  WALL  JET  AND  GROWTH  PARAMETERS  FOR  WALL  JET  DECAY 

11  C 

12  ******  ************ii********************<t***t*****#****w*******<i*** 

13  C 

14  COMMON  /CLOUDK/  QSMAX 

15  COMMON  /PROFIL/  RJ, ZBJ, 2HJ, ZMJ, UMJ, ZB, ZH, ZM, UM, CU, CY 

16  COMMON  /  UNITS/  lOUl, IOU4, IOU5, IOU6, IOU7 , lOUS, IGRAPH 

17  C 

18  C 

*****tfr*-***ir*****4r**i^4i****************'**********ft4r**A****1i^**A***** 

20  C 

21  c  — --- — - - - - - - - 

22  C  INITIALIZATION  OF  EXPONENTS 

23  C 

24  C  CZM,  THE  PROFILE  PEAK  VELOCITY  LOCATION,  WAS  INCREASED 

25  C  FROM  THE  GLAUERT  VALUE  (0.1944)  TO  THE  PRESENT  VALUE 

26  C  BASED  ON  DATA  AND  ASSUMPTIONS  PRESENTED  IN  FAA  REPORT 

27  C  FOR  VERSION  2.1,  MAY  1992.  THE  OTHER  COEFFICIENTS  WERE 

28  C  MODIFIED  WHEN  IT  WAS  DEMONSTRATED  THAT  CORRELATION  WAS 

29  C  IMPROVED  WITH  BOTH  MODEL  AND  FLIGHT  TEST  DATA. 

30  C 

31  C  VARIABLE  VALUE  FOR  V2 . 1  BEFORE  V2 . 1 

32  C 

33  C  EXU  -1.0  -1.143 

34  C  EXY  1.0  1.028 

35  C  CZM  0.20  0.1944 

36  C  QSMAXO  1.077  1.0  (SEE  NOTE -BELOW) 

37  C  - - 

30  C 

39  EXU  “-1.0 

40  EXY  -  1.0 

41  EXM  -  1.0+  EXU  +  EXY 

42  CZM  -  0.28 

43  CZB  -  2.8 

44  C 

45  C  - 

46  C  ITERATE  TO  FIND  INITIAL  RADIUS  OF  WALL  JET,  RJ 

47  C  - 

48  C 

49  TOL  «  l.OE-05 

50  RJ  -  2.0 

51  C 

52  QSMAXO  -  1 , 077 

53  C 

54  DO  100  1-1,20 

55  C 

56  C  - 

57  C  CALCULATION  OF  EQUIVALENT  JET  LENGTH 

58  C  - - 

59  C 

60  TK  -  H  +  (RJ  -  1.0) 

61  IDE  -  0.707*TR 

62  C 

63  C  - - - . . 
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64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 

109 

110 
111 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 


C  QSMAX  CURVEriT  ORIGINALLY  TO  FIG.  8,  USAAVLABS  TECHNICAL 

C  REPORT  68-52,  JULY  1968.  UPDATED  INFORMATION  FROM  FIG.  2, 

C  REPORT  DTNSRDC/ASED-79/04,  APRIL  1979. 

C 

C  OLD  TR  68-52  EQUATIONS: 

C 

C  QSMAXO  -  1.0 

C  IF{TDE.LE.4.0)  QSMAX  -  QSMAXO  +  (0.6  -  QSMAXO) /16 . 0*TDE* *2 

C  IF(TDE.GT.4.0)  QSMAX  -  2.4/TDE 

C 

C  RJNEW  COEFFICIENTS  ROUNDED  OFF  DURING  MAY  1992  CORRELATION 

C  EFFORT  WHICH  PRODUCED  IMPROVED  RESULTS  WHEN  SIMPLIFICATIONS 

C  WERE  INTRODUCED 

C 

C  OLD  EQUATION: 

C 

C  RJNEW  -  2.508078*(UB/UM) **(0.486) 

C  - 

c 

IF(TDE.LE.3.5)  QSMAX  -  1.08  -  0.025*TDE**2 
IF(TDE,GT.3.5)  QSMAX  -  2.7/TDE 
C 

UM  -  SORT (QSMAX) 

C 

RJNEW  -  2.5* (UB/UM) ** (0.5) 

C 

IF  (ABS  (RJNEW  -  RJ)  .LE.TODGOTO  200 
RJ  -  RJNEW 
C 

100  CONTINUE 
C 

WRITEdOUl,  10) 

10  FORMAT(  '******»*****»*************•*»****»*************',/ 

1  ITERATIONS  EXCEEDED  FOR  WALL  JET  INITIAL  RADIUS',/ 

2  ^  »**********************•****•*******************' ) 

C 

STOP  '  ' 

c 

200  CONTINUE 
C 

RJ  “  RJNEW 
C 

C  - 

C  VELOCITY  GROWTH  FUNCTION  CONSTANTS 

C 

C  TWO  CONSTANTS  WERE  ROUNDED  OFF  TO  SIMPLIFY  THE  FOLLOWING 

C  TWO  EQUATIONS  DURING  THE  CORRELATION  EFFORT  OF  MAY  1992 

C  FOR  VERSION  2 . 1 

C 

C  OLD  EQUATIONS: 

C 

C  UMB  -  ( (0. 3586*RJ**EXM* (UM*UN) * (UB*UN) ** (0 . 14) ) **  (0. 88) ) /UN 

C  ZHJ  -  0 . 654/ (UM/UMB) **2/RJ 

C  - 

Q 

UMB  -  ( (0 .36*RJ**EXM* (UM*UN) * (UB*UN)  **  (0 . 14) )  »• (0. 88) ) /UN 

ZHJ  -  0. 65/ (UM/UMB) •*2/RJ 
CU  -  UM/UMB*RJ** (-EXU) 

CY  -  ZHJ*RJ* • (-EXV) 

C 

C  - 

C  MAX  VELOCITY  AND  BOUNDARY  PARAMETERS  AT  RJ 

C  - 
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127  C 

128 

129 

130 

131 

132  C 

133 

134 

135  C 


UMJ  -  CU*RJ**(BXO) *UMB 
ZHJr  -  CY*RJ**(EXY) 

ZMJ  “  CZM*ZHJ 
ZBJ  “  CZB*ZHJ 

RETURN 

END 
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1 

2 

3 

i\ 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 
2  5 
26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


45 

46 

47 

48 
4  9 

50 

51 

52 

53 


C 

C 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c; 

c 

c 


SUBROUTINE  WJVEL (H, UN, UMB,  RVZ,  RADIUS, WSPD,  DELZ, ZMAX,  DXO, BDLAVM) 

SUBROUTINE  WJVEL  GENERATES  THE  VELOCITV  PROFILE  V(R,Z) 

AT  RVZ  FOR  THE  NON-INTERACTING  ROTOR  CASE 

CHARACTER*!  TEMCHAR 
CHARACTER*!  ICONT(5) 

CHARACTER*!  KEY, KKEY 
CHARACTEFv*l2  PTSFIL(4) 

CHARACTER*50  COMM (2) 

DIMENSION  ZZ(60) , VMF(60) , VMK ( 60) , VPF ( 60 ) , 

1  VPK(60) ,QM(60) ,QP(60) 

COMMON  /  CKEY/  KEY,KKEy 

COMMON  /CONSTS/  PI,  RHO,  Fl’SPKN,  RHOD2,  DRC 

COMMON  /INPUTC/  ICCNT, COMM, PTSFIL 

COMMON  /PROFIL/  RJ, ZBJ, ZHJ, 2MJ, UMJ, ZB, ZH, ZM, UM, CU, CY 
COMMON  /  UNITS/  lOUl , IOU4 , IOU5, IOU6, IOU7 , IOU8 , IGRAPH 


^A**AiA**AA*AA**i 


t*«*******«*****«***A*A***i 


'PROP.'IM'  PROVIDES  THE  VELOCITY  PROFILE  PARAMETERS 
OF  A  RADIAL  WALL  JET  (WITHOUT  INTERJECTION  PLANE) 


CALL  PROPRM(H,UMB,RVZ) 


DIMENSIONALIZE  VELOCITY  PROFILE  PARAMETERS 


38 

Ri'VZ 

RADIUS'RVZ 

,•9 

ZZb 

a 

ZB'RADIUS 

40 

ZZH 

V 

ZH*RADIUS 

41 

4,?  C 

ZZM 

ZM*RADIL'S 

43 

2ETAH 

ZH/ZB 

44 

ZETAM 

ZM/ZB 

OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILE  HEADER 


ICD  -  0 

CALL  HOMCLS(ICD) 

IF(IOU6.NE.IOU!)  WKvrE(IOU6,  '  ("  !' 


54  IFdOUe.EQ.  6)  WRITE  (1006,  93)  COMM  ( 1 )  ,  COMM  (2  ) 

55  93  FORMAT(  1  OX, A50 , / ,  1  OX, A50 , / / ) 

56  C 

57  RVZOUT  ■  RRVZ  t-  DXO 

58  C 

59  WRITE (10U6, 1000)  RVZOUT 

6C  1  000  FORMAT(  9X,  'SINGLE  ROTOR  VELOCITY  PROFILE  AT  KADIUO  ■■ 

61  1  F7,l,'  FT  ,/) 

62  C 

63  WRITE (IOU6, 1001)  ZZU, ZZH, ZZM 
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64 

1001  FORMAT (  15X, 'PROFILE 

BOUNDARY 

HEIGHT  - 

',F7.2,' 

FT' 

,/ 

65 

1  #15X,' 

HALF-VEL. 

HEIGHT  - 

FT' 

J 

66 

2  ,15X,' 

MAX-VBL 

HEIGHT  - 

',F7.2,' 

FT' 

./) 

67  C 

68  C  - 

69  C  INCREMENTS  AND  HEIGHT  ARE  FROM  DELZ  AND  ZMAX 

70  C  - 

71  C 

72  NPTS  -  IFIX(ZMAX/DELZ)  +  1 

73  C 

74  C  - 

75  C  BOUNDARY  LAYER  REGION  EXPONENT 

76  C  'AN'  IS  ACTUALLY  -  1. 0/7.0 

77  C  - 

7  8  C 

79  AN  -  0.142857142 

80  C 

81  C  - 

82  C  SHEAR  LAYER  REGION  EXPONENT,  TO  MEET  EDGE  CONDITIONS 

83  C  (FROM  FIGURE  7,  USAAVLABS  TECHNICAL  REPORT  68-52,  JULY  1968) 

04  c  - 

85  C 

86  ALPW  -  ALOGd.O  -  1 . 0/SQRT  (2 . 0)  ) /ALOG  (  (ZH  -  ZM)  /  (ZB  -  ZM)  ) 

07  C 

88  VN  -  UN 

89  VMN  -  UM 

90  C 

91  C  - 

92  C  CALCULATION  OF  THE  NON-DIMENSIONALIZED  MINIMUM  ALLOWED 

93  C  BOUNDARY  LAYER  THICKNESS  SO  THAT  THE  BOUNDARY  LAYER  CAN 

94  C  BE  ADJUSTED  IF  THE  ZM  POSITION  IS  PHYSICALLY  TOO  LOW 

95  C  (BDLAYM,  IN  FEET,  COMES  FROM  A  MENU  INPUT  PARAMETER) 

96  c  - 

97  C 

98  ZETAl  -  BDLAYM/ ZB 

99  C 

100  C  - 

101  C  OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILE  HEADER 

102  C  - 

103  C 

104  WRITE (IOU6,  1005) 

105  1005  FORMAT (  2X, ' HEIGHT' , 5X, ' MEAN  VELOCITY' , 7X, ' PEAK  VELOCITY' , 6X, 

106  1  'MEAN  Q' ,4X,' PEAK  Q',/, 

107  2  3X, ' (FT) ' , 5X, ' (FPS) ' , 6X, ' (KN) ' , 5X, ' (FPS) ' , 6X, ' (KN) ' , 5X, 

108  3  '  (PSF) ' , 5X, '  (PSF) ' , /) 

109  C 

110  C  - 

111  C  CALCULATE  THE  VELOCITY  PROFILE  POINTS  FOR  OUTPUT 

112  C  - - - - 

113  C 

114  LINES  -  0 

115  C 

116  DO  500  I  -  1,NPTS 

117  C 

110  LINES  -  LINES  +  1 

119  2  -  DELZ*FLOAT(I  -  1) 

120  2ETA  -  Z/ZB 

121  C 

122  I F ( ZETA . LT . ZETAM  . OR . ZETA . LT . ZETAl ) THEN 

123  C 

124  C  - - - - 

125  C  Z  IS  WITHIN  BOUNDARY  LAYER 

126  C 
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127  C 

128  C 

129  C 

130  C 

131  C 

132  C 

133 

134  C 

135 

136  C 

137 

138  C 

139 

140  C 

141 

142 

143  C 

144 

145  C 

146 

147  C 

148 

149  C 
15C 

151  C 

152 

1  Q  T 

X  <»'  vO  w 

154 

155  C 
15  6  C 

157  C 

158  C 

159  C 
16C  C 

161  C 

162  C 

163  C 

164  C 

165  C 

166  C 

167  C 

168  C 

169  C 

170 

17  1  c 

172 

173  C 

174 

175  C 

176 

nv  C 
178 

17  5 

18  0  C 
181 

18  2  C 
183 

19  4  C 
18  5 

136  C 
181 

18  8  C 
199 


NOTE  THAT  THE  BOUNDARY  LAYER  CALCULATIONS  NOW  USE 
THE  MINIMUM  THICKNESS  PARAMETER  AND  THE  PEAK  TO 
MEAN  VELOCITY  PARAMETER  IS  THE  MAXIMUM  VELOCITY 
HEIGHT  RATIO  (AT  ZM) .  ADDED  MAY  1992  FOR  V2 . 1 . 


V2M  =0.0 

IF (2ETAM.GT. 0 .0)  THEN 

VZM  =  (ZETA/ZETAM) **AN 
IF (ZETAl . GT . ZETAM) THEN 

VZMl  =  (1.0  -  ((ZETAl  -  2ETAM)/(1.U  -  ZETAM) ) * *ALPW) * *2 
VZ,-.  -  (Zi^IA/ZETAl)  **AN 

ENDIF 

VMTOPK  =  1.04653  +  0.373894*RVZ  -  0 . 0422525*KVZ*RVZ 
It  (VMTOPK.LT.1.2)  VMTOPK  =  1.2 
ENDIF 
GOTO  400 
ENDIF 


Z  IS  WITHIN  SHEAR  LAYER 

THE  PEAK  TO  MEAN  VELOCITY  RATIO  EQUATIONS  ARE 
SUBST^iJTlALLV  IMPROVED  OVER  THOSE  USED  PRIOR  TO 
MAY  1992.  EQUATIONS  ARE  NOW  USED  FOR  BOTH  THE 
MAXIMUM  VELOCITY  HEIGHT  (ZM)  AND  THE  1/2  VELOCITY 
HEIGHT  (ZH) .  VALUES  BETWEEN  ARE  INTERPOLATED  AND 
VALUES  ABOVE  ZH  USE  THE  ZH  RAT  10* ( ZETA/ ZETAH )  . 
THESE  2nd  ORDER  EQUATION  S'JBSTANTI ALLY  IMPROVED 
CORRELATION  WITH  .MODEL  AND  FLIGHT  TEST  DATA 
DURING  THE  MAY  1992  EFFORT  FOR  V2 . I . 


VZM  =0.0 

IF (Z . LE . ZB) THEN 

VZM  =  (1.0  -  ((ZETA  -  ZETAM)/ (1.0  -  ZETAM )  )  *  *  ALP W)  *  *  2 

IF ( ZETA . GE . ZETAH) THEN 

VMTOPK  -  (  1.48086  +  0.569177*RVZ  -  0 . 0 6 925 1  4  * RVZ * BVZ ) 
* (ZETA/ZETAH) 

rF(VMTOPK.  LT.  1  ,2)  \T1TOPK  =  1.2 
ELSE 

VMPKMX  =  1.C4653  +  0.373894  *RVZ  -  0 . 0 4 22 52 5  * RVZ * RVZ 
VMPK12  =  1.48086  +  0.565177*RVZ  -  "  25i4*RVZ*3VZ 
FRAC  =  (ZETA  -  ZETAM) / (ZETAH  -  ZETAM) 
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190  C 

191 

192  C 

193 

194  C 

195 

196  C 

197 

198  C 

199 

200  C 

201 

202  C 

203 

204  C 

205  400 

206  C 

207 

208  C 

209  C 

210  C 

211  C 

212  C 

213 

214  C 

215  C 

216  C 

217  C 

218  C 

219 

220 

221  C 

222  C 

223  C; 

224  C 

225  C 

226 

227 

228  C 

229 

230  C 

231  C 

232  C 

233  C 

234  C 

235  C 

236  C 

237 

2  38  C 

239 

240  C 

241 

242 

243 

244 

245 

246  C 

247  55 

248  C 

249  C 

250  C 

251  C 

252  C 


ir  (ZETAl  .GT  .  ZETWl)  THEN 

FRAC  -  (ZETA  -  ZETAl) / (ZETAH  -  ZETAl) 
ENDIF 

VMTOPK  -  FRAC*VMPK12  +  (1.0  -  FRAC) *VMP KMX 
IF (VMTOPK.lt. 1.2)  VMTOPK  -  1.2 
ENDIF 
ENDIF 
CONTINUE 
VZN  -  VZM*VMN 

DIMENSIONAL  HEIGHT 

2Z(I)  »  Z*RADIUS 

MEAN  VELOCITIES 

VMF(I)  -  VZN*VN 
VMK(I)  -  VMF(I) /FPSPKN 

PEAK  VELOCITIES 

VPFd)  -  VMF  (I) ‘VMTOPK 
VPK(I)  -  VPFd) /FPSPKN 

IF(VPK(I) .EQ.0.0)GOTO  55 


THE  EFFECT  OF  WIND  IS  TO  ADD  (DOWNWIND  SIDE)  OR  SUBXBACT 
(UPWIND  SIDE)  'XWK'  TIMES  THE  AMBIENT  WIND  VELOCITY  TO 
THE  HORIZONTAL  PROFILE  VELOCITY  (EMPIRICAL,  CH-53E  BASED) 


XKW  -  (-0.5*H)  +  2.5 

IF(XKW.LT.l.O)  XKW  -  1.0 

WSPD2  -  WSPD*XKW 
VMKd)  -  VMKd)  +  WSPD2 
VMFd)  •«  VMKd)*FPSPKN 
VPKd)  •  VPK(I)  +  WSPD2 
VPFd)  -  VPXd)‘FPSPKN 

CONTINUE 


DYNAMIC  PRESSURE 


E-C3 


SUBKOUTINE  WJVEL 


253  QM(I)  -  RHOD2*VM£'(I)  **2 

254  QP(I)  -  RHOD2*VPF(I) **2 

255  C 

256  IF(IOU6.EQ.IOUl)THEN 

257  IF (LINES .LT. 12) GOTO  450 

258  LINES  "•  1 

259  CALL  INKEy 

260  IF (KEY.NE. 'C' )GOTO  999 

261  WRITE (10U6,  1005) 

262  ENDIF 

263  C 

264  450  CONTINUE 

265  C 

266  C  - - - - - 

267  C  OUTPUT  THE  VELOCITY  AND  DYNAMIC  PRESSURE  PROFILES 

268  C  - 

269  C 

270  WRITE (IOU6,  1002)  ZZ ( I)  ,  VMF (I) , VMK ( I) ,  VPF  ( I ) , 

271  *  VPK(I) ,QM(I) ,QP  (I) 

272  1002  FORMAT(  F8.2,6F10.3) 

273  C 

274  500  CONTINUE 

275  C 

276  C  - 

277  C  WRITE  OUT  GRAPHICS  FILES  IF  SWITCH  IS  SET  BY  USER 

278  C  - 

279  C 

280  IF  (IGRAPH.EQ.  DTHEN 

281  C 

282  C  - - - 

283  C  OPEN  GRAPHICS  FILE 

284  C  - - - 

285  C 

286  OPEN(IOU8,FILE“PTSFIL(l) , STATUS”' NEW'  , ERR=2000) 

287  C 

288  WRITE (IOU8, 89)  COMM ( 1) , COMM (2) 

289  89  FORMA'K  1  OX, A50 , / , lOX, A50, // ) 

290  C 

291  WRITE(IOU8, 80)  RVZOUT 

292  80  FORMAT(  IX,  '  TITLE”’'VELOCITY  PROFILE,  DFRC  =',F5.1,'  FT,' 

293  *  '  GW  =  xxxxx  LB,  WAGL  =  xx . x  FT"') 

294  C 

295  C  - 

296  C  PRINT  OUT  MEAN  VELOCITY,  PEAK  VELOCITY,  AND  PEAK 

297  C  DYNAMIC  PRESSURE  PROFILES  VERSUS  PROFILE  HEIGHT  (AGL) 

298  C  - 

299  C 

300  WRITE (IOU8, 88) 

301  88  FORMAT (  IX, ' VARIABLES  =  X,HT') 

302  C 

303  WRITE (IOU8, 81) 

304  81  FORMAT  (  IX,  ' ZONE  T  ”  "MEAN  PROFILE,  KTS",  I  =  xx,  F=PGINT'  ) 

305  C 

306  D082I=  1,NPTS 

307  WRITE (IOU8, 83)  VMK(I),ZZ(I) 

308  83  FORMAT)  IX, F6 . 1 , IX, F6 . 2 ) 

309  32  CONTINUE 

310  C 

311  WRITE (IOU8, 84) 

312  84  FORMAT)  IX, 'ZONE  T  =  "PEAK  PROFILE,  KTS",  I-xx,  F=POINT' ) 

313  C 

314  DC  85  I  =  1,NPTS 

315  WRITE (IOU8, 83)  VPK(I),ZZ(1) 


E-84 
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316 

85 

317 

318 

C 

319 

86 

320 

321 

322 

C 

323 

87 

324 

C 

325 

c 

326 

c 

327 

c 

328 

329 

c 

330 

331 

c 

332 

333 

c 

334 

335 

c 

336 

c 

337 

c 

338 

c 

339 

c 

340 

c 

341 

c 

342 

2CC0 

343 

344 

345 

f 

346 

347 

348 

349 

350 

2001 

351 

r 

352 

999 

353 
35  4 
355 

c 

356 

357 

v-, 

CONTINUE 
WRITE (IOU8, 86) 

FORMAT (  IX, 'ZONE  T  -  “PEAK  Q,  PSF",  l=xx,  F=POINT' ) 

DO  87  I  -  1,NPTS 

WRITE  (IOU8, 83)  QP(1),ZZ(I) 

CONTINUE 


CLOSE  GRAPHICS  FILE 

CLOSE (lOUa, STATUS='KEEP'  ) 
ENDIF 

CALL  INKEY 
GOTO  999 


THE  ERROR  LOGIC  ALLOWS  FOR  THE  HANDLING  OF  FILE 
OPEN  ERRORS  BY  RETURNING  THE  USER  TO  A  MENU 


CONTINUE 

CALL  HOMCLS(O) 

WRITEdCUi,  2001) 

FORMAT (  ////,8X, 

'  ***  ERROR  *‘*  PLEASE  CHOOSE  A  NEW  OUTPUT  FILENAME' 

//,8X,'  TYPE  <RETURN>  TO  CONTINUE  ' , S) 

READ (lOUl, ' (Al) ' )  TEMCHAR 
KEY  =  ' P' 

CONTINUE 


RETURN 

END 


